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PREFACE 


The objective of this study was to investigate the 
feasibility of implementing an on-board aircraft navigation 
system based on using the ATS-I and ATS-III satellites, the 
Omega Position Location Equipment (OPLE) Control Center 
(NASA/GSFC Contract NAS 5-10248) as modified for an aircraft 
surveillance function in the VHF Satellite Navigation Experiment 
(NASA/GSFC Contract NAS 5-21079) , and a suitable aircraft 
terminal. 

The study was performed by first considering the problems 
associated with satellite navigation and how the specified 
components could be employed effectively in that context. From 
this effort the optimum system was defined. Next, the modifi- 
cations of the OPLE Control Center were described, which were 
necessary to configure a low-cost experimental system. In 
addition, an aircraft terminal and stationary ground reference 
terminals were designed. Although the resulting system was not 
optimum, the expected performance was sufficient to demonstrate 
at relatively low cost the principles developed in the study. 

The final task was to outline an experiment plan designed to 
demonstrate the navigational performance of the system, and an 
estimate of the cost of conducting the program was established. 

It was concluded that the on-board navigation system was 
feasible. With relatively minor OPLE Control Center modifica- 
tions it is expected that the system would provide aircraft 
navigation accuracy to approximately 1600 meters in a controlled 
experiment. With additional modifications the system would 
provide approximately 1500-meter accuracy over the continental 
United States in a configuration approximating an operational 
system. 
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SECTION I 
INTRODUCTION 

Texas Instruments was selected by NASA Goddard Space Flight 
Center to conduct a study of the feasibility of implementing an 
on-board navigation system based on using both ATS-I and ATS-III 
satellites, the Omega Position Location Equipment (OPLE) Control 
Center as modified for the VHF Satellite Navigation Experiment, 
an aircraft terminal, and other existing navigational aids. 

This report documents the results of the study. 

It was concluded that an on-board satellite navigation 
system comprised of the specified components is feasible. With 
relatively minor OPLE Control Center modifications such a system 
will provide aircraft navigation accuracy to approximately 1600 
meters in a controlled experiment. With additional modification 
the system will provide approximately 1500-meter accuracy over 
the continental United States in a configuration approximating 
an operational system. 

^ The content of this report is divided into three main 

sections. The first of these sections, which is the study 
report, considers the problem of satellite navigation and ob- 
jectively defines the optimum system. Constraints imposed by 
* the ATS-I, ATS-III, and the modified OPLE Control Center pre- 

clude implementation of the optimum system; however, the study 
clearly delineates the problems and error sources that arise in 
implementing an experimental system using recommended minor 
modification to the OPLE Control Center. 
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The second section describes modifications to the existing 
OPLE Control Center, the design of the aircraft terminal, and 
the design of ground reference terminals which would be used in 
an experiment intended to demonstrate the conclusions reached 
in the study. The final section outlines an experiment plan 
to effect the demonstration, and provides an estimation of the 
costs to be expected in conducting the program. 



SECTION II 


NAVIGATION SYSTEM ANALYSIS 


1. REQUIREMENTS FOR A NAVIGATION SYSTEM 

The ideal satellite navigation system will be accurate, 
have a wide coverage area, and be economically accessible to a 
large number of users. The last constraint is the most diffi- 
cult to achieve technologically, as illustrated by the fact 
that the bare minimum of user equipment required in a system 
such as the one developed in this report is a sophisticated VHF 
receiver system with a high stability clock synchronized to a 
ground control center (GCC) . The receiver and clock are 
necessary to determine the satellite (s) -to-aircraft ranges. 

Since the receiver and clock must be carried by all users, the 
dream of providing an inexpensive satellite navigation service 
to the private pilot is still unrealized. For commercial users, 
however, a satellite navigation system is potentially cheaper 
than, and of comparable accuracy to, existing inertial 
navigation systems. 

Wide coverage is essential to a satellite navigation system 
to make it economical. This follows from the fact that a major 
equipment cost is the cost of the ground control equipment, an 
expense that will most likely be shared among the users; there- 
fore, more users result in a lower cost per aircraft. Clearly, 
the maximum usage is obtained when the navigation service is 
provided in all areas frequented by commercial air traffic. 

Navigation accuracy is a function of equipment design. 
Categories of errors are thoroughly discussed and analyzed in 
this report. The results of the error study which are pre- 
sented in Section II-3, along with cost considerations, dictate 
the design of the recommended system. In addition to the above 
listed considerations, desirable features of a satellite 
navigation system are: ease of maintenance, reliability, and 
aircraft in-flight independence from ground equipment. This 
last feature refers to the need of an aircraft to remain in time 
synchronization with the ground control station without the 
necessity of ground interrogation of the aircraft. Texas 
Instruments believes that the system recommended in the 
following pages meets these requirements. 

1.1 Alternative Concepts 

There are 16 general configuration classifications for 
satellite navigation systems. A tree of' the various combina- 
tions is shown in Figure 2-1. In the following paragraphs the 
trade-offs involved in making each choice are briefly described 
and analyzed. 
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Figure 2-1. Tree of Possible Satellite .Navigation Configurations 
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1.1.1 Two vs. Three Reference Terminal Configuration 

Aside from the constraint that only two satellites will 
be available for the Dual Satellite Navigation Experiment, thus 
requiring that a third ranging tone source be ground based, the 
three reference terminal concept has been extensively studied 
and tested in the OPLE program in which Texas Instruments was a 
major contributor [1,2]. The technical advantage of the three 
reference terminal concept is that user position may be 
determined without allusion to altimeter data, whereas the two 
reference terminal configuration requires an external input, 
such as aircraft altitude data, to obtain a position location 
for the aircraft. Placement of a third reference station on the 
ground, while enabling implementation of the three reference 
terminal navigation concept, is impractical because all three 
references are simultaneously in view only over a small sector 
of the earth. In fact, the area of coverage is limited to a 
circle with a radius of no more than approximately 100 km for 
the 50°W and 150 °W longitude placement of the ATS-I and ATS-III 
satellites. For this reason range measurements will be obtained 
from just two satellite references and combined with on-board 
aircraft altitude data to determine position from one of the two 
points of intersection of the three spheres centered at the two 
satellites and the center of the earth. The second intersection 
will be in the southern hemisphere - sufficiently far away for 
the derived user position to be non-ambiguous . 

1.1.2 Ranging Tone Source Trade-Offs 

Ranging tones are the tones emitted from each reference 
at a constant and known phase which are used to obtain range 
information. They may be generated on the ground and transponded 
to users by the satellites or they may be generated at the satel- 
lites themselves. The advantage of placement of the source of 
ranging signals in the satellite is that the system becomes very 
difficult to jam because the satellite becomes an autonomous 
source, and therefore intentional jammers are forced to jam 
individual users. Anti- jamming features, however, are not a 
major advantage to non-military navigation systems, and con- 
sequently the additional cost and lower reliability are not 
warranted. It is, therefore, recommended that ranging tones be 
generated at a ground control center and relayed to the aircraft 
users via satellite 

1.1.3 Real-time Ephemeris Data vs. Pre-programmed Ephemeris 

Data 

The ATS-I and ATS-III synchronous satellite reference 
stations are not adequately stable in position to be used as 
absolute reference points. The earth-referenced position of 
each satellite varies in a daily cycle because of the orbital 



inclinations of each satellite with respect to the axis of earth 

rotation. This causes the subsatellite point, which is that ' 

point where the radius vector joining the center of the earth 

and the satellite intersects the earth surface, to trace out 

a figure "eight" during the course of one complete rotation of 

the earth. This movement may be predicted and compiled in some * 

form as ephemeris data. Two data formats are available for use 

in a satellite navigation system: (1) listing such as available 

in R204 Operational Predicts and World Map Routine , published 

by NASA, and T21 the simplified set of equations such as derived 

in Section II-3.1.1.2 of this study report. 

Ephemeris data is necessary for accurate position 
location because it provides the user with satellite position 
information of sufficient accuracy to utilize the range data 
derived from the GCC-transmitted ranging tones. The list 
format implies a need for a real-time ephemeris update via 
data- link through the satellite to the aircraft because the on- 
board data processor envisioned for the aircraft terminal will 
not have sufficient data storage capability to store the number 
of satellite coordinates compatible with the aircraft position 
location accuracy constraint and the aircraft around inde- 
pendence constraint. On the other hand, the simplified equation 
format does not require a large data storage facility. As shown 
in Section II-3.1.1.2, the ephemeris of each satellite can' be 
approximated by two simple sinusoids with eight adjustable con- 
stants to describe the phase relationships, amplitudes, bias, 
and frequencies of the latitude and longitude components of 
satellite movement. The lissajou pattern created by the 
equations is sufficiently accurate over a twenty- four hour 
period to eliminate any need for in-flight updating of the 
ephemeris. It will be shown later in this report that satellite 
altitude movement need not be considered. For all these reasons, 
pre-programmed ephemeris data is the recommended technique. 

1.1.4 Compensated vs. Uncompensated Ranging Signal 

The last alternative in the satellite navigation imple- 
mentation tree concerns the choice between ranging tone 
compensation with its associated equipment complexity and ex- 
pense and the uncompensated ranging tone technique which is less 
complicated to implement, but is also less accurate. Ranging 
tone compensation means that the phases of the ranging side f 

tones are adjusted according to some criteria at the source so 
the user may obtain a more accurate position location. It is 
effected by adding one or more ground-based reference terminals 
which tran'spond ranging signals back to the GCC via the satel- 
lites according to a pre-determined sequence. The GCC then * 

advances or retards the phases of the transmitted ranging tones 
according to how the measured phases of the returned signals 
compare with the expected phases as calculated for the particular 
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sector of the transponding reference, and knowledge of the 
approximate ephemeris information used at the aircraft in 
performing the position calculations. 

, The compensated ranging signal, or constant phase plane 
tecrinique, enables the ground control center to partially cancel 
the effects of ionospheric bias and undetected variations in 
satellite-GCC-aircraf t geometry. Changes in the electron 
density distribution of the ionosphere cause a variation of the 
index of refraction of this layer, which in turn chanqes the 
effective path length of VHF radio signals traveling through it. 
Without a constant phase plane technique to pre-distort the 
phase of the ranging signals to accommodate the variable delay 
introduced by propagation through the ionosphere, user aircraft 
will derive inaccurate position locations. Although the rate of 
change of the ionospheric index of refraction is relatively slow, 
the path length of a VHF signal traveling through it may vary by 
as much as a kilometer in a six-hour period. Since a one-km 
error in range translates to approximately 1.3 kilometers of 
aircraft position error, as shown in Section II- 3, the effect 
cannot be neglected. 

In similar manner, variations in the system geometry 
impact the accuracy of the aircraft position location. While 
the major diurnal satellite excursions will be largely accounted 
for by the use of ephemeris data (both on the ground to pre- 
distort the phase of the ranging tones , and in the air to cal- 
culate position) r the approximate form of the data to be used in 
this system makes it impossible to locate the satellite closer 
than within 35 km of its true position. Using the same ratio 
of range error to aircraft position error as calculated in 
Section II-3 (and used above in assessing the impact of 
ionospheric bias) this implies a worst-case aircraft position 
error of 47 kilometers. The results of computer analysis of a 
ranging system with a constant phase plane indicate that this 
worst-case error can be reduced to less than 1300 meters when 
located as far as 1000 kilometers from the nearest reference 
station (see Figure 2-11). Moreover, the typical error at this 
distance for a 35-kilometer satellite error is only 900 meters. 

As the aircraft moves closer to the reference station, this 
error decreases and becomes zero when the aircraft and the 
reference terminal are collocated. The need for the constant 
phase plane technique is clear. Although this feature entails 
more equipment complexity, it is a feasible method of obtaining 
the necessary position location accuracy over a large area. 

1.2 Selected Approach 

1.2.1 Summary of System Characteristics 

The recommended navigation system configuration which 



the Dual Satellite Navigation Experiment will approximate will 
enable any aircraft with the proper on-board equipment to deter- 
mine its own latitude and longitude to within 1500 meters of the 
true coordinates independent of external ground control equip- 
ment. The position computations will be handled by an on-board 
data processor which may be a separate unit or a small portion 
of a computer that is already on-board (as in the case of the 
larger commercial aircraft) . The problem of keeping the air- 
craft synchronized with the ground control center which 
generates the ranging tones has been resolved by the choice of 
an inexpensive, but highly stable, on-board oscillator. This 
oscillator may be set before leaving the ground and is 
sufficiently stable so that no in-flight clock update is 
required for the duration of the flight. 

The ground-based grid which enables the GCC to transmit 
ranging tones of the proper phase to all sections of the 
continental United States shall consist of one ground control 
center located in the South Central U.S. and a network of 10 
to 20 inexpensive reference terminals. The exact- number of 
reference terminals will depend on satellite elevation angles, 
and the desired accuracy of the position location at all points 
on the grid when the proposed system is implemented. Since the 
expense of the reference terminal is minimal compared to the 
cost of the ground control center, severe restrictions need not 
be placed on their total number. 

1.2.2 Surveillance Compatibility 

Aircraft position determination can be effected using 
the range data to a pair of satellites in conjunction with the 
altitude of the aircraft. Position determination on-board the 
aircraft is considered part of the navigation function while 
position determination at a ground station is commonly classi- 
fied as surveillance. In a broader sense, the surveillance 
function also includes provisions for a data link between the 
aircraft and the ground station. 

In both surveillance and navigation systems utilizing 
satellites, the aircraft receives the ranging tones transmitted 
from the ground station and relayed by the satellites. Because 
the transmitted tone package is common to both systems the two 
are compatible. For the recommended system, compatibility is 
further enhanced by the use of a 600-Hz ambiguity resolving tone 
which is ideally suited to the transmission of digital data 
(see Section II-2.1). 

Addition of the surveillance function to the navigation 
system would require that each cooperating aircraft be equipped 
with a transponder to relay the received ranging tones back to 
the satellite. The ground station would require the down-link 
receivers, and appropriate hardware and software for position 



^ determination and digital data processing- While the GCC 

function would be slightly more complex the additional aircraft 
equipment to provide for surveillance would be minimal . 

1.2.3 Extension to Three Satellites 

*. 

Two satellites in synchronous equatorial orbits have 
several limitations. They cannot be used to derive altitude 
information; they do not provide coverage in the polar regions; 
and while they derive excellent longitude information, latitude 
information is very poor in equatorial regions. In the North 
Atlantic region two satellites provide good coverage and 
reasonable accuracies (on the order of one nautical mile) . 

The addition of a third satellite can be used to improve 
coverage at higher latitudes while simultaneously improving 
accuracies at the low latitudes, or it can be used to extend 
coverage in a longitudinal direction. 

If a third satellite is placed in a synchronous in- 
clined elliptical orbit to improve coverage at higher latitudes 
and to enhance position location accuracies at the lower 
latitudes it must be remembered that the satellite will only be 
visible in the northern latitudes for twelve to eighteen hours 
a day depending on the orbit eccentricity. For example, with 
an eccentricity of 0.35, the satellite will be in the northern 
hemisphere approximately 18 hours while for an eccentricity of 
zero it will only be in the northern hemisphere for 12 hours. 

The addition of a fourth synchronous inclined elliptical orbit 
satellite gives twenty-four hour coverage for the same area of 
the northern hemisphere if properly spaced. The third and 
fourth satellites will not be available during the time span of 
the Dual Satellite Navigation Experiment, therefore their 
usefulness cannot be demonstrated by this experiment. For an 
operational system, however, added satellites are recommended 
because of the increased operational reliability and redundancy 
they provide in addition to the accuracy improvements achievable 
with inclined orbits. 

If the third satellite is placed in a synchronous 
equatorial orbit, redundancy is obtained for the same longi- 
tudinal coverage by placing the satellite at a longitude between 
* that of the other two satellites. The third satellite could 

otherwise be placed at a longitude east of the eastern most 
satellite (or west of the western most satellite) to increase 
the area of coverage that would be realized by using any two 
satellites. Position accuracies in the case of using three 
satellites are not significantly improved over that obtained 
with the original two satellites. 
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The system consisting of two equatorial synchronous 
satellites and one inclined synchronous satellite can also be 
used to obtain altitude information during the time the third 
satellite is in the higher latitudes. If the satellite is at 
25°N or higher, altitude accuracies will be reasonable good. 
Altitude measurement error is a function of range error multi- 
plied by a geometry factor. The best altitude accuracies are 
obtained when the aircraft is at a subsatellite point. It is 
necessary for this system to use satellite orbit and ephemeris 
in calculation of position and altitude; for the Dual Satellite 
Navigation system this correction is simply taken care of in the 
ground station data processor. 

It is important to recognize that extending the proposed 
two-satellite configuration to three or more satellites in order 
to extend the coverage area or improve accuracy and increase 
capabilities does not require significant changes to the Dual 
Satellite Navigation System concepts. 

2. POSITION LOCATION STUDY 

2 . 1 Range Measurement 

The time taken by a signal to propagate from a satellite 
to an aircraft is directly proportional to the range between 
them. Measurement of one way propagation delay, however, is 
not possible unless the time of transmission is known. Typical 
techniques that avoid this problem either measure the round 
trip signal delay that occurs for transponded signals or 
measure the time-of-arrival difference between the signals from 
two reference terminals at known locations. 

The measurement of time delay of ranging signals is 
common for essentially all range measurement techniques. There- 
fore, the following paragraph discusses alternate approaches 
used to measure time delay. 

2.1.1 Alternate Measurement Approaches 

Three basic approaches have been used to measure time 
delay by appropriate measurement of the parameters of a general 
ranging signal. The approaches, discussed below, are concerned 
with measurement of time, phase and frequency. 

2. 1.1.1 Time Measurement 

The measurement of time (as time-of-arrival) of a 
signal implies that the corresponding range is nonambiguous . 

This is true for a digital coded signal of appropriate code 
length, as for a pulse system when the period between pulses is 
long enough to be essentially nonambiguous. 
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2. 1.1. 2 Phase Measurement 


Measurement of phase implies that the measurement is 
made on the phase of a sinusoidal or tone signal. Phase 
measurements are inherently ambiguous as the identical answer 
is obtained for every cycle of the tone. Ambiguities can be 
resolved by using other tone phase measurements to obtain range 
information. Time delay is determined from the phase measure- 
ment by dividing by the tone frequency. 

Phase will be computed from quadrature samples of 
each ranging tone. Two quadrature samples of a sinusoidal 
signal of known frequency may be used to measure the phase of 
the location of the samples with respect to the zero crossing 
of the tone. For convenience, the in-phase sample is usually 
generated by the positive going zero crossing of the reference. 
To illustrate the principle, consider the sampling of the tone 
as shown in Figure 2-2. The samples shown are spaced 90 degrees 
apart. The location of the positive going zero crossing of the 
tone with respect to the sample (taken at the positive going 

zero crossing of the reference) is the angle 0 which is arbi- 
trary and unknown. For the signal shown, with arbitrary 
amplitude A, the sample values are: 


S I = A sin 9 ■ (2-1) 

S Q = A sin (9 + 90°) = A cos 0 (2-2) 

The angle 0 can then be calculated as follows : 

-1 

0 = tan x (2-3) 

S Q 

This technique therefore provides a measure of the phase of the 
sinusoid at the instant of Sj with respect to the reference. 

2. 1.1. 3 Frequency Measurement 

Measurement of the frequency of a received signal can 

be used to calculate the desired range. Since the desired 

range (R) can be represented as the integral of the range rate 

(R) plus an initial value of the range (R ) it can be calculated 

o 

using the following equation: 






I 

I 

I 



NOTE: Sj and Sq are quadrature 

samples 

6 is an arbitrary phase 
angle 




Figure -2-2. Quadrature Sampling of a" Sinusoid 
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where 


R = 
R = 


/ R dt + R 

o 



(2-4) 

(2-5) 


ancl c is the speed of light 

a) is the angular frequency of the received signal 
w t is the angular frequency of the transmitted signal. 

The difference in the frequency transmitted by the spacecraft 
and that received by the aircraft is the doppler frequency. 
Measurement of this parameter is discussed in Section II-2.4. 

2.1.2 Range Measurement Accuracy 

The theory of ranging signal measurement can be used to 
directly compare alternate ranging techniques if the measurement 
accuracy is expressed in terms of time delay. The basic theory 
of range measurement considers the concept of signal time 
measurement from the view of a maximum likelihood process which 
expresses ultimate measurement limitations in terms of signal 
and noise characteristics alone. Generally the value of noise 
considered in the analysis is a composite of many effects of 
which thermal noise is but one. Other "noise-like" effects are 
due to propagation disturbances and equipment implementations. 

The maximum likelihood approach was chosen because it is optimum 
in the sense that it minimizes the mean square measurement error. 
The theory of ranging signal measurement has been well covered 
in the literature. Woodard [3] did much of the original work, 
and Skolnik [4] , Wains tein and Zubakov [5] develop the theory 
further and present several interesting applications; the results 
presented here are based primarily on these developments. 

Given a demodulated received signal f(t) contaminated by 
additive white Gaussian noise. 


f(t) = S(t - t q ) + n (t) 


( 2 - 6 ) 
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the problem is to determine the true time of arrival (x Q ) which 

represents the propagation delay from the transmitter to the 

receiver. The optimum (minimum mean square error) estimate 

x of (x ) is found from the matched filter or cross-correlation 
o 

receiver that maximizes the message-signal cross-correlation 
function 


C (x) = /f(t) S(t - x) dt . (2-7) 

2 

From the cross-correlation function the variance (a ) in the 

x 

estimation of x of (x ) can be shown to be 

o 

a 2 = (seconds) 2 (2-8) 

T 0 2 (2E/N q ) 

where N q is the one-sided noise power spectral density and E 
is the signal energy given by 


and 



(2-9) 


( 2 - 10 ) 


where S(u>) is the Fourier Transform of S (t) . 

The parameter 8 is the second moment of the normalized 
signal energy and is thus a measure of bandwidth of the signal. 
For a general bandpass signal of the form 


S(t) = g(t) cos oj Q t 


( 2 - 11 ) 
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where g(t) is slowly varying with respect to ui Q , the Fourier 
transform pair is 


G(u - to ) + G ( u) + u ) 

g(t) cos u t > = s (to) (2-12) 


Then , 


where 


t °° 

to 2 |g(u) | 2 du 


S(u)* du 


n 2 2 
$ +04 

V o 


(2-13) 


u |G(u) I 2 du 


r °° 

j |G(u) | 2 du 

J _oo 


«u; 


(2-14) 


Then, for bandpass and tone signals of the form of 

2 

equation (2-11) the value of 8 , to a very close approximation 
2 

is w . 
o 

2.1.3 Accuracy of Alternate Techniques 

It is instructive to compare the fine range accuracies 

attainable with digital code ranging and with a fixed tone. The 

fine range measurement of a digital code is obtained from the 

measurement of phase of the clock frequency. The accuracy of 

the fine range measurement is given by equation (2-8) using the 

2 

appropriate value of 8 . To make a meaningful comparison, the 
signal energy and the baseband bandwidth must be the same for 
both cases. 

* The digital code ranging signal will have a minimum 

baseband bandwidth (B) requirement of twice the clock frequency. 
2 

The value of 8 will therefore equal the squared value of clock 

2 • 2 
frequency (u ) , where u = B/2. The minimum variance (o^) 
c c c 
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attainable for a digital code ranging system constrained- to a 
signal energy (E) and a bandwidth (B) can then be found from 
equation (2-8) , and the above value for the clock frequency to 


0 2 > __ — 4 (2 _ 15) 

3 2 (2E/N q ) 

The tone approach used for fine range measurement re- 
quires that the baseband bandwidth (B) , be approximately equal 
to the radian frequency of the ranging tone (w o ) . The minimum 
variance (a 2 ) attainable for a tone ranging system constrained 
to a signal energy (E) and a bandwidth (B) can then be found 
from equation (2-8) to be 

o 2 > —± (2-16) 

8 2 (2E/N q ) 

The measurement of ranging tones in a given bandwidth yields 
four times better fine range accuracy than is obtainable with a 
digital code of the same energy in the same bandwidth. The 
ranging tone concentrates all its available signal energy at the 
edge of the frequency band which gives a higher value of 6 than 
can be obtained for a signal with a distributed power spectrum. 

This suggests that a ranging tone be used for best 
resolution at a given signal-to-noise ratio. Choosing an 
appropriate ranging tone is the next step. The choice involves 
severed trade-off factors. The energy-to-noise power density 
ratio of the link between the satellite and the aircraft is the 
first constraint on the tone frequency. From Appendix A the 
worst cast effective carrier-to-noise power density of this link 
(ATS -I to A/C) is calculated to be 30.5 dB-Hz as shown in 
Table A-4. This includes the effect of all error sources. As 
the tone package will have three ranging tones (this choice is 
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justified by probability of ambiguity resolution as discussed 
below) and a carrier for a total of four tones among which the 
power must be shared, the amount of signal energy available to 
one tone will be approximately one-fourth of the total or 
24.5 dB-Hz. The next step is to determine an acceptable value 
for rms range error (a D ) . Figure 2-3 shows the relationship 

■K 

between rms phase deviation (a^) and energy-to-noise ratio for 
a one-second integration time. 

A phase detector integration time of one second will be 
used since this interval is a good compromise between obtaining 
a good phase value average and minimizing the effects of 
satellite movement on phase measurements. From the approxima- 
tion for in the linear region of the curve, the equation 
relating a 0 to a R is 


% 




a 


R 


£ 

a) <J> 


(2-17) 


From this expression, the equation relating the ranging tone 
frequency to E/N q and a R may be derived as 

f = T77 (2 — 18) 

2 it a R (2E/N Q ) 1/z 

where c = velocity of light 

Substitution of 24.5 dB-Hz for the worstrcase E/N and 200 

o 

meters as the worst case rms range deviation results in the 
selection of 10 kHz as a desirable frequency for the fine 
ranging tone. 

2.1.4 Ambiguity Resolution 

Ambiguity resolution is a lower accuracy range measure- 
ment that is necessary to obtain a measure of the total delay 
experienced by the ranging signal. It was shown in Section 
II-2.1 that the optimum ranging signal is a tone. Consequently, 
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STANDARD DEVIATION OF PHASE (a ) (RADIANS) 



Figure 2-3. Standard Deviation of Phase Measurement vs. 
Energy-To-Noise Power Density 
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the optimum ranging signal for ambiguity resolution of the fine 
ranging tone is a tone of lower frequency. The accuracy re- 
quired of the ambiguity resolving ranging system must be 
adequate to resolve, with high probability, the ambiguity 
inherent in the fine range measurement. 

If the highest tone frequency required for accurate 
range measurement is f^, then the next lowest tone frequency 
f^_l is determined by a number of factors. For example, if the 
minimum total transmitted power is the primary criterion, it 
can be shown that a tone frequency ratio R, defined by 

R = Vf^ (2-19) 

is approximately two for maximum probability of correctly 
resolving range ambiguity. If the lowest frequency tone re- 
quired for the coarsest ambiguity resolution is f ^ , then the 
total number of tones required (N) to meet this criterion is 
given by 

i°g[f h /f J 

N = 1 + — — ' (2-20) 

log R 

where R is 2 . The difficulty with this method is that a large 
number of tones are required which complicates the equipment 
design. 

To minimize receiver complexity it is desirable to 
reduce the number of tones required. This can usually be done 
while incurring only a very small penalty in the total power 
required. In the non-optimum case in which N tones are used 
from the lowest tone (f^) to the highest tone (f^) , the phase 
measurement variance allowed per tone can be determined once the 
probability (P) of the error in resolving all of the ambiguities 
is specified. In this case, the probability of error in the 
resolution (P ) of each of the (N-l) ambiguity steps is 
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s 


p 


N-l 


( 2 - 21 ) 


For each step, the probability of error in resolution is 


P 

s 


Probability 





( 2 - 22 ) 


where 



4- h 

variance of the phase measurement of the i 
frequency tone 


o 



variance of the phase measurement of the 
+■ 

(i+1; frequency tone 


and the probability distribution is approximately gauss ian. 

The variance of the phase measurements for each ambiguity tone 
are essentially equal and hence equation (2-22) can be written 
as 


V 


# 


P c = Probability [o.(R 2 + 1) 1/2 <ir] (2-23) 

S ]I 

or 

= Probability [o f (R 2 + 1) ^ //2 <tt ] (2-24) 

N-l r 

2 

Once the phase measurement variance due to all effects is 
known, and the error probability P and number of tones have 
been selected, then the tone ratio can be determined using 
tables of the normal probability function. 

A final factor that impacts the choice of the ambiguity 
resolving signal for the coarsest range measurements is that 
processing of low frequency tones introduces practical dis- 
advantages. The most significant of these is the narrow band- 
width filter required to demultiplex the low frequency tones in 
the frequency domain. These filters are limited in response 
and have considerable time delay. As a result, it is desirable 
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to utilize an alternative technique other than tones for the 
resolution of the coarsest range measurements. 

In a satellite navigation system, the choice of 
ambiguity resolving tone frequencies is dictated by the phase 
measuring technique, and the desire to minimize equipment com- 
plexity by maximizing equipment commonality. Because such a 
navigation system would also be compatible with a satellite 
surveillance system, as discussed in Section II-1.2, one am- 
biguity resolving tone will be recovered from a 600-Hz digital 
data channel which is incorporated in the system to satisfy this 
constraint. Ambiguity resolution of the 600-Hz tone would be 
provided by measuring the time of arrival of the data frame 
synchronization pulse with respect to the on-board frequency 
standard. This technique makes maximum use of the accurate on- 
board oscillator, yet eliminates the need to process low 
frequency tones . 

The need for additional ambiguity resolving tones be- 
tween 600 Hz and the fine ranging tone near 10 kHz can be 
determined from the results of equation (2-24) . This equation 
can be expressed by requiring that 

ff f [R 2 + 1] 1/2 1 f (2-25) 

where F is a confidence factor used to meet the allowable 

probability error for each ambiguity resolution. If a total 

-4 

probability of error (P) of 10 is chosen, then in a three-tone 

system, the probability of error at each ambiguity resolving 
-5 

step is 5 x 10 . The F factor in this case can be obtained 

from tables of the normal probability function [6] to be 4.06. 

The tone ratio R in this case is approximately 4. Hence the 
allowable phase measurement variance for ambiguity resolution is 
approximately 10 degrees on the 600-Hz tone, the fine ranging tone 
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near 10 kHz, and the intermediate tone near 2.5 kHz. In fact, 
the phase measurement variance on the fine ranging tone and the 
600-Hz data channel will be appreciably better than this. This, 
therefore, represents the requirement for the intermediate tone. 

The exact frequency of the fine ranging tone near 10 kHz 
and the intermediate tone needed to insure that the fine ranging 
tone ambiguity can be resolved will be determined from the con- 
straints of the sampling system used in the phase measurement. 

A number of constraints apply to the choice of the sampling 

frequency. First, the sampling frequency (f ) must be of a high 

s 

enough rate to provide two quadrature samples of each ranging 
tone. Second, f must be greater than twice the bandwidth of 
the filters used to present a band limited signal to the sampler 
(in order to prevent aliasing as described in Section II-2.1.5). 

These requirements may be met by choosing an inter- 
mediate tone frequency of 2.4 kHz and a sampling frequency of 
9.6 kHz. The sample rate of 9.6 kHz has the added advantage 
that it is possible to derive sample pulses from this frequency 
which will provide two quadrature samples of a frequency which 
is nine times the fundamental frequency of 600 Hz or 10.2 kHz. 
Hence, the three ranging tone frequencies are 600 Hz, 2.4 kHz, 
and 10.2 kHz. 

The fine grain tone at 10.2 kHz generates ambiguities 
in phase measurement every 98 ysec or every 29.9 km for a one- 
way range measurement. The 2.4 kHz tone increases the distance 
between one-way range ambiguities to 416 ysec or 126.8 km. The 
data subcarrier at 600 Hz further increases the distance between 
ambiguities to 1.66 ms or 507.0 km. The resolution of range 
provided by the period between successive frame sync pulses on 
the digital data channel will be on the order of 300,000 km 
between ambiguities as the message period is designed to be 1 
second. Together, the four time intervals provided by three 
ranging tone frequencies and the repetition rate of the digital 



code will be more than adequate to resolve any range ambiguities 
that the navigation experiment will encounter. 

2.1.5 Implementation 

The demodulated orderwire signal spectrum containing 
ranging tones and digital data is reproduced in Figure 2-4. 



FREQUENCY (kHz) 


Figure 2-4. Demodulator Orderwire Signal Spectrum 


The dotted line at 600 Hz represents the phantom data 
subcarrier that is recovered by processing of the digital data 
, signal. The 600-Hz signal is used both as a range ambiguity 

resolution tone and also for bit synchronization for the digital 
data. Prime emphasis, however, will be on the use of the 600-Hz 
tone for ambiguity resolution. Measurement of the phase of the 
10.2-kHz, 2.4-kHz and the 600-Hz tones will be done by accumu- 
lating binary samples of individually filtered signals, for each 
signal burst, as shown by the block diagram of Figure 2-5. The 
values of the accumulated quadrature samples will serve as in- 
puts to the computer to determine the phase of each tone. The 
range will then be calculated by resolution of ambiguitites . 

The ambiguities in the 600-Hz signal will be resolved by the 
computer by the successive knowledge of bit synchronization, 
message framing and system timing. 

lb The filters shown in Figure 2-5 are necessary to present 

a band-limited signal to the sampler and analog-to-digital con- 
verter. Band limiting is required since sampled data systems 
operating at a sampling rate less than twice the band-limited 
. bandwidth can be contaminated by a phenomenon termed aliasing. 

Signals and noise outside a frequency band that is half the 
actual sampling rate (termed the Nyquist bandwidth) will be 
folded (aliased) into the Nyquist band and cause distortion. As 
filter bandwidths are reduced, however, the filter time delay 
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Figure 2-5. Block Diagram of Phase Measuring Equipment 




increases. The filter time delay causes a delay in acquisition 
of the signals because of filter transient effects. 

A good compromise between sampling rate requirements 
and acquisition delay is made by selecting a 500-Hz filter 
bandwidth. The filter time constant, to a first approximation, 
is 2 milliseconds which permits the attenuation of essentially 
all significant transient effects within six to ten milliseconds 
The sampling rate required for a 500-Hz noise bandwidth filter 
is 1000 samples per second. Since the 500-Hz filters will not 
totally reject all unwanted signals, the sampling technique was 
configured to assure that the accumulated samples taken for 
measurement of each tone did not contain other undesired signals 
from the orderwire spectrum. This can be assured if the tone 
frequencies and integration periods are selected such that: 

1. Harmonics of the sampling frequency do not beat 
(hetrodyne) other tones or data to dc (frequency 
equals zero) 

2. Beat frequencies which do occur fall at zeros 
of the integrator sync function. 

The measurement technique chosen for implementation generates 
periodic samples at a 9.6-kHz rate. The following paragraphs 
describe the process in which the signals are sampled, converted 
and accumulated. Ranging tone sample timing and sample 
identification is shown in Figure 2-6 to aid the following 
discussion. 

10. 2 -kHz Tone 


* 


f 


Every fourth sample taken at the 9.6-kHz rate is used to 
integrate the 10.2-kHz signal. If the first sample chosen is 
arbitrarily designated sample 3^, the sample sequence to be used 
for the 10.2-kHz tone phase measurement is (S^, S,., Sg , • • 

S 357 ) . Samples (S^, Sg, , . . .) constitute the set of in- 

phase samples, while samples (S^, S 21 7 * • • ) are the set 

of quadrature samples. 


The sense (sign) of adjacent samples is negative. 
Accumulation therefore requires that the magnitude of in-phase 
samples be added in the 10.2-kHz in-phase accumulator. Quadra- 
ture sample magnitudes will likewise be added in the 10.2-kHz 
quadrature accumulator. The sum of the samples in the in-phase 
accumulator for 10.2-kHz signals is A^ and the sum of the 
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gr 


samples in the quadrature accumulator is A^q, 


where 


V 


and 



\ I S [1+8 (n-1)] I 
n-1 



45 

n l x ' S [5+8 (n-1)] ' 


9 


2.4-kHz Tone 

The sample sequence to be used for the 2.4-kHz tone 
phase measurement is (S 2 , S^, S 1Q , S^, s ig/ s 19 * • • • S 354 

S 355 ) * 

Samples with even numbered subscripts constitute the 
set of in-phase samples while samples with odd numbered sub- 
scripts are the set of quadrature samples. The 2.4-kHz in- 
phase accumulator will contain the sum A 2I and the quadrature 
accumulator will contain the sum A 2 q, where 


45 

A 2I = I S [2+8 (n-1) ] I 

n=i 

and 

45 

A 2Q = E . I S 13+8 (n-1) J I 
n— l 


600-Hz Tone 

The sample sequence used to measure the phase of the 
600-Hz tone is (S^, Sg , S^ 2 , * ' * S 360^ ' The i n ~Ph ase sample 
set is (S g , S 24 , . . . Sg 5 g) . Alternate sample values are 

negatives of each other which then requires that the accumulator 
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sum magnitudes of samples. Alternatively the first sample can 
be added, the second subtracted, the third added, and so on. 
The 600-Hz in-phase accumulator will contain the sum Agj, and 
the quadrature accumulator will contain the sum Ag^, where 


and 



l l S [8+8(n-l) ] I 

n=i 



n f x f S [4+8 (n-1) ] I 


% 


The values of the approximate sums contained in each 
accumulator determine the phase of the three tones. Once the 
phase of the 600-Hz tone is measured the phase of the 2. 4 -kHz 
tone will be measured such that ambiguities involving whole 
cycles are removed. That is, the zero crossing of the 2. 4 -kHz 
tone nearest the zero crossing of the 600-Hz tone is selected. 

The same process is effected with the 10.2-kHz signal. The 
correct epoch is determined by the 600-Hz tone zero crossing 
and the accurate measure of that time is given by the measured 
zero crossing of the 10.2-kHz tone. If the zero crossing of the 
600-Hz tone corresponding to the end of the last bit of the last 
character of the digital message is chosen, the time delay of 
the signal with respect to the transmitted signal (and hence 
range) can be calculated. 

2.2 Position Location Calculation 

There exists several coordinate systems, each with its 
own set of equations describing the satellites and aircraft 
geometry, which can be used in computing the position of the 
aircraft. The set of coordinates which is most convenient, 
however, in the sense that the associated equations can be 
solved in a straightforward manner is the orthogonal earth- 
centered set oriented as shown in Figure 2-7. In this 

configuration three parameter measurements are required: the * 

range from each of two satellites to the aircraft and the air- 
craft altitude above the surface of the earth. Additionally, 
the input parameters defining the position of both satellites 
and the distance from the center of the earth to the surface *> 

are required in the solution for aircraft position. Although 
the aircraft position solution is developed using two satellite 
reference locations, the solution technique is applicable in 
determining aircraft position location utilizing any pair of 
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Figure 2-7. Two-Satellite Navigation Configuration 


known reference locations from which straight-line range. . 
measurements to the aircraft have been made.. For example , one 
reference location could be on the ground within line of sight 
of the aircraft while the other location could be a single 
satellite. The problem may be concisely stated as follows: 

GIVEN: R1 = measured range between aircraft and 

satellite 1 

R2 = measured range between aircraft and 
satellite 2 

RSI = distance from earth center to 
satellite 1 

RS2 = distance from earth center to 
satellite 2 

RH = distance from earth center to aircraft 

B1 = direction cosine vector to satellite 1 

= Bll 
B12 
B13 

B2 = direction cosine vector to satellite 2 

= B21 

B22 
B2 3 

FIND: BA = direction cosine vector to the aircraft 

= BA1 
BA2 
BA3 

SOLUTION: The components in Bl (the direction cosine 

column vector) are defined as shown in 
Figure 2-7 to be: 

Bll = cosine 

B12 = cosine cxj 

B13 = cosine 


2-28 





* 


« 


If the orthogonal coordinate system X^, X 3 , X 3 is aligned 
such that the X^ axis corresponds with the North/South axis 
of the earth and the X^ axis passes through the Greenwich 
meridian (0° longitude) , then the vector B1 is expressed in 
terms of satellite latitude and longitude as 

Bll = cosine (Latitude) cosine (Longitude) 

B12 = cosine (Latitude) sine (Longitude) 

B13 = sine (Latitude) . 

Vector B2 is defined in a similar manner with respect to 
satellite 2, and vector BA is similarly defined with respect 
to the aircraft. Once the components of the vector BA are 
found, the latitude and longitude are calculated by substitution 
in the equations 


LAT 

= arc sin (BA3) 

(2-26) 

L0N 

= arc tan (BA2/BA1) 

(2-27) 


The derivation of the equations necessary to find the components 
of the direction cosine vector BA begins with the application 
of the law of cosines to the two triangles created by the earth 
center, the aircraft, and the two satellites. For this geometry 

Rl 2 = RH 2 + RSI 2 - 2*RH*RS1* (B1 T *BA) (2-28) 

R2 2 = RH 2 + RS2 2 - 2*RH*RS2* (B2 T *BA) (2-29) 

T 

where Bl indicates the transposed column matrix Bl. This 
matrix is defined by 


B1 T = |B11 B12 B13 I 


(2-30 


The asterisk is used in this case for conventional multiplica- 
tion. Solving equations 2-28 and 2-29 for B1 T *BA and B2 T *BA 


gives 


B1 T *BA 

B2 T *BA 


2 2 2 

RH + RSI - Rl * 

2*RS 1 *RH 

2 2 2 

RH -I- RS2 - R2 

2*RS2*RH 


(2-31) 


(2-32) 


Every quantity on the right hand side of equations 2-31 and 
2-32 is known, so two of the three necessary equations can be 
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derived from the expansion of the dot products on the left hand 
side of equations 2-31 and 2-32. Expansion of these products 
yields : 


B11*BA1 + 

B12*BA2 

+ B13 *BA3 = 

Dl 

(2-33) 

B21*BA1 + 

B22*BA2 

+ B23*BA3 = 

D2 

(2-34) 


The third equation necessary to solve for the aircraft 
direction cosine is derived from the fact that the direction 
cosine vector is a unit vector. This gives 

BA1 2 + BA2 2 + BA3 2 = 1 (2-35) 


The solution to equations 2-33, 2-34, and 2-35 is complicated 
by the fact that equation 2-35 is nonlinear. The equations. 


however , 

can be solved 

l in the following manner 



Step 1. 

Eliminate 

BA2 

: from equations 2-33 and 

2-34. 

Then , 


DEL4 

= 

BA1*DEL1 + BA3*DEL2 



or 

BA1 

= 

(DEL4 - BA3*DEL2) /DELI 


(2-36) 

where 

DEL4 

= 

D1*B22 - D2*B12 


(2-37) 


DEL 2 

= 

B22*B13 - B12*B23 


(2-38) 


DELI 

= 

B22*B11 - B12*B21 


(2-39) 

Step 2 . 

Eliminate 

BA1 

. from equations 2-33 and 

2-34. 

Then, 


DEL5 

= 

-BA2*DEL1 + BA3*DEL3 



or 

BA2 

= 

(-DEL5 + BA3*DEL3) /DELI 


(2-40) 

where 

DEL3 

= 

B13*B21 - B23*B11 


(2-41) 


Step 3. Substitute equations 2-36 and 2-40 into equation 2-35. 


BA3 2 


(-DEL5 + BA3 *DEL3) 

2 

(DEL4 - BA3*DEL2) 

DEL! 


DELI 


(2-42) 


which simplifies to: 

BA3 2 *SQD - 2*BA3*SUM1 + SUM2 = 0 (2-43) 

where SQD = DELI 2 + DEL 2^ + DEL3 2 (2-44) 



# 


SUM1 


= DEL2*DEL4 + DEL3*DEL5 
SUM2 = DELI 2 + DEL4 2 + DEL5 2 


(2-45) 

(2-46) 


Step 4. Solve for BA3 by using the quadratic equation 


BA3 


SUM1 




SUM1 - SQD*SUM2 
SQD 


(2-47) 


BA1 and BA2 may be determined by substitution into 
equations 2-36 and 2-40. Finally, aircraft latitude and 
longitude are determined by application of equation 2-26 and 
2-27. It is worthwhile to note that when the satellites are 
positioned directly over the equator so that B13 = B23 = 0, 
the calculation of the direction cosine BA is simplified to 
solving two simultaneous equations in two unknowns. Equations 
2-33 and 2-34 reduce to 


Bll* 

BAl + 

B12*BA2 = Dl 

(2-48) 

B21* 

BAl + 

B22*BA2 = D2 

(2-49) 

the solution : 




BA1 = 

Dl 

D2 

B12 

B22 

Dl*B22 - D2*B12 

(2-50) 

Bll 

B21 

B12 

B22 

Bll*B22 - B21*B12 


and 


BA2 


Bll Dl 
B21 D2 


Bll B12 
B21 B22 


D2*B11 - D1*B2 1 
B11*B22 - B21*B12 


(2-51) 


Equation 2-35 is then solved for BA3 , and the aircraft latitude 
and longitude are found as stated above. 

2.3 Sector Organization for a Navigation System 

Navigation using range measurements to satellite trans- 
ponders in which the constant phase plane technique is used 
requires the use of a reference terminal at a known location. 

As was discussed in Section II-l.l, the constant phase plane 
technique relieves the navigation system of being required to 
continually supply updated ephemeris data to the aircraft. The 
function of the reference terminal is to provide the ground 
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control center (GCC) with the round trip range information 
needed to properly adjust the phase of the ranging signal to 
implement the constant phase plane concept. 

As the objective of a satellite navigation system is to 
provide coverage over the widest possible area, multiple 
reference terminals will be required. Each reference terminal 
will provide the GCC with ranging phase information that will 
enable the GCC to generate a set of ranging phase planes that 
are appropriate at each reference terminal for the assumed 
satellite ephemeris data. This control of the ranging phase 
planes in relation to the various reference terminals will be 
performed in a time multiplexed manner. Thus, the ranging tone 
transmitted at any given time will only be accurate for that 
sector containing the reference terminal that is controlling 
the GCC tone transmissions. Aircraft operating in other 
sectors will obtain position locations which are degraded 
according to how far removed they are from the phase adjusted 
sector. By time multiplexing, however, aircraft in all sectors 
will have one opportunity during each cycle through all 
reference stations to obtain the optimum position location for 
their area. 

The phase control information required at the GCC will be 
obtained from the various reference terminals by using trans- 
ponding reference terminals that relay the ranging signal back 
to the GCC via each satellite in time division multiplexed 
channels common to all reference terminals. The time sharing of 
the GCC tone phase control among sectors will be organized such 
that the phase correction measured in the first time slot for 
sector k will be used to update the stored estimate for sector 
k derived from earlier measurements. Transmission of the 
ranging phase applicable to sector k will then be made in the 
second time slot. During the second time slot, the phase 
correction that applies to sector (k+1) is being updated. 

The time slot duration allocated to each sector is 
regulated by a number of factors. These include: 

1. The number of reference terminals to be interrogated. 
i 2. Sector transfer time 

3. Phase detector integration time 

4. The acquisition time of the GCC phase-lock-loop 
receiver 

5. The maximum allowable error introduced by configura- 
tion changes between the sector k update calculation 
and the subsequent transmission of that update to 
sector k. 



d 


* 


The number of sectors depends on the area of earth 
coverage adjacent to each reference station that enables a user 
to fix his position with acceptable accuracy. This is discussed 
in detail in Section II-3.1 where it is shown that a reference 
terminal separation of 400 km would provide adequate position 
location accuracy. In the recommended navigation system, 
reference terminals will be interrogated sequentially by the 
ground control center via the digital data link. This technique 
insures that the appropriate terminal transponds in the correct 
time slot as determined by the GCC clock, and avoids the prob- 
lem posed by the alternative sector transfer control method 
which requires updating reference terminal clocks to counter- 
act long term drift effects. Sector control by interrogation 
requires that the GCC wait for a return to confirm that the 
message is properly received. Approximately 0.5 second will be 
needed to complete the interrogation of a reference terminal 
from the ground control center via the satellite. 

Phase detector integration time relates to the available 
C/N q at the GCC and is the primary constraint affecting the 
minimum time between successive time slots. With an available 
C/N q of 24.5-Hz, the recommended integration time is one second 
(as shown in Section II-2.1.3). 


Phase- lock- loop receiver acquisition time has a great 
impact on the width of the guard bands surrounding the phase 
detector integration time band. When one reference terminal 
transponder shuts off and the next one starts, the signal at the 
GCC may be totally absent for as long as 40 ms. Hence, the GCC 
will lose signal lock and will have to reacquire it again. 

Using a fast acquisition phase-lock-loop technique, acquisition 
can be easily acquired in 50 ms. 


Errors introduced by configuration changes between the 
phase correction measurement for a given sector and the sub- 
sequent transmission of the correction during the next time 
slot establish the maximum time slot duration. The rate at 
which the navigation configuration changes can be computed by 
finding the slope of the curve portraying the diurnal variation 
of the range from the satellite to the ground control center. 

A graph of this function for the ATS-I satellite is shown in 
Figure 2-10 in Section II-3.1. Modeling the range variation as 
a sine wave gives 

R (t) = d sin tot (2-52) 

o 

where to = 7.26 x 10 -5 radians/second 

d = 200 km 


= time in seconds . 
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t 



The rate of change of range is found from the time 
derivative of R s (t). This results in 

dR (t) 

= wd cos wt (2-53) 

dt 

The maximum value of the derivative occurs when the cosine is 
unity. Hence the maximum rate of change of range is 

wd = (2 x 10 5 m) (7.26 x 10~ 5 sec" 1 ) 

- 15 m/sec 

Assuming a maximum error due to this factor of 75 meters, 

, no more than five seconds may elapse between the completion of 
the phase correction measurement, and the completion of the 
phase correction transmission. Taken together, these con- 
straints on the slot time duration indicate that the optimum 
slot organization is two one-second guard bands surrounding a 
one-second integration time band for a total of three seconds . 

2.4 Range-Rate Measurement 

2.4.1 Basic Technique 

The measurement of range-rate between the satellite and 
cooperating aircraft provides data that can be used to deter- 
mine aircraft velocity and range. To implement the measurement, 
the received carrier frequency from the satellite would be 
measured with respect to the stable on-board oscillator. The 
difference between this measured frequency and that known to 
exist without doppler is then a measurement of the doppler 
frequency which can be converted to range-rate. This technique 
assumes that the satellite transmitter frequency is known. 

This will be true if a coherent satellite transponder is used, 
or if the frequency is relayed to the aircraft from the, ground 
over the data link. For simplicity in this discussion, it will 
be assumed that a coherent satellite transponder is used. 

The method by which the frequency, as received at the 
aircraft, is measured with respect to the on-board oscillator 
is shown in Figure 2-8. The stable oscillator is used as the 
reference in the counter and also as the source for the multi- 
plier in the receiver down- converter. The multiplier ratio is 
assumed to be M, while the actual oscillator frequency is 

(f + Af) , and is in error by Af from its assumed value of f . 
o o 

The frequency that is presented to the counter (f ) is thus 



SATELLITE 



r 


Demodulator 

Output 


Frequency 

Measurement 



Figure 2-8. Doppler Measurement System 
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given by 


f c - M(f o + “> - f RA (2 - 54> 

where is the frequency received by the aircraft. This 

frequency is given by 


where 


f RA + 

f TS = satellite transmit frequency 
R a = aircraft-to-satellite range-rate 
c = velocity of light. 


(2-55) 


If the frequency measured by the counter is F, then it is re- 
lated to the true frequency f at this point by 

c 

f c = F ( f Q + Af)/f o (2-56) 


The following relationship can then be derived between the 
system parameters and the oscillator error. 

f + Af 

1 + R a /c = [M - F/f Q ] (2-57) 


To simplify this expression it is desirable to define 

F = F + AF 
o 

where F = count expected with no doppler and an exact 

oscillator 

AF = count increment measured due to doppler and 
oscillator error 

Then, 


Mf - 
o 


F 


o 



(2-58) 
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and Equation (2-57) can be written in the form, 

1 + R a /c = [1 + Af/f Q ] [1 - AF/f Tg ] (2-59) 


This can be written approximately as 

R a c = Af/f Q - AF/f TS (2-60) 

The effect of oscillator instability upon the range- 
rate measurement is now apparent. The error in range-rate is 
given by 


Range-Rate Error = Af c/f (2-61) 

g 

An oscillator error of only one part in 10 will result in a 
velocity error of 6.7 mph. The measurement of accurate 
velocity thus requires a very stable oscillator. 

2.4.2 Oscillator Up-Date 

The oscillator stability requirement for accurate range 
rate measurement is sufficiently severe as to require some 
method whereby the oscillator drift can be corrected, and the 
oscillator frequency up-dated. Not only is this technique de- 
sirable for a satellite navigation system, but it is also 
necessary in certain cooperative collision avoidance systems 
(CAS) being proposed at this time. 

A technique is proposed by Texas Instruments to imple- 
ment on-board oscillator frequency up-date. This technique 
assumes that the aircraft satellite navigation equipment is also 
part of an aircraft satellite surveillance system. This being 
the case, it will be shown that the doppler frequency measure- 
ment made on-board the aircraft can be used at a ground station, 
together with a measurement of the received aircraft trans- 
mitter frequency as relayed by the satellite, to compute the 
aircraft oscillator frequency at the ground station. This data 
can then be relayed back to the aircraft and any desired com- 
pensation can be made. 

A simplified block diagram of a single oscillator trans- 
ponder is shown in Figure 2-9. This arrangement also provides 
approximate doppler frequency compensation of the aircraft 
frequency as received at the satellite. The equation relating 
the aircraft transponder transmit frequency (f TA ) and the 
receive frequency (f^) is 

f TA " <" + N > f X - f RA (2 - 62 > 
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where f x is the on-board oscillator frequency, and M and N are 
multiplier factors. The value of f x is not known a priori and 
may be computed using this method. The frequency f x is equal 
to the nominal oscillator frequency f plus the error Af. The 
IF frequency in the aircraft terminal is measured as for the 
range-rate measurement. If the count obtained is F, then the 
following relationship applies, 

Mf X - f RA - F V f o < 2 - 63 > 

The measurement F is relayed to the ground station. 

If it is assumed further that the relay satellite is 
coherent and remodulating, then other frequencies in the system 
can be computed. These are listed below. 


Hi 

£ 

H 

f TS 

[1 

+ R A /cl 

(2-64) 

f = 

TS 

f TG 

S 1 

[1 + R g /c) 

(2-65) 

Hi 

s 

II 

f TA 

[1 

+ V o! 

(2-66) 


where 

f TS = satellite-to-aircraf t transmit frequency 

f T g = ground station-to-satellite transmit 
frequency 

fps = satellite-from-aircraf t receive frequency 
S 1 = satellite translation frequency ratio 
R = range-rate of satellite-to-ground station 

VJ 

9 

and c = velocity of light 

The frequency that is measured at the ground station is 
r the channel frequency that exists after demodulation of the re- 

ceived satellite carrier signal. This channel carrier frequency 
Af_.„ is derived in the coherent satellite by down-converting 

the received aircraft transmission with respect to the satellite 
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synthesizer derived frequency. In a coherent satellite trans- 
ponder, all synthesizer frequencies are derived from a reference 
frequency transmitted by the ground. The channel frequency as 
measured at the ground can thus be written as 

Af RG = f RS - f TG S 2 [l + V Cl (2 - 67) 

The constant S 2 relates the satellite aircraft receiver band 
down-conversion frequency to the synthesizer reference frequency 
received in the satellite. 

In the Equations (2-62) and (2-67) the unknowns are 
R, , R„ and f„. The term R-, will be derived by making doppler 
measurements, as described above, with respect to a reference 
terminal located at the ground station. In this case, R equals 
R,_. When R_ has been computed, the Equations (2-62) to (2-67) 
can be solved for R A and f x and hence the on-board frequency can 
be measured. The results of this measurement would then be 
relayed to the aircraft so that its oscillator frequency can be 
updated. 


If a coherent translating satellite is used, then the 
aircraft oscillator frequency can still be computed. The 
equations involved, however, will be slightly different. If 
a multiple oscillator translating satellite transponder is 
used, then additional measurements would be necessary to elimin- 
ate the frequency uncertainty of these oscillators in the 
range-rate equations. 

3. POSITION LOCATION ERROR ANALYSIS 
3.1 Error Sources 

The recommended ranging system is designed to obtain an 
effective ranging accuracy of 1500 meters. In actuality, the 
total range measurement error is caused by the combined effects 
of several independent error sources. The purpose of this 
section is to describe the effects and derive values of the 
errors caused by each source. 

Error sources have been separated into four categories 
as follows: 

(1) Geometric errors introduced by the diurnal position 
variations of the satellites 
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(2) Propagation errors consisting of ionospheric and 
tropospheric effects, multipath influence, and un- 
certainty in the speed of light 

(3) Measurement errors caused by system noise and 
ground processor inaccuracies 

(4) Delay errors caused by uncompensated delays of the 
aircraft, reference terminal, satellite transponders, 
and the ground processor. 

Ionospheric and tropospheric propagation errors are given for a 
worst-case elevation angle of 10 degrees. Compensation for a 
large portion of the absolute range errors will be accomplished 
by the use of a differential measurement technique which 
utilizes range measurements to a reference terminal. Range 
error values presented for each error source are essentially 
worst-case values and are considered conservative, with the 
values seldom being exceeded. Several error values reflect 
improvements and innovations recommended by Texas Instruments. 
The exact nature and consequence of each error source are dis- 
cussed in the following paragraphs. 

3.1.1 Geometric Position Location Errors 

The position of most synchronous satellites is not 
completely stationary with respect to the earth because of 
satellite orbital inclination of 2 to 3 degrees. During the 
course of one 24-hour day, the range between a point on the 
earth's surface (e.g., the ground station) and the satellite 
can vary considerably. Figure 2-10 shows the variation of the 
range from Dallas to ATS -I during the period between 1500 hours 
on April 2, 1970 and 1500 hours on April 3, 1970. From a peak 
range of approximately 39,480 km at 0800, the range decreases 
300 km 12 hours later when it reaches the minimum value of 
39,180 km. By applying the position location equations of 
Section II-2.2, the result of a range change unsuspected by the 
user can be demonstrated. Because it simplifies the position 
location equations, it is assumed that the satellites are 
directly over the equator at longitudes of 50°W and 150°W. 

Using these satellite positions and Equations (2-26) , (2-27) , 
(2-35) , (2-48) , and (2-49) , as shown in Section II-2.2, the 
maximum range gives an aircraft position of 37.63°N, 99.14°W. 

The minimum range produces the aircraft position of 40.71°N, 
101.27°W. (For these calculations, the second satellite is 
assumed to be stationary, and the value of D2 in Equation (2-49) 
is arbitrarily chosen as 0.50.) The distance between the minimum 
and maximum range locations is approximately 400 km. It is 
apparent that without some method to compensate for this large 
geometrically introduced aircraft position error, the accuracy 
of any satellite navigation system is severely compromised. 
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Correction for these variations can take a number of 
forms , the most promising of which is to control the phase of 
the transmitted ranging tones in such a manner as to hold the 
signal phase constant at a reference receiver. This technique 
can be applied to navigation by producing a constant and 
stationary phase plane in the vicinity of the satellite, thus 
creating an effectively stationary satellite reference station. 
It can be implemented by using a reference terminal located at 
the control center, or any other judicious location on the 
ground, and adjusting the phase of the ranging signals generated 
by the satellites so as to keep a constant phase at the refer- 
ence terminal on the ground. 

While this technique provides increased accuracy when 
compared to the uncompensated system, the improved performance 
is degraded as the distance between the user and the reference 
terminal is increased. This deficiency can be corrected by 
providing multiple reference terminals to extend the area of 
coverage or, in other words, reduce the maximum distance a 
user can be from the reference terminal. By time multiplexing 
the ground control center between reference terminals, the 
equipment needed at each terminal can be reduced to a single 
transponder. The major complication introduced into the system 
by adding simple transponding reference terminals which are 
not collocated with the ground control center is that the con- 
stand phase plane transmission must be a round-trip path from 
the GCC to the reference station through the satellite and back 
so that the signal processing can be done at the GCC. This 
doubles the transport delay time in the phase control loop com- 
pared to the case when the reference terminal and the GCC are 
in the same location, and means that the phase at the reference 
terminal must be inferred from the average phase shift over the 
two-way link. The accuracy of the phase adjustment depends on 
the stability of the system configuration (i.e., the rate of 
change of ionospheric bias and the rate of change of satellite 
position) , the phase measurement accuracy, and the delay before 
the computed phase adjustment is actually transmitted to the 
sector. 


3. 1.1.1 Computer Analysis of Geometric Error 

Given the magnitude of aircraft position location 
error resulting from uncompensated satellite movement, and a 
constant phase plane technique to counteract this affect, the 
next step is to determine how completely the constant phase 
plane technique cancels geometrically introduced aircraft 
position location errors. This has been done by computer 
simulation of the satellite navigation system. 

The computer model first precisely locates the air- 
craft in the zone around the ground control center situated at 
Dallas, Texas. Then the two satellites are perturbed from the 
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positions assumed for each by the aircraft, and the aircraft is 
located by satellite ranging using the equations developed in 
Section II-2.2. Next, the true aircraft position and the cal- 
culated position are compared and the maximum, minimum, and 
root mean square aircraft position location errors are cal- 
culated as a function of distance and direction from the ground 
control center. The satellite perturbations are categorized by 
altitude, longitude, and latitude. Figure 2-11 is a sample of 
the output produced by the program. The first two columns give 
the perturbation of each satellite in kilometers. These 
perturbations are made about the assumed satellite positions 
with the positive sign indicating movement away from the earth 
surface in the case of altitude, or increasing latitude and 
longitude in the conventional sense. The variable RGCC is the 
distance of the aircraft from the ground control center in 
kilometers. The next three columns list the maximum, minimum 
and root mean square aircraft position error for each radius. 
Inside the parentheses following EMAX and EMIN is a number in- 
dicating the direction of each error. This number shows the 
compass heading as referenced to the GCC, as illustrated in 
Figure 2-12. 

The computer output reproduced in Figure 2-11 is of 
particular interest because the altitude and longitude cases 
represent the worst-case diurnal position excursions for the 
ATS -I satellite. The 50 -km satellite latitude perturbation for 
both satellites, while not the maximum diurnal position excur- 
sion for either satellite, shows that the worst-case position 
error for satellite latitude perturbations occurs when both 
satellites are perturbed in the same direction. This contrasts 
with the altitude and longitude cases where the worst case 
occurs when the satellites are perturbed in opposite directions. 

Another point to notice is that the errors produced 
by satellite perturbations are symmetrical. That is, moving 
ATS -I in a positive direction and ATS-III in a negative direc- 
tion results in the same aircraft error as does moving ATS-I 
the same amount negatively and ATS-III the same amount 
positively. Finally, the worst-case altitude error produces a 
maximum aircraft error of only 190 meters at a distance of 
1500 km from the ground station. This is sufficiently accurate 
so that ephemeris data to correct for altitude perturbations 
will not be required. 

The error results for the most disadvantageous satel- 
lite perturbations will now be plotted on log-log axes to aid 
in the determination of a mathematical equation to describe the 
variation of aircraft position location error as a function of 
aircraft distance from the GCC. Choosing the worst-case satel- 
lite error perturbations insures that the aircraft position 
location error equations will not assess the magnitude of the 
geometrically introduced position errors too optimistically. 

Root mean square values will be used to derive the error 
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Figure 2-11 (a). Position Error Simulation Program -Altitude Perturbation (km) 
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Figure 2-11 (b). Position Error Simulation Program - Longitude Perturbation (km) 
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Figure 2-11 (c). Position Error Simulation Program - Latitude Perturbation (km) 
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Figure 2-12. Computed Error Direction Key 


equations. The rationale for choosing ERMS instead of EMAX or 
EMIN lies in the fact that the root mean square value is the 
standard deviation of the error, giving a measure of the 
fluctuation of the error about its mean value. Like the 
analogous electrical quantity, rms volts, the rms aircraft 
position location error is the "ac" component of the error. 
Figure 2-13 confirms that ERMS is a reasonable estimate of the 
typical aircraft position location error as it clearly matches 
the behavior of both EMAX and EMIN in terms of slope and 
linearity, but falls closer to the maximum error reflecting the 
prevalence of the larger error values for a given radius. 

The linearity of the aircraft position error as a 
function of distance from the ground control center when plotted 
on log-log axes is shown in Figure 2-14 where four values of the 
parameter Sa (the magnitude of the worst-case satellite latitude 
error perturbation) are plotted. From the figure, it is clear 
that multiplication of the satellite perturbation by a constant 
results in multiplication of the aircraft position location 
error by the same constant. The linearity of this relationship 
combined with the fact that satellite altitude, longitude, and 
latitude are independent makes it possible to write the follow- 
ing equation for aircraft position error as a function of 
worst-case satellite perturbation, and distance of the aircraft 
from the reference terminal (which is located at the GCC in this 
simulation) . 


E = .001 S 2 R 1,065 + .021 S. R 1 * 023 + .034 SR 1,052 (2 ~ 68) 

A A <p 
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Aircraft Position Error (Meters) 


0 


1000 



Figure 2-13. Aircraft Position Error vs Distance From GCC 

For Worst-Case Satellite Altitude Perturbation 
of 50 km 
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RMS Aircraft Position Error (km) 



Figure 2-14. Aircraft Position Error vs Distance from GCC 
For Satellite Latitude Perturbations 
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* 


where 


E - aircraft rms position error in meters 

S = magnitude of satellite worst-case perturba- 
tion in km with subscripts : (A) indicating 
altitude perturbation, (X) indicating 
longitude perturbation, and (<f>) indicating 
latitude perturbation. 

R = distance of aircraft from ground control 
center in km. 


The exponents of R account for the slope of the 
straight line approximating the distribution of data points 
plotted on the log-log axes . This equation depends on the 
position of the satellites. In this case, both satellites were 
placed over the equator with one at 50°W and the other at 150 °W. 
To modify the equation for a different satellite placement re- 
quires that the new coordinates be inserted in the computer 
simulation and a new set of error curves generated. It is 
important to note that the error predicted by this equation re- 
sults from the most detrimetal combination of satellite 
perturbations and may be as much as a factor of three more 
pessimistic than the typical rms error resulting from a more 
optimum conjunction of satellite position errors. 

The equation for aircraft error magnitude as a 
function of satellite position error and distance from the GCC 
was used to plot error contours showing aircraft error as a 
geographic function of position in order to determine the 
location of secondary reference stations. The data for the 
maps shown in Figures 2-15 and 2-16 was obtained by re-formatting 
and sorting the output of the aircraft navigation simulation 
program. The striking feature about these maps is the way that 
the error falls in bands of latitude. This means that reference 
stations need to be placed only on one line of constant longi- 
tudes, north and south of the ground control center to adequately 
compensate the phase of the ranging tones over a longitudinal 
distance of 1500 to 2000 km. The other conclusion which can be 
drawn from the maps is that the bands of constant aircraft error 
are more densely packed south of the ground control center than 
north of it, particularly in the case of satellite-latitude 
generated errors. This suggests that the ground station which 
generates the ranging tones is best located at the extreme 
southern end of the area of navigation coverage. The maps are 
drawn for the worst-case satellite error perturbations and for 
satellite excursions of only 50 km. In the case of satellite- 
latitude perturbation induced aircraft errors, 50 km is only a 
fraction of the expected satellite excursion. Since ephemeris 
data will be used (as described in the following section) to keep 
each user more up to date on satellite movement, basing the 
reference station placement on error contours generated with 
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Error between 
50 and 100 m 


Error between 
200 and 500 m 



Error between 
500 and 1000 m 


Figure 2-15. Aircraft Position Error Contours For 
Satellite Longitude Errors of 50 km 
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Error less than 100 m 



Error between 
100 and 200 m 



Error between 
200 and 500 m 



Error between 
500 and 1000 m 


Error between 
1000 and 1500 m 


Figure 2-16. Aircraft Position Error Contours for Satellite 
Latitude Errors of 50 km 
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50 km satellite errors is valid. If the ephemeris data pro- 
vides the user with satellite positions with a 50-km accuracy, 
reference stations must be placed no further than 350 km apart, 
south of Dallas, or 400 km apart north of Dallas (if the ground 
control station is located at Dallas) in order to hold the air- 
craft error resulting from satellite geometry to less than 
500 meters. As a result of this error constraint, only 10 
simple transponding reference stations would be required to 
provide coverage over the whole region of the United States 
if both satellites were placed over the equator at 50°W and 
150 °W longitude. 

3. 1.1.2 Ephemeris Data 

Geometrically introduced aircraft position errors can 
also be lessened by providing users with ephemeris data des- 
cribing the position of each satellite. Two alternative 
approaches exist. The ephemeris data can be in the form of a 
table, giving the satellite latitude and longitude for a list 
of times , or it can be a pair of equations approximating the 
latitude and longitude of each satellite.! The table approach 
is potentially more accurate because the increment between each 
time that satellite coordinates are provided may be reduced 
without limit. The disadvantage of use of tables is that, a 
table of ephemeris data accurate enough to pinpoint a satellite 
within 10 km of its true position quickly approaches an im- 
practical size if constrained by small storage capacity of the 
on-board data processor. Hence, two equations for each 
satellite will be derived to provide ephemeris data. 

The constants governing the exact form of the 
equation will vary depending on the particular behavior of the 
satellite involved. The process of setting up the equations is 
quite straight forward, however. The diurnal position vari- 
ation of a satellite is well approximated by a two dimensional 
lissajous figure. The frequencies and phases of the two 
orthogonal sine waves which are components of such a lissajous 
figure describe the diurnal behavior of the latitude and longi- 
tude coordinates. Figure 2-17 is a picture of the path of the 
subsatellite point of ATS-III for one day in July, 1970. 

ATS-III was still drifting eastward then, at the rate of 0.06° 
per day. The drift is easily accounted for in the ephemeris 
equation for longitude, although it is not expected that the 
satellites will be drifting when the experiment is carried out. 
The two sine wave approximations are 


1. Ephemeris data will not be provided for satellite altitude 
because altitude variations have an insignificant affect 
on aircraft position error (see Section II-3.1.1) 
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Latitude (Degrees; 




























<j> = LAT = 0.94 cos 2ir(t/144 + 1/6) 


(2-69) 


A = LON = 0.04 sin 2ir (t/144 + 1/6)- t + 64.24 (2-70) 

where t is marked off in units of ten minutes. For the more 
typical figure-eight satellite path the latitude sine wave 
component would contain a sine wave oscillating at twice the 
rate of the longitude frequency oscillation in order to account 
for two maximums and two minimums that occur every cycle. The 
accuracy of this approach is calculated to be +0.05° of latitude 
or longitude which is approximately +35 km at the satellite. 
Implementation of two sets of these equations on the aircraft 
data processor, requiring only five constants each, is easily 
within the capabilities of the machine memory. As the satel- 
lite ephemeris varies from the behavior predicted by these 
equations over the span of a week, provision will also be made 
to update the equations by inserting new constants. 

3,1.2 Propagation Errors 

The propagation effects fall into two classes: those 
that occur along the satellite-to-aircraft path due to the 
variation of refractive indexes of the troposphere and iono- 
sphere, and those that are due to reflection of the electro- 
magnetic wave by the earth's surface. Since the indexes of 
refraction of both the troposphere and ionosphere vary with 
height above the earth's surface, the direction of the wave- 
front is changed as it passes through the medium, and a range 
error is introduced because the actual transmission path is 
longer than the line-of-sight path. In order to gain adequate 
coverage, a fairly large sector of the earth's surface is 
illuminated by the satellite antenna, and some of the energy 
incident upon the earth will be reflected toward the aircraft 
and interfere with the direct-path signal. These two classes 
of propagation effects are discussed in detail in this section. 

3. 1.2.1 Multipath Effects 

As described in Appendix C, multipath interference 
introduces both amplitude and frequency modulation to the 
direct path signal. The net result is amplitude fading and a 
time-varying phase shift that introduces a ranging error. 

Bergemann and Kucera [7,8] , have obtained amplitude 
fading distributions from 135.6-MHz transmissions between ATS"-! 
and ATS-III and aircraft equipped with a Dome and Margolin 
SATCOM antenna, and Kucera [9] has obtained fading spectra for 
these transmissions. Since their results are directly applic- 
able to this program, they were used in analyzing the multipath 
effects in Appendix C. 
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The aircraft signal fading parameters are summarized 
in Table 2-1. For satellite elevations between 0 and 30 degrees 
the fading is parimarily specular, and the fading frequency is 
generally less than 0.3 Hz. For higher satellite elevations, 
the fading is parimarily diffuse although some specular fading 
may be present. The bandwidth of the diffuse fading spectra is 
only a few Hertz. 



Table 2-1. 

Multipath Fade Parameters 


Satellite 

Elevation 

Type of 
Fading 

RMS* 

Worst-Caset 

Fading 

Bandwidth 

0° - 30° 

Specular 

-2 dB 

-7 dB 

<0.3 Hz 

30° - 50° 

Diffuse 

-1.6 dB 

-3.8 dB 

<5 Hz 

O 

O 

r* 

o 

o 

in 

Diffuse 

-1.1 dB 

-2.3 dB 

<5 Hz 


* The diffuse fade margins are 90 percentile values 
t These are 99 percentile values 


The differential tone demodulation process described 
in Appendix C minimizes the specular multipath phase error be- 
cause the magnitudes of the acquisition/reference (A/R) tone and 
ranging tone phase errors will be approximately equal. The 
diffuse phase error, however, will increase somewhat because 
the A/R tone and ranging tone phase noise will be additive. The 
multipath ranging errors for ranging tone frequencies of 600 Hz, 
2.4 kHz and 10.2 kHz are given in Table 2-2. 

Table 2-2. Multipath Ranging Errors 



Satellite 

Elevation 

Type of 
Fading 

Nominal 

600Hz 

Ranging 

2.4kHz 


ft orst-Case Ranging Erro: 
600Hz 2.4kHz 10.2kHz 

* 

0 - 30° 

Specular 

0 .76km 

0 .62km 

0 .38km 

2 . 25km 

3.4km 1.7km 


30° - 50° 

Diffuse 

9.7 km 

2.4 km 

0 . 56km 

- 

- 

♦ 

50° - 70° 

Diffuse 

6.8 km 

1.7 km 

0 . 39km 

— 

— — 
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3. 1.2. 2 Ionospheric Effects 


As described in Appendix B, the ionosphere reduces the 
group velocity of an electromagnetic wave, and thereby increases 
the path length relative to free space. This path length - or 
ranging - error is a function of the integrated electron density 
and satellite elevation, and occurs on both the GCC-to-satellite 
and satellite-to-aircraft links. Mean ionospheric range errors 
for several conditions of satellite elevation and integrated 
electron density are given in Table 2-3, where I = 2 x 10 16 

2 

electrons/m corresponds to a night-time minimum electron 
• 17 2 

density, I = 3 x 10 electrons/m corresponds to a typical day- 

17 2 

time maximum value, and I = 8.6 x 10 electrons/m is a worst- 


case value. It is evident from the figures in this table that 
ionospheric corrections must be introduced in order to avoid 
substantial ranging errors . 


Table 2-3. Mean Ionospheric Range Error* 


Satellite 
Elevation 
Angle (e) 

I = 2x10 

16 elec 

T~ 

m 

I = 3x10 

17 elec 
—T~ 

m 

r - o c ini’ 7 elec 
I = 8.6x10 y- 

in 

One-way 

Path 

Error 

Two-way 

Path 

Error 

One-way 

Path 

Error 

Two-way 

Path 

Error 

One-way 

Path 

Error 

Two-way 

Path 

Error 

e = 20° 

-0 . 16km 

-0 .32km 

-2 . 6km 

-5 . 1km 

-6 . 5km 

-13km 

e = 45° 

-0 . 12km 

-0.25km 

-2km 

-4km 

-5. 3km 

-10.7km 

e =60° 

-0 . 11km 

-0 .23km 

-1 . 8km 

-3 . 8km 

-4.9km 

-10km 


* Satellite Elevation at GCC = 25° 


Irregularities in the structure of the ionosphere also 
introduce signal amplitude fluctuations or scintillations. The 
scintillations are produced by regions with irregular boundaries, 
and both the magnitude of the scintillation and time of occur- 
rence are quite variable. In general, mid-latitude scintilla- 
tions are neither as large nor occur as frequently as do 
scintillations at equatorial and high latitudes [10] . 

Mid-latitude VHP scintillations were measured by 
Aarons and Whitney [11] . Their data gave a 90 percentile 
scintillation fade margin of 1.3 dB and a 99 percentile 
scintillation fade margin of 9.4 dB . The distribution of the 
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periods of the scintillation fades is given in Table 2-4 where 
it can be seen that the scintillation fading period will 
generally exceed one second. Since the fading period is so 
large, compensation for the accompanying ranging error can be 
readily achieved. 

Table 2-4. Mid-Latitude Scintillation Fading Periods [11] 
Fading Period Percentage of Occurrences 


t-3 

to 

< 5 sec 




43 

6 

sec < T 
— s 

< 

14 

sec 

29 

15 

sec < T 
— s 

< 

59 

sec 

28 


3. 1.2. 3 Tropospheric Effects 

Since the index of refraction of air (n) is nearly 
unity, varying from an average of about 1.00032 at the surface 
to 1.000005 at 100,000 feet, the N unit, where 

N = (n-1) X10 6 

is commonly used for refractivity studies. The variation of N 
with altitude is given by [12] 

N(h) = N + (h-h) AN, h <h<h +lkm 

S S o o 

N (h) = N 1 exp [ln{ (N^^/105) (h-h s -l) }/(h s -8) ] , h s +lkm<h<9km 
N (h) = 105 exp [0.1424 (9-h) ] , h>9km 

AN = -7.32 exp [5 . 577xl0 -3 Nl 

w 

where h = height in kilometers 

h s = height of the surface above MSL in kilometers 

N = surface value of refractivity 
s 

and N = N , + AN 

1 s 
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For elevation angles (i)j) greater than 10 degrees, the 
tropospheric refraction error is given to within three percent 
[12] by 


AR tropo = ( 10_6 /sin^) N(h) dh * 

The range error determined from the above equation is plotted 
in Figure 2-18 as a function of satellite elevation with 
respect to the aircraft altitude for N g = 400. For an aircraft 
flying at 20,000 feet or higher, the maximum tropospheric 
range error will be less than 6.5 meters. Since this error is 
much smaller than the multipath and ionospheric ranging errors, 
it may be neglected. 

3. 1.2.4 Velocity of Liqht Effect 

I 

Since the basic side tone range measurement process of 
this satellite navigation system is dependent upon conversion 
of round-trip time delay in the range signals to an equivalent 
aircraft range, a range error component proportional to the 
error in the velocity of signal propagation is introduced. The 
best accepted value of uncertainty in the speed of light (signal 
propagation constant) is +300 meters/second [13] . The uncom- 
pensated portion of the range measurement path for this error 
source can be as large as twice the Earth's radius with an un- 
compensated path length of 12,740 km corresponding to a nominal 
propagation time of 0.04 second. Therefore, the resulting 
range uncertainty due to speed of light error is +13 meters, 
which is considered to be a random error because of the proposed 
differential range measurement technique. 

3.1.3 System Noise and Measurement Errors 

3. 1.3.1 System Noise Effect 

System noise effect is that range error contribution 
attributed to the combination of received and self-generated 
input noise imposed with the ranging tones at the output of the 
aircraft receiver. There is no direct range error caused by 
input noise to the satellite transponder because phase com- 
parisons or measurements of the ranging tones are not made in 
this equipment. It is necessary, however, to consider the 
indirect error effect related to the retransmission of noise 
superimposed on the tone signals. The net result of total system 
noise received by the various communication links is to present 
a 30.5 dB-Hz signal-to-noise power density tone package to the 
aircraft receiver. Therefore, the worst-case signal-to-noise 
power density per tone becomes 24.5 dB-Hz; this value was used 
for the worst-case system noise error effect calculations 
developed in Section II-2.1.3 for a one-second signal integration 
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Figure 2-18 



period. For this worst-case signal-to-noise power density con- 
dition , the rms value of the noise-related ranging- error for the 

10.2 kHz tone is 200 meters. 

3. 1.3.2 Ranging Tone Processor Effect 

The fine grain ranging signal will be processed and 
measured digitally. The accuracy of digital processing can be 
an order of magnitude better than the corresponding analog 
implementation. The detail design of the processor is discussed 
in Section II-2.1; however, accuracy and processor requirements 
are presented in the following paragraphs. 

The general processor block diagram is shown in 
Figure 2-19. The input signal f(t) (consisting of signal and 
noise effects) is filtered and sampled using pairs of quadrature 
samples. Each sample value is converted to a digital word. The 
digital words are accumulated for the one-second ranging signals 
that are present and usable. 

To determine the accuracy of the processor it is 
necessary to consider the quantization effects of both the 
analog-to-digital converter and the accumulators . Analog-to- 
digital converter accuracy is typically stated in terms of the 
number of bits of the converter. An n-bit analog-to-digital 
converter is assumed to be accurate to within one-half the 
quantization level. (That is, the accuracy l/2 n+1 ) . 

The accumulator word size requirement is determined 
first. Range measurement errors for the range phase should not 
exceed 0.1 degree to assure a negligible error contribution due 
to the instrumentation inaccuracy. Noise effects are related 
to phase errors by the expression 


n 


(2E/N q ) 


T 7* 


radians 


,~3 


(2-71) 


Therefore, for a n = 0.1° or 1.75 x 10 radians, the energy- to- 
instrumentation noise ratio will be 55 dB. 


It may be shown that the required quantization 

interval (A ) in the accumulator is given by the expression 
s 


S o (t o } 


(2E/N q ) 


1/2 


(2-72) 


* 
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where S Q (t ) = P ea ^ signal amplitude that may be sampled 

2E/N = energy-to-noise power density ratio 
° (In this case, signal energy-to- 

instrumentation noise power density ratio) 

Then, it can be shown that A = 1/573. Therefore, the 

s 

number of bits (n) in each accumulator is 

l/2 n = A (2-73) 

O 

or n = 10 


Since the signal energy-to-noise power density ratios are 
smaller than the above (55 dB) instrumentation noise effect, the 
quantization level allowable for signal measurement is larger 
than that for the accumulator. 

The analog-to-digital converter (A to D) size can be 
determined by investigating the effect of the quantizer on the 
mean and variance of sampled signal. It is desired to choose 
an A to D converter size that will cause a minimum degradation 
of these parameters. Assuming a noise effect that is dis- 
tributed in a Gaussian manner, it can be shown that, the A to D 
converter effect will be negligible (less than .001 percent) for, 

A < o (2-74) 

where 

A is the A to D converter quantization interval 

a is the standard deviation of the noise. 

The minimum required energy-to-noise ratio of the 
fine-grain ranging tone is approximately 24.5 dB-Hz. The 
maximum energy-to-noise ratio may exceed this value by 10 dB 
which gives the smallest value Of a to be expected; this is the 
worst case and yields 



(2-75) 


2 

The value of a is the final measurement variance after samples 
m 

have been accumulated. The relationship between the variance 

2 

of the independent samples ( a ) and the final measurement 
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variance is given as 



°n/L 


(2-76) 


where L is the number of accumulated samples. 

To complete these considerations it is necessary to 
determine the number of independent samples (L) that can be 
accumulated during the one-second period. Since the minimum 
filter noise bandwidth will be 500 Hz the Nyquist sampling rate 
will exceed 1000 samples per second, or a minimum of 1000 samples 
for the one-second accumulation period. Therefore, 

2 _ 2 1 

m n 1000 (2-77) 


which is a processor signal-to-noise ratio improvement of 30 dB. 

A maximum energy-to-noise ra,tio of 4.5 dB-Hz is ex- 
pected at the A to D converter input. The corresponding value 

of a in radians is 
n 

— — radians (2-78) 

2.38 


(2-79) 


then 

2 n > 2.38 (2-80) 


n 


(2E/N q ) 


tti 


Since 


A = — 


c a. 

n — h 


which requires a value of n = 2. 

Therefore, for the signal-to-noise ratios, bandwidths , 
and total processor accuracy required a two-bit A to D con- 
verter used in conjunction with a 10-bit accumulator word is 
adequate to digitally process the fine grain ranging signal 
with a negligible error. In practice, the A to D converter will 
be larger and the accumulator word will be longer. This pro- 
vides increases performance margins and improved operational 
reliability. 

3.1.3. 3 Clock Timing Error 

Each of the three reference tones used at the aircraft 
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in measuring the phase of the received ranging tones is derived 
by dividing down the aircraft clock frequency. The ranging 
error caused by an unsynchronized clock is directly proportional 
to the accumulated phase difference between the aircraft and the 
ground station clocks . Clock errors result from long term 
oscillator drift, temperature instability, and initial frequency 
and phase offsets remaining uncorrected at the time of synch- 
ronization. For a crystal clock, the instantaneous frequency 

error (Aw ) at any time (t) following clock synchronization is 
c 

given by 


A "c “ “c £ LTS 1 + “c £ CAL (2 - 81) 

where 

w c = nominal clock frequency 

e LTS = lon< 3 term stability 

= Aw /w /24 hour period 
c c 

e = measure of the frequency offset after 
synchronization 

= Aw/w c at t = 0 

It is seen that the long term drift produces a clock frequency 
change which is a linear function of time. 

Equation (2-81) can be integrated with respect to (t) 
to determine the accumulated clock phase error. 

? A 0 n 

A9 c ' V2 s LTS t 2 + e CAL t + (2-82) 

c 

where 

A9 c = the instantaneous clock phase error 

A9„ = the initial clock phase error 
Oc c 

The ranging error is given by 

? A9 o 

AR = c ( 1/2 e LTS t + + ‘ ) (2-83) 

where 
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AR = the ranging error 

c = speed of light 

= the nominal ranging tone frequency 

A0 O = the initial phase error in the reference 
tone 

= A8_ ♦ a) /a) 

Oc r' c 

To minimize the effect of clock errors on the total system 
ranging accuracy, the allowable clock-related ranging error was 
set at 100 meters. 

It would be convenient if the update interval (the 
time between successive synchronizations of the aircraft clock 
with the ground clock) required to maintain the accuracy be 
longer than any flight duration that might take place in the 
area of coverage. This implies an update interval on the order 
of one day. 

Using a precision, temperature compensated, crystal 
oscillator with e LTS = 5xl0 -10 and e CAL = lxlO -10 and which can 
be set to within +0.1 degrees, the 100-meter limit would be ex- 
ceeded in slightly less than 50 minutes. For this time period 
the greatest source of phase error is the initial frequency 
offset. It is therefore apparent that some method must be used 
to compensate for clock errors since it is difficult to achieve 
much better drift rates and offsets. 

Compensation for frequency drift can be effected once 
this parameter has been accurately measured. In one technique 
[14] a positive frequency drift is partially compensated in the 
synchronization process by inserting a positive phase shift 
equal to the largest allowable error, and adjusting the clock 
frequency to obtain a negative frequency offset. By the time 
the positive frequency drift has exactly cancelled the negative 
frequency offset the aircraft clock lags the ground clock by the 
maximum allowable phase error. At this point, the positive 
frequency drift predominates and begins to advance the phase 
until it again leads the ground station clock by the maximum 
allowable phase error. When this occurs the clock must be re- 
synchronized. The time interval, during which the clock remains 
within the maximum phase error limit, is determined by the 
maximum allowable phase error and the long term stability. 
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By applying the technique to a clock with the same 
parameters listed in the previous paragraph, the, allowable range 
error would be exceeded after eight hours. Even though this 
technique provides considerable improvement over an uncompen- 
sated system it still does not meet the 24-hour criterion. 

An alternate technique is available , and does provide 
superior compensation for clock errors. Since an on-board 
computer is available, the accumulated phase error given by 
Equation (2—82) can be calculated and used to modify the 
measured phases of the received ranging tones. The only inputs 
required by the computer are the measured values of the initial 
frequency and phase offsets, the long term stability, and the 
elapsed time from synchronization. Since the accumulated phase 
error is removed before ranging calculations are made , the 
requirements for an extremely high long term stability and low 
frequency and phase offset are not nearly as critical as for the 
other techniques. The following clock specifications are con- 
sidered adequate and can be met by most general purpose pre- 
cision crystal oscillators, 

. . . 7 

long term stability < 3.0 parts in 10 /day 

8 

frequency offset < 1.0 part in 10 

inital phase offset <_ 0.1° 

In practice, the frequency standard would serve as 
the clock for the elapsed- time counter. The output of this 
counter would be available as an input to the on-board computer. 
The phase of the aircraft reference oscillator would be set by 
using phase measuring techniques similar to those used in the 
ranging tone processor, and adjusting the oscillator frequency 
up or down until zero phase is measured. The accuracy of the 
phase synchronization is the accuracy of the phase comparator 
in the ranging tone processor, which is designed to be +0.1 
degrees. 


The initial frequency offset would be measured after 
each synchronization for entry into the on-board computer. The 
frequency drift remains relatively constant and does not need to 
be measured each time the clock is synchronized. The implemen- 
tation of this technique is therefore straightforward and will 
maintain the clock-related ranging error over, a 24-hour period 
within the 100-meter limit without resorting to extremely high 
stability cloaks. 

3.1.4 System Delay Errors 

3.1.4. 1 Aircraft Terminal Uncompensated Delay Effect 

The effect of uncompensated time delay and phase shift 
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variations expected to occur in the aircraft terminal is an 
equivalent range error of 60 meters in the ranging path from 
the ground station-to-aircraf t . The resulting error is assumed 
to be random because the nominal values of time delays and 
phase shifts are to be removed by calibration. The balance of 
the range error effects caused by signal amplitude and frequency 
changes, temperature variations, input voltage variations, etc., 
contribute "random" errors. 

3. 1.4. 2 Satellite Transponder Uncompensated Delay Effect 

The nominal value of the satellite transponder group 
delay will be removed through simultaneous aircraft and refer- 
ence terminal range measurements. The residual or uncompensated 
time delay which is caused by variations in the signal amplitude 
and frequency is conservatively estimated to be 0.1 yseconds 
corresponding to a range error of 30 meters. 

3. 1.4. 3 Sector Interrogation 

The constant phase plane adjustment which is trans- 
mitted by the ground control center during a given time slot 
applies to only one sector. While the GCC is transmitting the 
constant phase plane correction for sector k, it is measuring 
the return from sector k+1 to determine the required correction 
to be transmitted in the next time slot. Because the GCC up- 
dates the constant phase plane of a sector with no feedback 
signal from that sector as the phase is being transmitted, an 
error is introduced as a result of satellite movement during 
the transmission. Section II-2.3 shows that this error can 
have a maximum value of 15m/second (calculated for ATS-III) . 

The duration of the ranging tone transmission for each sector 
is three seconds, implying a maximum error of 45 meters. This 
translates to an aircraft position error of approximately 60 
meters at 35° North Latitude. 

3.2 Error Sensitivity Analysis 

3.2.1 Differential Error Sensitivity Factors Affecting 

Position Location 

Error sensitivity factors are calculated for the dual 
satellite navigation configuration by taking the partial 
derivative of each of the aircraft-to-satellite range equations 
and taking the rms sum of the appropriate errors. The following 
equations express the measured range from the aircraft to the 
satellite in terms of all the delays appearing in each segment 
of the ranging tone navigation link. The parameters are as 
shown in Figure 2-20. 
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The following equations describe the path delays used in this 
analysis. 


1 — 1 
aS 


2(d 1 + d 3 + d sl + 

d + d . + d ) 
a ml c 

(2-84) 

R rl 

= 

2d 1 + 2d 2 + d sl + 

d~ + d + d 
s2 r g 

(2-85) 

R a2 

= 

2 (d 6 + d s3 + d 4 + 

d a + d m3 + d c> 

(2-86) 

R r2 

= 

2d 6 + 2d 5 + d s3 + 

d s4 + d r + d g 

(2-87) 

d 3 

= 

i[R _-R . +2d 0 +d ~+d +d -2d -2d , -d„, -2d ] 
2 al rl 2 s2 r g a ml si c 

(2-88) 

d 4 

= 

■i[R _-R 0 +2d c +d .+d +d -2d -2d 0 -d„--2d ] 
2 a2 r2 5 s4 r g a m3 s3 c 

(2-89) 

where 



• 


d l 

= 

range from ground 

station to satellite 1 


d 2 

= 

range from satellite 1 to reference terminal 
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= derived range from satellite 1 to aircraft 

d^ = derived range from satellite 2 to aircraft 

d 5 = range from satellite 2 to reference terminal 

dg = range from satellite 2 to ground station 

R , = twice the one-way range measurement from GCC to 

aircraft via satellite 1 

R 1 = two-way range measurement from GCC to reference 

r terminal via satellite 1 

R - = twice the one-way range measurement from GCC to 

a aircraft via satellite 2 

R r2 = two-way range measurement from GCC to reference 
terminal via satellite 2 

d a = range delay effect of aircraft terminal 

d r = range delay effect of reference terminal 

d = range delay effect caused by on-board clock error 

d, = forward signal range delay effect of satellite 1 
transponder for R a measurement 

d gl ' = forward signal range delay effect of satellite 1 
transponder for R measurement 

d g2 = return signal range delay effect of satellite 1 
transponder for R & measurement 

d g2 ' = return signal range delay effect of satellite 1 
transponder for R r measurement 

d g3 = return signal range delay effect of satellite 2 
transponder for R a measurement 

d g3 , = return signal range delay effect of satellite 2 
transponder for R r measurement 

d . = forward signal range delay effect of satellite 2 

transponder for R a measurement 

d g ^' = forward signal range delay effect of satellite 2 
transponder for R r measurment 

d = range delay effect of ground processor for R a 
g measurements 
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d ' = range delay effect of ground processor for R 

^ measurements r 

d m ^ = range error effect due to multipath in satellite 1 
to aircraft link 

d^ = range error effect due to multipath in satellite 2 
to aircraft link 


From the above equations the sensitivity factors re- 
lating satellite-to-aircraft range error to the various 
parameters of uncertainty are calculated. The results are 
tabulated in Table 2-5. The calculations used to determine the 
sensitivity factors are made according to the following analysis: 


Table 2-5. Derived Satellite-to-Aircraft 

Range Error Sensitivity Factors 


Error Source 


Sensitivity Factor 

System Noise 


0.707 

On-board Clock Error 


1.0 

Terminal Delay 


0.707 

Ground Processor Delay 

0.707 

Satellite Transponder 

Delay 

1.0 

Along Range Satellite 

Position 

0.008 

Cross Range Satellite 

Position 

0.13 

Atmospheric Effects 


1.41 

Satellite-to-Aircraft 

Multipath Effect 

0.5 


NOTE: Delay uncertainty for differential configurations in- 

clude only short-term variation of delay that occurs 
during time intervals of 1 minute . 


From Equation (2-88) , sources of measurement error are 

R , and R , . These two sources are assumed to cause equal and 
al rl 

independent errors. Therefore, the total error sensitivity 
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factor is the rms sum of the factors for the two sources. The 
error sensitivity of to is 


3d 

Tr 


3_ 

al 


1/2 


and for the sensitivity is 



1/2 


(2-90) 


(2-91) 


Assuming independent range measurement errors, the error 
sensitivity to range measurement is given by 




0.707 


(2-92) 


Therefore, the error sensitivity to range measurement is 0.707 
which is the first entry to the Table 2-5. 


In a similar manner the error sensitivity factor due to 
terminal delay uncertainty is calculated to be 0.707. This 
value is based on the assumption that the uncertainties in both 
aircraft and reference terminal delay (d and d ) generate 

cl 3T 

errors which are independent and of Gaussian distribution. 


The range error sensitivity to variations of the ground 
processor delay occurring between the aircraft and the reference 
terminal measurements is given by 



(2-93) 


The sensitivity to satellite delay involves the sum of four 
independent delay factors corresponding to four one-way passes 
through the satellite transponder. Hence, for equal value of 
uncertainty for all four passes the range sensitivity factor 
to satellite delay becomes 
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The calculation of range error sensitivity to satellite 
position error for the differential navigation configuration is 
complicated by the fact that this error sensitivity is not con- 
stant over the area of coverage, but rather is a function of the 
geometry between the aircraft, reference terminal, and satellite. 
In particular, this error sensitivity is shown in the following 
paragraphs to be proportional to the angle formed by lines join- 
ing the satellite to the positions of both the aircraft and the 
reference terminal. 

Satellite position uncertainty is defined to be the 
difference between the assumed position and the actual position 
of each satellite. This position uncertainty is then broken 
into two components, the first one being the projection of 
satellite position error along the satellite-to-aircraft range 
line, and the second one being the projection of satellite 
position error across the satellite-to-aircraft range line. 
Assuming differential range measurements, the geometric relation- 
ship showing the residual or uncompensated satellite-to-aircraft 
derived range error for the along-range projection of satellite 
position is shown in Figure 2-21. The resulting derived range 
sensitivity factor to the along-range projection is expressed 
by the following equation 



F, 

sa 

[2 sin (a/2)] 

(2-95) 


= 

0.008 


for 

a = 

o 

Lf) 

• 

f" 


and where 

F_ is 
sa 

the sensitivity 

factor between derived range 


and satellite position uncertainty in the along-range projection, 
and a is the worst-case differential angle at the satellite 
between the reference and aircraft terminals. 

The geometry relating the effect of satellite position 
uncertainty in the cross-range direction to satellite-to- 
aircraft range error is shown in Figure 2-22. The effect of 
satellite position uncertainty in this direction is not can- 
celled by the differential navigation configuration and the 
sensitivity factor is expressed by the following equation. 

F so = sin a 

= 0.13 for a = 7.5° (2-96) 
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-21. Differential Range Error Geometry for Along- Range Position Uncertainty 
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Figure 2-22. Differential Range Error Geometry for Cross-Range Position Uncertainty 



* 


where F is the sensitivity factor relating aircraft-to- 
s c 

satellite range error to satellite cross-range position error, 
and a is the differential angle from the satellite between the 
aircraft and reference terminals. 

The sensitivity factor for the atmospheric effect is 
defined to include the range error due to propagation un- 
certainty through both the ionosphere and the troposphere. The 
listed sensitivity factor of 1.41 is based upon the assumption 
that independent errors due to atmospheric effects occur in 
both the satellite-to-aircraf t and satellite-to-reference 
terminal range paths (d^ and d 2 for satellite 1) . The listed 
value for this sensitivity factor does not take into account 
common mode rejection of the atmospheric bias error which is 
achieved in the differential navigation configuration. This 
rejection is expected to significantly reduce the sensitivity 
factor when the separation distance between the reference and 
aircraft terminal is less than 400 kilometers. 


3.2.2 Position Location Accuracy 


The technique used to derive the aircraft position 
location accuracy is summarized by the following equations ex- 
pressing standard deviations. The rms range error (a n ) in the 
aircraft-to-satellite 1 link is related to range error 
sensitivity factors (F.^) and the individual rms parameter 
errors (ck) by the expression 


a 


R1 



UF U ) 





1/2 


(2-97) 


The expression for the range error (o^) 
aircraft-to-satellite 2 link is 


o 


R2 



[(F 2i ) 2 (o.) 2 


1/2 


(2-98) 
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The aircraft latitude position error is related 

to the latitude geometric dilution of precision factors for both 
satellites ( G 1LA p and G 2LAT ) by 


LAT 


I 1/2 


(a Rl ) (G 1LAT ) + ( °R2 ) (G 2LAT J 


(2-99) 


The aircraft longitude position error ( a L0 N^ * s given 
by a similar expression which is 


LON 


-I 1/2 


(0ri) + ^°R2 ) (G 2L0N ) 


( 2 - 100 ) 


Finally, the aircraft's total position error (<*po S ) i- s 
the rms sum of the errors in latitude and longitude. This 
expression is 

1/2 


POS 


(a LAT } + ^LON 5 


(2-101) 


The next required step in the determination of the 
relative position accuracies is to establish an error allocation 
for the various parameters that are involved in the fine-range 
measurement. Table 2-6 summarizes the error allocation for the 
10.2-kHz ranging tone. The allocated projections of satellite 
position error with a constant phase plane adjustment and 
approximate ephemeris data can be calculated from the following 
typical values of satellite position parameters. 

Satellite position error in latitude = 

Satellite position error in longitude = 

Satellite position error in altitude = 

Satellite differential angle between 
earth's center and aircraft = 

(maximum for a 10-degree satellite 
elevation angle) . 

The rms value of errors in the derived satellite-to- 
aircraft range were calculated by multiplying the individual 
sensitivity factors by the respective allocated errors. The 
results are tabulated in Table 2-7. 


2 km 
1 km 
0.3 km 

8.5 degrees 
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Table 2-6. Error Allocation for Basic Range 
Measurements (10.2-kHz tone) 


Error Source 

RMS Error (meters) 

Range Measurement Error caused by Noise 
(Tone E/N q = 24.5 dB-Hz) 

200 

On-board Clock 

(Phase error = 1.2°@10.2 kHz) 

100 

Terminal Delay Uncertainty 
(Delay Error = 0.4 ysec) 

60 

Ground Processor Delay Uncertainty 
(Delay Error = 0.1 ysec) 

30 

Satellite Transponder Delay Uncertainty 
Delay Error = 0.1 ysec) 

30 

Satellite Position uncertainty 
Along-Range Projection 
Cross-Range Projection 

445 

2200 

Ionospheric Effects 

(Satellite elevation angle = 10°) 

60 

Satellite-to-aircraft Link Multipath Effect 
(Phase Error = 1°) 

90 


Table 2-7. RMS Errors in the Derived Satellite- 
to-Aircraft Range 


Error Source 

RMS Error (meters) 

System Noise 

141 

On-board Clock 

100 

Te rmin al De lay 

43 

Ground Processor Delay 

22 

Satellite Transponder Delay 

30 

Along-Range Satellite Position 

4 

Cross-Range Satellite Position 

290 

Ionospheric Effects 

85 

Satellite to A/C Multipath Effect 

45 

Total rms Error 

354 
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SECTION III 


EXPERIMENTAL SYSTEM DESIGN 
1. INTRODUCTION 

The preceding design study has shown that a system can 
be realized which will determine position location at a remote 
point on earth's surface to within 1500 meters by measuring a 
set of ground originating signals emanating from two synchron- 
ous satellites. The system can be used aboard an aircraft if 
altitude inforamtion is provided . The baseline system 
proposed herein will demonstrate the capabilities of an 
experimental system utilizing existing equipment to a maximum, 
and will 'illustrate extendability to a more sophisticated system 
for operational status. 

The system will be designed around the ATS-I and ATS-III 
VHF synchronous satellites using the ground control center 
(GCC) incorporated in the VHF Satellite Navigation Experiment. 
Because of these two considerations, some system parameters are 
pre-determined which result in constraints on the design. 

These constraints must be used as the basis for operational 
analysis and for performance prediction. 

Of primary consideration is the satellite system made 
up of ATS-I, ATS-III, and the Earth, which has the following 
major characteristics: 

(1) The satellites are synchronous but due to orbital 
inclination, their positions vary with time. 

(2) Satellite operation frequencies are: 

Transmit - 135.6 MHz 
Receive - 149.22 MHz 
Bandwidth - 100 kHz 

(3) Satellite transmit power is fixed, which limits 
the signal-to-noise ratio at the Earth's surface. 

Implications of these characteristics are quite broad and will 
be covered as they fit in context. 

The addition of the GCC to the system further reduces 
the number of design choices and limits system performance 
because : 

(1) The GCC-generated ranging tones are 426 Hz, 707 Hz, 
and 941 Hz. 
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( 2 ) 


The GCC is limited in phase measurement accuracy. 




(3) The VHF link parameters limit the signal-to-noise 
performance of the system, as discussed in detail 
in Appendix A. 

The preceding factors illustrate that the optimum 
navigation system discussed in Section II cannot be realized 
with the ATS-I , ATS-III, and GCC as now configured. Neverthe- 
less, these components are sufficient to demonstrate the 
principle and feasibility of on-board navigation by ranging to 
two satellites, and the following base-line system is proposed to 
effect this demonstration. 

2 . BASELINE SYSTEM 

Certain operational aspects of the navigation system 
must be selected prior to describing the general electronic 
design. These aspects are discussed in the following paragraphs. 

The study has shown that a simple phase plane correction 
cannot completely compensate for the satellite position 
uncertainty; however, a combination of two techniques can com- 
plement the phase plane correction to attain adequate accuracy. 
First, satellite motion can be partially compensated for by a 
very simple approximate equation. Second, the geometric errors 
can be reduced by placing reference terminals at the center of 
geographic sectors to time share the GCC. The first technique 
requires a software function to be implemented in the GCC com- 
puter and in the remote aircraft terminal computer. (Throughout 
this discussion the computer and other equipment at the navi- 
gation site will be termed aircraft (A/C) terminal) . The 
second technique requires a selection of an appropriate number 
of sectors based upon the end requirement of demonstrating 
adequate navigational capability and extension to operational 
status. In an operational system, the number of sectors would 
be determined by a different criteria such as satellite coverage, 
maximum allowable time between position fixes, minimum time to 
obtain an accurate fix, and others. For this experiment, it 
is expected that three sectors would be adequate to perform 
the various tests. 

Figure 3-1 is a diagram of the resulting system com- 
ponents. The multiplicity of RF links is apparent, and for 
simplicity throughout the discussion the following representa- 
tion will be used to designate a link: 

ATS-I to A/C 

or (Transmitter) to (Receiver) 
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Figure 3-1. Major Elements Of A Dual Satellite Navigation System 




That is, the transmitter will be the first element mentioned. J 

Channelization is an additional aid in the description of the 

various link frequency spectra. The satellites are entered 

(receive) at 149.22 MHz and transmit at 135.6 MHz band center 

frequencies. For convenience their 100-kHz bandwidths are * 

divided into 40 equal 2.5-kHz intervals or channels, numbered 

from 1 to 40 beginning with 1 as the lowest frequency 2.5-kHz 

channel. This is illustrated in. Figure 3-2. Figure 3-3 then 

shows the channels used and/or occupied in each of the links. 

Six major implications are embodied in Figure 3-3. 

(1) The channels selected are grouped within a 50-kHz 
bandwidth so that the GCC may be operated with 
its narrow (50 kHz) bandwidth receiver IF filters. 

(2) The GCC antenna gains are approximately 20 db 
down in the direction of the opposite satellite 
such that any radiation in, or reception from, 
that direction is negligible. This is represented 
by the figure in that no GCC-to-ATS-I channel 14 
appears at ATS-III, etc. 

(3) The antennas at the reference terminals are 
assumed narrow beam even though the frequency 
multiplexing would conceptually allow illumination 
of both satellites. The directivity is specified 
to improve the equivalent radiated power of the 
link and to prevent unnecessary power sharing 

at the satellites. 

(4) The selected channels, are arbitrary and are chosen 
as shown only to give a reference point for 
discussion. They may be relocated any time up to 
final design of the hardware, within the 
constraint of the 50-kHz GCC receiver IF bandwidth. 

(5) A pilot tone will be used with ATS-III to allow 
tracking out the jitter known to be present in 
the transponder. 

(6) The tone package is transmitted from the GCC to < 

the satellites (s ) , then to the A/C and reference 

terminal, from the reference terminal back to the 
satellite (s) and finally back to the GCC. 

Reviewing briefly, the baseline system is to consist of 
a GCC which controls the system operation by applying phase 
correction to the ranging tones such that the phase received at 
the A/C can be used to determine position with minimal error. 
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In addition, three reference terminals to return the signals to 
the GCC for measurement and subsequent phase correction cali- 
bration, and an A/C terminal to effect the phase measurements 
and position determination while airborne constitute the total 
system component complement. 

3. SYSTEM OPERATION AND CONTROL 

3.1 Cycle Description 

3.1.1 General 

The experimental system consisting of the components 
stated in the preceding section operates in a cyclic fashion. 

This arrangement was shown diagrammatically in Figure 3-1 and 
can be referred to here for clarity. The primary function of the 
GCC is to measure phase perturbations caused by satellite 
movement and/or atmospheric diurnal variations. From this 
measurement the GCC then calculates a transmit phase such 
that the perturbation effects are minimized at the center 
of the sector; i.e. , at the reference terminal. The overall 
result is that phase compensation at the GCC corrects for the 
phase errors caused by these perturbations that otherwise would 
occur in the navigational sector. 

A simplified block diagram of the GCC system for one 
tone is illustrated in Figure 3-4. As shown, the transmitted 
phase is directly controlled by the computer. The signal 
returned from the reference terminal is coherently compared to 
the original master clock and digitized so that the phase shift 
over the link can be calculated by the computer. It is 
corrected (or biased) for ephemeris, added as a sector offset 
constant, and stored for later transmission. 

3.1.2 Sequence and Control 

System sequence and control will be synchronized by 
stable clocks which control the time division multiplexing 
among the sectors. During the time a phase corrected 
navigation tone package is being transmitted to a sector, the 
reference terminal in the next sector in sequence will 
transpond the signal. This is done so that the GCC can compute 
the applicable phase shift as near as possible to the time it 
will be applied. Figure 3-5 shows the operating time relation- 
ships for the three reference terminals. As the GCC transmits 
the corrected phase signals for sector 3 the signal is 
transponded from the sector 1 reference terminal. In the next 
time slot the corrected phase as just measured for sector 1 
is transmitted and transponded from sector 2. Next, sector 2 
corrected phase is transmitted and sector 3 transponds , with 
the following step completing the cycle. 
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Figure 3-4. GCC Functional Diagram 




















3.1. 3 


Sector Period 
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The duration of transmission to each sector has an 
upper limit dependent upon a number of factors : 

(1) Number of sectors 

(2) Time between position fixes 

(3) Phase measurement time required 

(4) Circuitry settling time 

(5) Transmission time guard band 

The number of sectors limits the duration only if the number 
is so large that acceptable time between position fixes 
would be exceeded. Assuming that the maximum allowable time 
is three minutes between required position updates [15J, then 
10 sectors would allow 18 seconds of tone transmission per 
sector, or 50 sectors would allow 3.6 seconds per sector. 
Obviously for the experiment with only 3 sectors no time 
limitation problem will exist. For an operating navigation 
system the number of -sectors will depend upon allowable phase 
error at the sector edge and the total area of coverage. Here 
it is seen that the experiment results will define the 
sectorization parameters, and therefore will provide a solid 
base for trade-off decisions. 

The amount of time required to integrate for noise 
error reduction is plotted in Figure 3-6. The point (X) shows 
that for a carrier-to-noise power density of 30 db-Hz, with 
a two-second integration period after the detector, the rms 
phase error would be 0.9 degree. The experimental system is 
expected to have closer to 40dB-Hz C/N q for which two seconds 
of integration reduces the rms phase error to 0.3 degree. 

The integration period will be selected so that noise error 
does not dominate, thereby assuring other less controllable 
factors are the system limitations. 

The circuit response times fall into three categories: 
(1) time of transmission, (2) settling time, and (3) acquis- 
ition time. Differences in transmission time from sector to 
sector and settling times for the 40-Hz bandwidth filters are 
on the order of tens of milli-seconds , whereas the phase- 
locked loop acquisition cycle is four seconds. For the three - 
sector system this could be tolerated. It is planned, however, 
to subsitute new fast acquisition phase-lock-loop receivers 
developed by Texas Instruments to verify system operation 
timing. 
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Finally, the guard band time is selected to accommodate 
the differences in satellite-to-sector distances, the maximum 
change in satellite apparent range, and the sector timing 
inaccuracy. From the above considerations, a sector period of 
five seconds has been selected for the experimental system. 

3.1,4 Navigation Signal 

From the preceding section, it is seen that for one 
third of the system cycle a phase-corrected tone package will be 
available at each sector for determination of position.. During 
the other two-thirds of the cycle a tone package corrected for 
a distant sector will be available for range computation. The 
range calculated during this portion of the cycle will be in 
greater error than would be attained with the proper sector 
signal. 


The frequency spectrum of the tone package prior to 
the GCC transmission is as shown in Figure 3-7. This spectrum 
is identical to the OCC self-check subsystem set of tones. 


f t = 149.22125 MHz 
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Figure 3-7 


GCC To ATS-III Ranging Tone Spectrum 



Conversion to computer controlled phase is, therefore , ‘ simply 
a replacement of the OCC Self-Check Amplitude and Interval 
Selector circuit with a Controlled Phase Oscillator circuit. 

A block diagram of the Self-Check Subsystem is shown in Figure 
2-19 of the OPLE Maintenance Manual [16], m which the Amplitude 
and Interval Selector is shown as one printed wiring board. 

The conversion to computer controlled phase would simply re- 
quire removal of this card and connection of the three-tone 
Controlled Phase Oscillator, as shown in Figure 3-8, in its 
place. The resulting circuitry will be sufficient to generate 
the tones for one satellite. To generate the tones needed for 
the second satellite the circuits described will be duplicated, 
except where signal or function sharing is possible. 

3.1.5 On-Board Measurement 

The signals as received at the A/C during the assigned 
time interval for the sector in which it is situated will be 
corrected such that the on-board assumed satellite position 
will provide the best range calculations. To assure reception 
of both satellites the A/C antenna will be broadbeam with the 
following electronics split into two channels. One channel 
will accept the ATS-I signals and the other the ATS-III signals. 
The twelve sine and cosine function measurements (3 tones, 2 
satellites) which will be the system outputs will then be 
digitized and used by the computer along with altitude obtained 
from other on-board instrumentation to compute position. The 
resulting position will then be displayed. 

3.2 Ground Control Center (GCC) Functional Description 

Fundamental design of the GCC is the OPLE Control Center 
(OCC) as modified for the VHF Satellite Navigation Experiment, 
and as shown in simplified block diagram form in Figure 3-9. 
Figure 3-10 is the same block diagram with the additions for 
the Dual Satellite Navigation Experiment super-imposed. These 
diagrams can be compared to most easily view the system 
additions. The predominant changes are the additional transmit 
channel and- the controlled phase oscillators with the computer 
control. 

The modifications result in a system consisting of two 
sets of VHF transmit-receive equipment controlled by a single 
computer with one side assigned to the ATS-I satellite and the 
other to ATS-III. The primary function of the GCC is phase 
control of the ranging tones. Since both sets operate in 
identical fashion, discussion of one side is sufficient to des- 
cribe both. The ultimate output is a time multiplexed phase- 
adjusted tone package, which is transmitted to the satellite so 
that it can be relayed to and for use in the sector to which 
the time slot is allotted. The total sequence from the 
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Figure 3-8. Three-Tone Controlled Phase Oscillator 
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Figure 3-9. VHF Satellite Navigation Experiment GCC Simplified Block Diagram 





















beginning proceeds in the following manner: The computer senses 

the sector clock output and selects from memory the tone phase 
information stored there. After turning the system on, the phase 
in memory for the first sector would be zero, resulting in an 
uncorrected navigational signal. The phase from memory is 
transferred to the controlled phase oscillators which generate 
squarewaves at the tone frequencies and with proper phase shift. 
The squarewaves are inputs to the self-check subsystem which 
reject the higher order harmonics, leaving sine wave tones. In 
the -check subsystem the tones are also translated up, 

summed with the A/R tone, and translated to the 500-kHz inter- 
mediate frequency. The next signals are fed into the VHF 
transmitter where the tone package is translated to the 
149.2-MHz transmit frequency and placed at the center of the 
assigned channel. This signal is routed through the diplexer 
to the antenna for transmission to the satellite. The satellite 
then retransmits the signal, at 135.6 MHz into the entire area 
of coverage. The satellite-relayed signals are utilized in two 
different processes. In one, users in the sector to which the 
transmission is intended make their phase measurements and 
calculate position. In the other, the reference terminal in 
the next sector in the sequence transponds the signal, shifting 
it to a different channel, and transmits back to the satellite 
at 149.2 MHz. The satellite transponded signal then arrives 
back at the GCC at 135.6 MHz in a different channel than 
originally transmitted. The receive portion of the GCC down- 
translates the tones and samples them with quadrature signals 
in exactly the same way that the OCC receivers now operate. 

The computer then calculates the received phases, compares them 
to the transmitted phases , and calculates the new phases to be 
transmitted to the next sector. The phase data is stored in 
memory to be applied at the next signal from the sector clock. 
Succeeding sectors are handled in the same manner. 

3.3 Reference Terminal Functional Description 

In order to correct the signal phase a reference terminal 
transponder is placed mid-sector to return the tone package to 
the GCC. Since the geographic position of the reference 
terminal is known the phases received at the GCC describe the 
range perturbations due to satellite movement and atmospheric 
change. These may then be compensated by altering the trans- 
mitted phase. The reference terminal requirements for this 
function are simply to transpond the ranging tones during the 
sector time allotted with minimal phase shift of the tones. 

The block diagram of the reference terminal shown in 
Figure 3-11 is precisely the VHF Satellite Navigation 
Experiment reference terminal with the addition of the clock to 
control time of transmission. Note that the phase-lock-loop 
receives the A/R tone (and ranging tones) at all times, and the 
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multiplexing is carried out on the transmit side. This 
effectively eliminates consideration of acquisition time. 

3.4 Aircraft Terminal Functional Description 

In the system being described, the term Aircraft 
Terminal means much more than a translational transponder, as 
used in the VHF Satellite Navigation Experiment. The A/C 
terminal becomes the phase measurement, the range calculation, 
and the position computation equipment, implying of course, 
that the A/C need not retransmit the signal for ground 
measurement and calculation. Figure 3-12 is a simplified 
block diagram of the A/C terminal. There are two receiver 
sides, one to accept ATS-I signals and the other for ATS-III 
signals. The VHF broadbeam antenna receives both satellite 
signals which are pre-amplified and power-split to drive the 
two sides. Considering one side, the signal is down -converted 
to IF in the receiver and phase-lock converted further in the 
following stage. The resultant baseband tones are sampled by 
in-phase and quadrature signals with the sampled components 
feeding an integrator. The integrated samples are then 
digitized and input to the computer for phase and range cal- 
culation. The altitude information for the third sphere of 
position is multiplexed into the computer along with the phase 
detector outputs. Calculated position is then displayed as 
latitude and longitude. 

Compensation for accumulated on-board clock phase errors 
will be accomplished by implementing the technique described in 
Section II-3.1. During the experiment the clock will be synch- 
ronized on the ground every 24 hours. The frequency offset 
will be measured and stored in the computer along with the 
frequency drift rate. The elapsed time indicator, which will 
be an integral part of the clock, will be reset each time the 
clock is synchronized. 
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Figure 3-12. Aircraft Terminal Simplified Block Diagram 


















SECTION IV 


EQUIPMENT RECOMMENDATIONS 


1. INTRODUCTION 

The equipment recommendations in the following paragraphs 
reflect an effort to use as much as possible the existing hard- 
ware and design that was developed for the OPLE Experiment and 
the VHF Satellite Navigation Experiment. Primarv emphasis in 
the discussion, however, is given to the new key designs unique 
to the Dual Satellite Navigation system requirements. The 
three equipment groups - GCC , Reference Terminal, and A/C 
Terminal - are discussed in separate paragraphs followed by a 
concluding paragraph in which the expected system performance 
is summarized. 

2. GCC EQUIPMENT 

2.1 Unmodified Equipment 

Much of the GCC as it now exists, unmodified, will be 
used. The computer and peripheral devices; i.e., ADC, magnetic 
tape recorder, and teletype, have the capacity and capability 
to perform the functions required for the system. New pro- 
gramming is the only requisite to bring the computer system up 
to operating status. One statement of reservation is in order 
at this point: The computer system has had problems with 

reliable operation in the past. As a result the experiment 
must allow schedule and cost for the purpose of computer main- 
tenance and to cover reruns if a problem occurs mid-experiment. 

The entire RF sections now in use will be employed as is. 
These include all parts starting at the two VHF antennas and 
positioners, through to the phase-locked converters. There are 
needed components missing, however, these will be discussed 
later. 


The digital control electronics also will be used with- 
out modification, but a large portion of it will be idle as is 
the case in the VHF Satellite Navigation Experiment. The 
operation of this subsystem can be set up by the control signals 
and needs little special effort. 

Various portions of the low frequency analog electronics 
are also to be used. In particular, the self-check subsystem 
tone filters and associated converters will be implemented to 
condition the ranging tone package prior to transmission. The 
other analog circuits to be used are simple buffer amplifiers, 



filters, shaping circuits, etc., which interface the parts of 
the system. 

Finally, the computer-addressed buffer registers which 
constitute the data flow electronics will be reconnected where 
needed to permit their usage in transferring the phase infor- 
mation to the proper oscillator. 

2.2 Modified Equipment 

Two areas of the GCC require modification to permit use 
in this system, the most important of which is the phase 
measurement subsystem. Recent measurements in the VHF Satel- 
lite Navigation Experiment have indicated that the rms error of 
the present circuitry is approximately 1.5 degrees, which at a 
maximum tone frequency of 941 Hz corresponds to 1.3 km rms 
range error. This error magnitude is considered excessive. 

The linear integrated circuit-operational amplifier technology 
available today makes possible improvement in the circuit 
functions. Figure 4-1 is a schematic of the present phase 
measurement circuit for one tone. The four major circuit func- 
tions represented include: 

(1) Passive Bandpass Filter 

(2) Symmetric Limiter 

(3) Active Bandpass Filter 

(4) Quadrature Sample, Integrate, and Dump 

The latter three of these can definitely be improved. The 
symmetric limiter is easily replaced with a high-gain comparator 
operational amplifier followed by low-offset clamp switches to 
accurate plus and minus voltage references. Such a limiter 
ensures clipped signals at minimum signal-to-noise ratio and 
suppresses second harmonic components. This simplifies the 
active bandpass filter and prevents wide amplitude deviations 
from appearing at the sampler input. 

The active bandpass filter can be synthesized using 
operational amplifiers with high feedback ratios to insure 
gain stability. This results in an improved transfer function 
stability, controlability , and accuracy. 

Last, the quadrature (sine - cosine) sampling and the 
integrate and dump functions can be improved. The dynamic 
range and impedance requirements for the series modulating FET 
switch (see Figure 4-1, Q g and Q g ) limit the sampling accuracy 
possible with this method? An improved circuit results if the 
FET is enclosed in the feedback loop of a zero-order hold 
circuit. The sample and hold function, which this creates, 
will give less switching error, and after integration is the 
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same as the gate and integrate for one-half period intervals 
which describes the present circuit. In addition, the passive 
RC integration can be improved by replacing it with a capacitor 
feedback operational amplifier Miller integrator. Thus, the 
RC time constant can be decreased so that the 1-megohm resistor 
can be reduced to a more reasonable value. 

2.3 Duplicated Equipment 

The GCC as it is in operation in the VHF Satellite Naviga- 
tion Experiment does not have the capability of transmitting to 
two satellites simultaneously. Conversely, the receiving ele- 
ments are all in existence including two VHF antennas. In order 
to supply the second transmit capacity the diplexer, transmitter 
subsystem, and self-check subsystem will be duplicated. 

The diplexer, frequency synthesizer, and VHF power ampli- 
fiers in the GCC are purchased items and will be duplicated by 
buying equivalent units. The rest of the transmitter subsystem 
is circuitry which can be component for component reproduced 
from drawings that are on file. These circuits will be tested 
and adjusted to procedures that are proven, and interconnected 
with confidence in wire list accuracy. There will be some 
design effort required to determine what circuits can be deleted 
and to compensate in loading, etc., for their absence; however, 
the overall duplication effort is minimal. 

The necessary effort to reproduce the required portions 
of the self-check subsystem is similar to the effort required 
for the transmitter subsystem as described above. An exception 
is that no major element need be purchased. 

2.4 Additional Equipment 

The preceding paragraphs of this section have dealt with 
use of the GCC essentially as it is. There are, however, two 
new additional requirements. These are the computer controlled 
variable phase circuits (Phase Shift Subsystem) which control 
the transmitted tone phase, and a totally new set of computer 
programs . 

2.4.1 Phase Shift Subsystem 

The Phase Shift Subsystem will interface with the GCC at 
the control electronics data output register as shown in Figure 
4-2, which for simplification shows the block diagram for one 
tone only. Total complexity would include two more circuits 
which differ only in the divider ratios and the connection to 
different data output registers. As shown in the block diagram, 
the data from the computer is transferred (under software 
control) into the data output register. After a short delay a 
transfer pulse is generated in the gate circuit to transfer the 
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Figure 4-2. Block Diagram of Phase Adjust Subsystem for One Tone 
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data to either the delay register for ATS I or ATS-III delay as 
determined by the state of bit EB15, which is also computer 
controlled. This transfer also gates the clock off at the 
divider for the Phase Shifted signal (f /4>) and sets the 
divider to state zero. On the next positive transition of the 
reference frequency, at zero phase, the delay register begins 
to count down. When it reaches zero the divider clock is 
enabled and the phase shifted signal commences. Thus, in less 
than 2 cycles of the tone its phase has been shifted by an 
integral number of periods of the master clock. The clock 
frequency of 3.946308 MHz has been selected to give the frequen- 
cies and minimum phase shift increments as shown in Figure 4-3. 

Figure 4-4 and 4-5 are logic diagrams of the period 
generators [17] for the 941-Hz and 707-Hz signals. The 426-Hz 
divider is similar, but requires three more circuit packages. 

Each divider is used twice, first to generate the tone references 
in quadrature to drive the sample gates, and second, to be preset 
and gated to generate the shifted tones which are filtered and 
shifted for transmission by the self-check subsystem. Figures 

4- 6 and 4-7 represent the logic diagrams of the Delay Registers 
and the Delay and Gate circuit. Operation of the period gen- 
erator dividers is described in the referenced report [17]. The 
delay registers are 15-stage shift registers made up of three 

5- bit binary counters, with decoding of the zero state used to 
disable the clock and to enable the divider. Finally, the Delay 
and Gate senses the EB15 and the DRNX gated with the address 
signals from the GCC to produce the pulses that transfer the data 
to the proper register, ATS-I or ATS-III, and clears the divider 
to the beginning state. 

2.4.2 Software 

In Section II. 3. 2 it was shown that the system 
sequence is under computer control, and since this operation is 
different than either OPLE or the VHF Satellite Navigation 
Experiment, new computer software is required for this system. 

The GCC general software organization has been covered in the 
"Design Study Report for Omega Position Location Equipment 
Control Center" [18] and is directly applicable to this system. 

Computer control is under an "executive" program 
that waits for an interrupt signal to initiate a subroutine. The 
executive program is described in the aforementioned document; 
therefore, only the needed new subroutines will be discussed 
here . 

Sequence of operation is triggered by a clock signal 
on one of the interrupt lines which indicates either time to 
output new phase data or time to read and process measured phase. 
The two subroutines needed then are the Signal Phase Measurement 
and the Phase Data Output subroutines which will be described 
below. 
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Figure 4-3. Ranging Tones Dividers 
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Figure 4-4. 941-Hz Divider 
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Figure 4-5. 707-Hz Divider 
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Figure 4-6. Delay Register 
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Figure 4-7. Delay and Gate Circuit 



Signal Phase Measurement (SPM) 


At the end of the phase measurement integration period 
the digitized signals are transferred into a buffer register and 
the interrupt line associated with the sector being received is 
pulsed. The phase data (sine and consine) are reduced to a phase 
angle, which is used with the transmitted phase to compute the 
phase at the reference terminal. From this information, and a 
calculation of the assumed satellite position and range, a new 
transmit phase for each tone is produced. This data is converted 
to number of cycles of the master clock, keyed in the 16th bit 
as to application for ATS-I or ATS-III, and stored in the appro- 
priate memory location. 

Phase Data Output (PDO ) 

The beginning of each sector transmission will be initiated 
by one of three interrupts to the computer. This essentially 
asks for the phase information for the appropriate sector to be 
output to the controlled phase oscillator. The subroutine will 
control the data flow so that the proper memory locations are 
output in sequence until all of the tone phases have been ad- 
justed. Control is then returned to the Executive Program to 
await the next SPM interrupt. 

In total, then, six interrupts will be required to perform 
the control required. 

3. REFERENCE TERMINAL EQUIPMENT 

Each Reference Terminal (R/T) is simply a transponder 
which returns the signals to the GCC, via the satellites, to 
provide a known range reference against which the system pertur- 
bations caused by satellite position uncertainty and ionosperic 
range bias can be compensated. 

A block diagram of the R/T is shown in Figure 4-8 which 
is the basis for the following discussion of the circuits and 
hardware involved. 

The link parameters derived in Appendix A assume 0-dB 
gain antennas at the A/C and R/T. This is not a constraint at 
the stationary R/T, however, and it is proposed that two rela- 
tively inexpensive 8 to 10-dB VHF helical antannas be used. 

This is done to improve the C/N at the GCC in that the gain of 
the antennas improves both the STS-to-R/T and the R/T-to-ATS 
links, while also minimizing satellite power sharing with 
unnecessary signals intended for the other satellite. 

The associated diplexers will be the same as those now 
used in the VHF Satellite Navigation Experiment. 
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Figure 4-8. Reference Terminal Block Diagram 













On the receive side the pre-amplifiers, first converters, 
and the phase-lock tracking converters will also be the same as 
existing equipment with slight modifications. In the transmit 
electronics the bandpass filters, balanced modulators, sideband 
filters, summers, and up-converters will be very similar to those 
in the OCC self-check subsystem. 

The control clock and gated pre-amplifiers perform the 
time multiplexing function as required for sector switching. 

The clock will be sufficiently stable so that synchronization 
may be accomplished by methods such as the following: (1) a 

telephone line from the GCC, (2) a WWV receiver at the R/T site, 
or (3) a synchronizing signal transmitted from the GCC. It is 
assumed that a WWV receiver will be available at each R/T site, 
so the second method is planned here. 

In review, the reference terminal with the exception of 
the antennas, clock, and the gated preamplifiers on the transmit 
side is made up of existing proven circuitry, and is therefore, 
a straight-forward design. 

4. AIRCRAFT TERMINAL EQUIPMENT 

4.1 Introduction 

In Section III-3.4 a brief description of the functions 
to be performed by the aircraft equipment was presented. In this 
section a more detailed description will be presented, and will 
include specific equipment recommendations. 

In specifying the aircraft equipment the primary criterion 
is to use as much equipment as possible that has been designed 
and is operationally proven. In addition, a particular airborne 
qualified computer which' may be made available by Texas 
Instruments on a loan basis during the experiment will be recom- 
mended. This approach will allow maximum concentration of effort 
on experimental goals rather than hardware design considerations. 
In the VHF Satellite Navigation Experiment (NAS 5-21079) and 
the Omega Position Location Experiment (NAS 5-10248) considerable 
equipment has been designed, built, and tested for some functions 
which are identical to those required in this experiment. This 
equipment will be recommended with little or no modification. 

The aircraft equipment can be divided into four main 
subsystems: (1) receiver demodulator, (2) phase measurement, 

(3) data processor, and (4) timing and control. The signal is 
received by the receiver demodulator subsystem where it is de- 
modulated by phase locking to the A/R tone in the ranging tone 
package. Since the phases of the incoming signal tones contain 
the desired information, one of the primary constraints on the 
receiver demodulator subsystem is that it preserve this phase 
information for presentation to the phase measurement subsystem 



» 




♦ 




for subsequent measurement by comparison with the on-board 
reference tones. Since the phase measurement subsystem must 
accommodate two sets of three non-harmonically related tones, a 
sampled data system implementation was selected. The data 
processor subsystem will be used to accumulate the phase data 
and calculate the latitude and longitude, while control of the 
aircraft system will be performed by the timing and control 
subsystem. Although a high stability crystal oscillator will be 
used, compensation must be made for any accumulated clock phase 
error resulting from a frequency offset or drift. By removing 
the accumulated phase error from the measured phase angle, (see 
Section II-3.1), the system will perform as if the on-board clock 
were closely synchronized with the ground clock. 

4.2 Receiver Demodulator Subsystem 

4.2.1 General Description 

The aircraft will simultaneously receive signals from 
both the ATS- I and ATS-III satellites. As discussed in Section 
II-2.3, the error caused by satellite movement restricts the 
maximum time the correct signals will be available for measure- 
ment (sector time interval). With this limit, both phase 
measurements must utilize as much of the interval as possible 
for signal integration and processing. Some of the time will 
be required for instrument settling; however, the A/R tone will 
be present continuously, which excludes PLL acquisition time 
from consideration. 

The block diagram in Figure 4-9 shows the receiver 
demodulator subsystem. This subsystem receives the signals 
from ATS- I and ATS-III. down-converts them to baseband, and 
filters out the six ranging tones. The ranging tones are then 
sampled by the phase measurement subsystem to be discussed in 
the next section. 

Since the two sections of the subsystems are identical 
for ATS-I and ATS-III except for the different channel frequencies 
of the tone package, it will be sufficient to discuss the re- 
ceiver demodulator for one satellite. As shown in Figure 4-9, 
the incoming signal is first down-converted to a 70-MHz IF, 
then further down-converted to 10.7 MHz using a phase-lock loop 
converter. The baseband output of the phase detector is then 
filtered using three 50-Hz bandwidth filters to separate the 
three tones. The tones are processed next by the phase measure- 
ment subsystem. 

4.2.2 Equipment Considerations 

All of the circuits enclosed by dotted lines have been 
designed and packaged for airborne equipment in the VHF satel- 
lite Navigation Experiment now being performed by Texas Instruments. 
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n The remaining circuits are currently being used in the ground 

control center of the above mentioned experiment and will be 
repackaged for airborne use with only slight modification. 

, 4.2.3 Self-Check Section 

r 

It is essential for the aircraft equipment to have the 
capability of completely checking the operation of the on-board 
system without requiring signals from the GCC. To perform this 
function, a self-check section of the Receiver Demodulator 
subsystem is included. The block diagram in Figure 4-10 shows 
the elements of the self -check section. In operation, the 
reference tones simulate ranging tones by modulating a self- 
generated A/R tone. These tones are then up-converted to the 
receive frequency channels of approximately 135.6 MHz. When 
the self-check switch is "on" this signal is switched into the 
two receivers. The resulting output on the display will be a 
fixed quantity which can be used not only to indicate system 
operation, but to provide some quantitative measure of its 
performance. 

4.3 Phase Measurement Subsystem 

4.3.1 General Description 

The primary function of the on-board equipment is to 
measure the phases of the ranging tones, as received from both 
satellites, with respect to an on-board generated reference. 

The receiver demodulator subsystem discussed in the previous 
section separates the six tones and presents them to the phase 
measurement subsystem for phase comparison with the reference 
tones. The equipment used to make the phase measurement is 
discussed in this section. 

The tone phase will be measured by a sampled data digital 
system. Each tone frequency (426 Hz, 707 Hz and 941 Hz) will 
be sampled at a different rate. The sampling rate is derived 
from the stable master clock. It can be shown that by sampling 
the received tone with one sampling signal that is in phase 
with the reference tone and another which is 90° out of phase 
with the reference tone, the phase of the received signal with 
respect to the reference can be found. In addition, the Nyquist 
sampling theorem states that a signal of bandwidth B need only 
be sampled at the rate of 2B. In this case, in which the tones 
are separated by a 50-Hz bandwidth filter prior to sampling, a 
sampling rate of only 100 Hz is required; however, the 50-Hz 
filters are not perfect and a sampling rate of 200 Hz was chosen 
to prevent aliasing. 

4.3.2 Sampler Timing 

Since the three tone frequencies are non-harmonic a 
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common sampling rate could not be chosen. This situation 
* necessitated the selection of a separate sampling system for 

each tone. The same frequency tone from both satellites may be 
sampled by the same sampler; however, since the tones do not 
have to be sampled on each cycle the sampler can be switched 
)■-' between them. A timing diagram of the sampling of the two 

941-Hz tones is shown in Figure 4-11. The sampling rate is 
derived from a clock which is four times 941 Hz so that both 
the in-phase and the quadrature phase samples can be decoded. 

Note also that the quadrature samples are taken 270° out of 
phase rather than 90°. This is to allow more settling time 
for the sample and hold and analog-to-digital converter. The 
quadrature sample will be inverted to obtain the 90° sample from 
the 270° sample, and its sign changed either in the computer or 
in the ADC. 

4.3.3 Subsystem Detailed Description 

A block diagram of the phase measurement subsystem for 
the 941-Hz tone is shown in Figure 4-12. The equipment for the 
other two tones is identical except for sampling rates. The 
counters used to generate the multiplexing and sampling signals 
are shown in this diagram for clarity; however, they will be 
located with the master clock so that they can be preset to the 
GCC timing. This will be discussed in Section IV-4.6. The 
sampling and digitizing are performed in the phase measurement 
subsystem while the integrating is performed in the data 
processor. 

The number of bits required by the accumulator and 
digitizer is discussed in paragraph 11-3.1. The number of bits 
in the accumulator remains the same for this system (10 bits) ; 
however, since the system proposed for this experiment will be 
using the same ranging tones as currently instrumented in the 
GCC, there will be a different number of bits required by the 
ADC. The number of bits in the ADC is dependent upon the 
signal-to-noise power density ratio expected at the input to the 
aircraft receiver. From the link calculations in Appendix A 
the worst-case E/N = 24.5 dB-Hz. In specifying the ADC require- 
ments the best E/N° should be chosen, and the best E/N that can 
be expected is less than 40 dB-Hz. 

» Since the tones will pass through a 50-Hz bandwidth filter, 

the signal-to-noise ratio at the ADC is 23 dB. This corresponds 
to a o = 1/20.1 radians. Since the ADC half scale must equal 
+ ir/4 radians, A = ir/2 n . From Section II-3.1 it is seen that 

* A = n/2 n < 1/20.1 

or 2 n > 63.2 

or n = 6 
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Figure 4-11. Sampling Signals for the 941-Hz Tones 
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After the sample has been digitized a signal must be 
sent to the computer to indi cate that digitization is complete. 
This is the DATA 1 AVAILABLE signal shown in Figure 4-12. Also, 
an identity code generated in the sampling decoder is read in 
with the data to identify the data word; i.e., satellite 1, 

941 Hz, in-phase, etc. 

4.3.4 Data Handling 

One problem that occurs because the three tone frequencies 
are asynchronous is that more than one ADC can have completed 
its cycle and be waiting to be read by the data processor at the 
same time. To solve this problem the first word that the pro- 
cessor reads in will tell the processor which ADC to read. As 
can be seen in Figure 4-13, bits BO, Bl , and B2 indicate which 

ADC is complete. The processor then reads each completed ADC 

by sending out a read command with the address of the ADC to be 
read. The address decoder receives the read command, gates in 
the proper A/C data and resets that DATA AVAILABLE flip flop. 
After the data is received by the processor it checks bits BO 
to B2 to see if another ADC is to be read. If true, the aboye 
process is repeated. If not true, a clear data command is sent 
and the processor waits for the next signal on the sense line. 

4.4 Data Processor Subsystem 

4.4.1 Required Tasks 

The Data Processor subsystem is required to perform 
several tasks which are summarized under general headings below. 

(1) Read and accumulate phase data 

(2) Read altimeter data and time of day 

(3) Calculate satellite constants from ephemeris 
function 

(4) Calculate latitude and longitude 

(5) Output data to the display 

(6) Determine if system is in self-check mode 

This section will describe what equipment is necessary to perform 
the tasks listed above and then describe and justify the choice 
proposed for this experiment. 

4.4.2 Hardware Requirements 

Some of the considerations in selecting a data processor 
are word length, speed of operation, input/output, and data/ 
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Figure 4-13 (a). Data Input Gates and Word Decoder 
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program storage. Examining these one at a time will serve to 
determine the equipment requirements. The choice of word length 
can be based upon processor considerations such as number of 
operation codes and directly addressable words in memory, or on 
data word size. Since the processor does not need to have 
exceptional programming flexibility the word size will be chosen 
based on the data. Each word to be transferred in is 13 bits 
long (6 bits data, 4 bits identity code, 3 bits ADC identity). 
This would indicate that a 16-bit processor would be appropriate 
allowing some margin for increased data resolution if necessary. 

A 16-bit word also allows a floating point number to be repre- 
sented with 22-bit accuracy with two words. This accuracy has 
been proven sufficient for making latitude and longitude cal- 
culations with programs now in operation at Texas Instruments. 

One method to estimate the amount of time required to 
perform a particular program is to count the number of mathe- 
matical operations and then multiply by the number of computer 
operations required for each. This will give a rough total 
number of computer operations. The computer cycle time may be 
deduced knowing the amount of time available for the calculation. 
In this case, assuming all calculations will be performed in 
floating point arithmetic, there will be approximately 350,000 
operations. The calculation will be performed when the ground 
station is transmitting phase data for the sectors other than the 
one in which the aircraft is located. This means that the time 
available for calculation is on the order of seconds. If one 
second is chosen for calculation time, a two-microsecond cycle 
time processor will be sufficient. The calculation will be 
complete in 700 milliseconds allowing 300 milliseconds to output 
the data. 

The input/output requirements for this data processor 
are not very demanding. In the discussion of the phase measure- 
ment subsystem it was mentioned that the sampling rate will be 
200 Hz. The 941-Hz tones will be sampled in-phase and 270° 
out of phase in the same cycle; therefore, there will be a 
minimum time of 785 microseconds between 941-Hz samples. Also, 
it is possible to have data samples ready from the two other 
samplers during this time. The input rate must, therefore, be 
greater than one 16-bit word per 261 microseconds; however, 
since the computer must perform several operations between data 
inputs, the time for inputting will be reduced. For 30 computer 
operations between data word inputs the time allowed to transfer 
one word into the computer is 200 microseconds. For data output 
there are three items: latitude, longitude, and time of day. 

If latitude and longitude are outputted in degrees, minutes, and 
seconds, and time of day in hours, minutes and seconds, the 
required output is 6 words after each calculation. 



Another consideration is the interrupt and sense line 
requirement. Since there is more than ample time to make the 
position calculation the data processor may be completely 
dedicated to data input and then completely dedicated to cal- 
culation of position. In other words, an interrupt capability 
is not required. A number of sense lines are required, however, 
so that the data processor will know when to read data and when 
to start calculating. 

It is obvious that such features as direct memory access 
and channel controllers are not necessary. In fact, while 
airborne all that is required for input and output is a 16-bit 
data bus and buffer registers. For programming the data pro- 
cessor on the ground it is suggested that some type of I/O 
device such as a teletype be used. This will be especially 
useful for an experimental system in which changes in the pro- 
gramming will be required more frequently than in an operational 
system. 


The final point to be discussed is the data and program 
storage requirements. The data storage for this system is 20 
words for each calculation. If it is desired to store the in- 
termediate values of certain variables used during the calcula- 
tion there are 120 words per calculation. For an experimental 
system it is desirable to be able to store the intermediate 
values for many calculations for playback and analysis on the 
ground after each flight. Obviously some compromise must be 
made concerning what will be stored. The suggested approach in 
the experimental system is to use as much storage area as is 
available for storing data which may be used in post real-time 
analysis. To accurately estimate the amount of storage for the 
program would require the program to be written. Although the 
exact program for this system has not been written a similar 
program was written to make the same position calculation at 
the ground control center in the VHF Satellite Navigation 
Experiment. This program took less than 4K of memory. In 
conclusion, it is recommended that an 8K memory for this experi- 
ment be available for program changes and storage of data for 
ground-based analysis. 

To summarize the data processor requirements Table 4-1 
is included below. 

Table 4-1. Data Processor Requirements 


Word Size 

Cycle Time 

Input/Output 

Interrupts 

Sense Lines 

Storag* 

16 bits 

2 u sec 

16-bit data 
bus + buffer 
registers 
(airborne) 
teletype 
(ground) 

None 

4 

8K 
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4.4.3 Equipment Considerations 

There are several alternatives in choosing a data pro- 
cessor to meet the requirements presented in the previous 
section. First, a data processor could be designed specifically 
for this experiement. Second, a non-airborne qualified computer 
can be modified to enable it to withstand the environment ex- 
pected in this experiment. Third, a qualified airborne com- 
puter can be chosen. The third choice has many advantages in 
that its performance is guaranteed for an airborne environment 
and design of a new data processor will not be required. A 
major problem with airborne computers is cost. Rigid vibration, 
thermal, RFI, and altitude specifications necessitate high cost 
and make the purchase of an airborne computer unreasonable for 
this experiment. However, there is a possibility that for the 
duration of this experiment Texas Instruments will be able to 
furnish the CP-967 (V)/UYK computer on a loan basis. 

4.4.4 CP-967 (V) /UYK Computer 

This computer is in operation for the Naval Ship Systems 
Command and its performance has been demonstrated. Although 
it is not packaged for airborne use it is small enough to fit 
on board many aircraft. In addition, software has been devel- 
oped for this computer, such as assembler, teletype I/O rou- 
tines, and fixed and floating point arithmetic. Some of the 
characteristics of this computer are summarized in the table 
below. 


•f 


f 


Table 4-2. CP-967 (V)/UYK Digital Data 

Computer 

Volume 

9 . 6 cubic feet 


Weight 

215 pounds 


Power Requirement 

115 +20 VAC, 400 Hz 

+5%, sin' 

Input Power 

400 watts average 


Size 

3 2 . 8 1 " H x 21 . 25"W x 

24 . 00"D 

Environment 

MIL-E-16400F Amendment 2 

Shock 

MIL-S-9 01 Grade A 


Vibration 

MIL-STD-167, Type 1 


Interference 

MIL-I-16 910 


Word Size 

16 bits 


Cycle Time 

2.0 microseconds 


Storage 

8 K to 64K 


Input/Output 

ASR 32 



phase 
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4.4.5 Flow Charts 


The functions of the data processor can best be expressed 
by a flow chart. For clarity a very general flow chart will be 
shown and then each block will be broken down into a more detailed 
flow chart. The flow charts are shown in Figures 4-14 through ^ 

4-18. 

4.4.6 Display 

As previously mentioned, the outputs to be displayed are 
latitude and longitude in degrees, minutes, and seconds, and 
time of day in hours, minutes, and seconds. This will require 
a storage of six 16-bit words. The display can be implemented 
with six 16-bit shift registers and 20 BCD decoder lamp drivers 
or with the serial type display now being tested at GSFC in 
conjunction with the VHF Satellite Navigation Experiment. 

4.5 Aircraft Antenna 

The aircraft antenna recommended for use in this experi- 
ment is the Dome and Margolin DMC33-2 VHF Airborne Satellite 
Communication Antenna. Extensive tests have been performed 
with this antenna in conjunction with investigative work now 
being performed by Texas Instruments on VHF sidetone ranging 
using the ATS-I and ATS-III satellites. The experience gained 
from using this antenna will be very useful in this experiment. 

Some of the characteristics are tabulated in Table 4-3. 

Table 4-3. DMC33-2 Satellite Communications 

Antenna Specifications 


Frequency Range 
Operating Frequencies 
Bandwidth at each Frequency 
VSWR 

Radiation Patterns 
Az imuth 
Polarization 
Power Handling (CW) 

Gain (Referred to a CP 
Low Angle Average 
Minimum 

High Angle Average 
Minimum 


135 to 150 MHz 

135.6 and 149.2 MHz 
1 MHz 

< 2:1, typically 1.5:1 

Omnidirectional 
LH circular or RH circular 
500 W 
) 

+1 dB 
Isotropic 
+ 3 dB 
Isotropic 


«• 
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Figure 4-14. General Flow Chart 
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(TWO SATELLITES) 


Figure 4-16. Calculate Satellite Constants Flow Chart 
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Figure 4-16 (cont.). Calculate Satellite Constants 

Flow Chart 


4-32 










m = SATELLITE (1 to 2) 
n = TONE (1 to 3) 


(TWO SATELLITES) 


Figure 4-17. Calculate Range Flow Chart 
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Figure 4-17 (cont'd.) Calculate Range Flow 

Chart 
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Figure 4-18. Calculate Aircraft Latitude and Longitude Flow 

Chart 











Figure 4-18 (cont'd.) Calculate Aircraft Latitude 
and Longitude Flow Chart 
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Table 4-3. DMC33-2 Satellite Communication 

Antenna Specifications (Continued) 

Impedance 50 ohms 

Weight 50 lbs 

4.6 Timing and Control Subsystem 

4.6.1 General Description 

The Timing and Control Subsystem coordinates the actions 
of the Phase Measurement Subsystem and the Data Processor Sub- 
system. The heart of this subsystem is the master clock which 
was selected based on the requirements discussed in Section 
II-3.1. The master clock controls the sector decoder, all of 
the sampler decoders in the phase measurement subsystem, and 
the elapsed time from synchronization counter. In addition, 
the time of day clock, sector select switch, and self -check 
switch are located in this subsystem. A block diagram is shown 
in Figure 4-19. 

4.6.2 Clock Updating 

In this experiment the master clock will not be updated 
while in flight. All updating will be performed on the ground 
at the GCC since the phase error prediction technique will 
provide the required ranging accuracy over a 24-hour period. 

In addition to the updating of the master clock, the counters 
which control the samplers in the phase measurement subsystem 
must be updated on the ground at the GCC. This is required 
because those counters which count down the ranging tone refer- 
ence frequencies from the master clock must be synchronized with 
the corresponding counters at the GCC. If this were not so, 
there would be a fixed delay between the on-board reference and 
the reference at the GCC. Since the entire system is based on 
these two being identical the delay would cause an incorrect 
phase measurement. In addition, the sector decoder must also be 
updated at the GCC to insure that the aircraft will be able to 
recognize the proper sector for making phase measurements. 

The fact that the master clock and the two decoders 
must be synchronized at the GCC poses some operational difficul- 
ties. First of all, the clock and counters must be portable so 
that they may be carried from GCC to aircraft while in operation 
Secondly, every precaution must be taken so that the synchroniza 
tion is not disturbed while in flight. In order to accomplish 
this a battery and power supply must be part of the timing and 
control subsystem. Further, RFI shielding must be provided to 
protect and clock and counter synchronization. 




4-38 


Figure 4-19. Timing and Control Subsystem 




4.6.3 Sector Selection 


As discussed in previous sections, there will be three 
different sectors in this experiment. Each time the aircraft 
moves into a different sector it must perform the phase measure- 
ments in a different time slot. This will be under the control 
of the sector decoder and the sector select switch. As shown 
in Figure 4-20, the output of the sector decoder is three signals 
SI, S2 and S3. SI is high when the GCC is transmitting signals 
for sector one, S2 is high when the GCC is transmitting signals 
for sector two, etc. When in sector one, the aircraft switch 
must be in position one (sel 1) in order to allow the samplers 
to sample when the GCC is transmitting sector one signals. 

Note that there is some dead time between the sector signals. 

This is to prevent the on-board system from sampling during the 
step change in phase from the GCC when changing sectors. 

4.6.4 Time of Day Clock 

A clock which reads the time of day must be available on 
board because the calculation of satellite constants from the 
ephemeris data requires the time of day to determine satellite 
position. The clock will have a BCD output of hours, minutes, 
and seconds which is fed to the data input gates of the Data 
Processing Subsystem. 

4.6.5 Master Clock 

The master clock chosen is the Motorola S1081A Oscillator. 

. -9 

This clock has the required long term stability of 1 x 10 per 

-9 . 

day and 5 x 10 over its temperature range. In addition, it 

-9 , , 

may be adjusted to within + 1 x 10 , and stabilizes to within 

-9 . 

1 x 10 in only 30 minutes. 

5. EXPECTED SYSTEM PERFORMANCE 

As stated earlier the experimental system is not the 
optimum configuration. It is generally true that a specifically 
tailored design will perform better than a rebuilt adapted old 
system. It is also true, however, that operation and concepts 
can be verified without attaining the ultimate in performance. 
Consequently the evaluation of the expected system performance , 
in the areas in which it may fall short in technological capa- 
bility, may require more data points for statistical inter- 
polation or carefully designed experimental procedures to 
circumvent the problem. 

Table 4-4 is a tabulation of errors caused by various 
portions of the hardware, as obtained from the study of equipment 





capability and expected design performance. It is emphasized 
that these numbers are estimates. These values are based upon 
the study results, similar design experience, and in some cases 
computer simulation. 

I 

Table 4-4. Range Errors Attributable to Hardware 
Error Source RMS Range Error (meters) 


GCC Processing and Control Errors 

Phase Measurement, 0.8° rms 
(2-seconds of integration 
time at 40 dB-Hz) 

Quantization 

Time Delay Between Sectors 

A/C Processing Errors 

Phase Measurement, 0.4° rms 
(2-seconds of integration 
at 40 dB-Hz) 

Quantization 

Altitude Uncertainty 

(Equivalent for 200-meter 
altitude error) 

Reference Terminal Delay 
Uncertainty 

Uncompensated Satellite 
Transponder delay 


In addition to the errors attributable to hardware 
performance, other ranging errors are encountered as a result 
of the general system component configuration. These errors 
are listed in Tables 4-5 and 4-6. 


Table 4-5. Multipath Ranging Errors 


Satellite 

Elevation 

Type of 
Fading 

Nominal Ranging 
Error (km) 

426Hz 707Hz 941Hz 

Worst Case 
Ranging Error 
426Hz 707Hz 941Hz 

0°-30° 
30°-50° 
50° -70° 

Specular 

Diffuse 

Diffuse 

0.76 0.73 0.72 

13.8 8.1 6.2 

9.5 5.7 5.3 

3.8 3.6 3.7 


710 

178 

75 

355 

30 

100 

120 

30 



Table 4-6. Ranging Errors Attributable to System Configuration 


Error Source 


RMS Range Error (meters) 


Ionospheric & Tropospheric 
Perturbations 

120 

Uncompensated Satellite 
Position Uncertainty 
(Worst case: Sector 

Boundary) 

500 

Satellite-to-A/C Multipath 


(0°-30° nominal) 

720 

/ 

(30°-50° nominal) 

6200 


The multipath errors for the 0° to 30° satellite 
elevation angles apply to the experimental system. In order to 
provide a larger area of coverage in an operational system 
satellite elevation angles greater than 30° will be encountered. 
The use of a higher frequency ranging tone (10.2 kHz) will tend 
to cancel the increased multipath effects. 

This combination of all errors listed in the above tables 
gives the total root mean square error in range that may be 
expected for the experiment. Positioning the ATS-I and ATS-III 
satellites at 50°W and 150°W longitude provides a sector of 
satellite coverage centered at the mid point of the continental 
United States. The minimum satellite elevation angle in this 
sector is 10° and the maximum is approximately 30°. Therefore, 
the lower value of multipath error listed in Table 4-6 is 
applicable to the experimental system. The sum of all root 
mean square errors is formed by taking the square root of the 
sum of the squares and is equal to 1210 meters of range error. 

At the latitudes of interest (approximately 35°N) this range 
error translates to 1600 meters of radial position location 
error . 
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SECTION V 


EXPERIMENT OUTLINE 

This study has shown that in principle it is feasible to 
perform on-board navigation using the ATS-I and ATS-III as re- 
lays for the sidetone ranging signal generated at the modified 
OPLE Control Center. The primary objective of an experiment 
following the study would be to demonstrate this feasibility; 
i.e. , to demonstrate that the navigation function can be 
reliably performed in a predictable manner. 

In Section III of this report a system design was 
developed around three major unique features: 

(1) Real-time differential sidetone ranging in a 
Range/Range mode 

(2) Real-time compensation for satellite position 
uncertainty and ionospheric range bias 

(3) Sector time-slot organization for system control 

Presented in the following paragraphs is a general outline of 
an experiment which will illustrate the feasibility and 
practicability of a system with these features, thereby serving 
as a testbed for more advanced satellite ranging navigation 
systems not constrained by the specific characteristics of the 
OPLE Control Center, ATS-I and ATS-III. 

1. SYSTEM ENGINEERING TESTS 

As a first step certain static system parameters must be 
verified. Principal among these is the verification of the 
expected C/N Q for all communication links. In the course of the 
VHF Satellite Navigation Experiment it was found that the ATS 
transponder noise characteristics differed from the predicted, 
which required modification of the design signal power sharing 
within the satellite transponders. For the Dual Satellite 
Navigation Experiment the option of modifying the ground control 
center (GCC) and the reference terminal (R/T) transmitted power 
levels must be maintained to optimize the satellite transponder 
power sharing and maximize the probability of usable signal 
acquisition . 

Because the ranging system is basically a timing system, 
the calibration of time delays through the system components 
must also be verified at the beginning of the experiment. These 
include the delay through the R/T transponder and the delay 
through the receiving system at the GCC for use in subsequent 
R/T range calibration, the aircraft receiving system delay to 



calculate the fixed range bias for use in the on-board computer 
subsystem, and the GCC transmitter subsystem delay for use in 
predicting the phase of the compensated transmitted tone. 

The principal feature of the aircraft equipment per- 
formance to be verified in the engineering tests is the 
stability of the on-board clock. When the drift characteristics 
of the clock are predictable, the phase and frequency may be 
"over compensated" when initialized to drift through the usable 
region during the period of testing. It is anticipated, however, 
that the on-board computer will be programmed to compensate for 
a well-behaved clock error even though this technique requires 
very high precision in the measurement of initial frequency 
offset. In either case, to preclude the requirement for a 
costly clock the stability characteristics of the clock must 
be verified at the beginning and periodically throughout the 
experiment. 

2. REFERENCE TERMINAL RANGE CALIBRATION 

In the second group of tests the capability to perform 
real-time compensation for satellite position uncertainty and 
ionospheric range bias is verified. In addition to ascertaining 
proper sequencing among the R/T sectors and GCC signal acquisi- 
tion from each R/T, it is in these tests that the GCC algorithm 
for range computation and compensation with the satellite 
ephemeris update program will be validated. Since the geo- 
graphical positions of the R/T's will be known, the applicable 
equipment time delays will have been measured, the satellites' 
true positions will be available from published ephemeris data, 
and an estimate of the ionospheric range bias may be made, the 
expected GCC- to- R/T range can be calculated and compared with 
the measured results. These measurements and corresponding 
transmitted phase adjustments may then be monitored through the 
diurnal ionospheric changes and, with cognizance of the true 
satellite ephemeris, the proper operation of the real-time 
constant phase compensation may be verified. 

3. STATIC AIRCRAFT LOCATION TESTS 

These tests represent the first evaluation of the navi- 
gation capabilities of the system and immediately precede the 
navigation demonstration. It is here that the position location 
algorithm in the aircraft terminal is verified with the terminal 
placed at a known location. A two-step plan will be followed: 

Measurement and computation of range from each 
satellite, verifying the operation of the range 
measurement subsystem and validating the constant 
range compensation performed by the GCC 


( 1 ) 




(2) 

Combination of both satellite ranges and altitude 
to verify the position location computation 

c 

During these tests the predictability of the effects of various 
error sources will also be determined. These sources include; 


(1) 

Aircraft altitude uncertainty determined by simu- 
lating changes in altitude 


(2) 

On-board clock timing error, determined by off- 
setting clock phase 


(3) 

Received S /N influence, determined by varying the 
GCC-transmitted signal power and/or the aircraft 
antenna/receiver gain 


(4) 

Satellite position uncertainty, determined by 
modifying the precision of the satellite ephemeris 
predict in the aircraft and GCC algorithms 


(5) 

Time delay between sectors, determined by modifying 
the sector timing intervals 


(6) 

Data processor error, determined by reprogramming 
the on-board computer for various digitizing levels 
or data word lengths 


4. NAVIGATION DEMONSTRATION 

In the last series of tests the aircraft terminal will be 
mounted on-board the aircraft for flight demonstration of the 
system to culminate the essential tests. Position determination 
accuracy will be established by comparison of the satellite 
navigation measurement with a prearranged external navigation 
reference . 

In addition to the capabilities demonstration, two 
additional sources of position error will be investigated during 
in-flight tests; 

(1) Multipath 

The gross effect will be measured by performing tests 
in the same hour with the aircraft grounded and with 
the aircraft airborne. 

(2) Geometric Dilution of Precision (GDOP) 

The effect of reference terminal-to-aircraf t separa- 
tion distance will be examined, as well as the effect 
of aircraft transition into an adjacent sector. In 
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addition, both the aircraft and a mobile reference . 

terminal will be moved a considerable distance to 
verify the expected gross GDOP effect. 

5. EXTENDED EXPERIMENTATION 

With a successful experiment incorporating the essential 
tests outlined above, the Dual Satellite Navigation System as 
developed in this study can form the framework for considerable 
additional investigations. A list of potential extensions 
would include : 

(1) For further testing the GCC may be again modified 

to incorporate the optimum higher frequency (10.2 kHz) 
fine ranging tone discussed in Section II of this 
report. 

(2) With the higher measurement precision permitted by 
the 10.2-kHz ranging tone more extensive multipath 
investigations may be conducted, including examina- 
tion of the effect of multipath spreading at various 
climb and descent rates, and the influence of terrain 
(land/water) at various aircraft velocities. 

(3) The system may be combined with the VHF Satellite 
Navigation System in an integrated VHF navigation/ 
surveillance system. In essence, this may be 
accomplished merely with the installation of a VHF 
transmitter in each aircraft. 

(4) Additional navigation modes may be investigated; 
e.g.. Range Rate, Range Difference. 

(5) In-flight on-board clock up-dating techniques may be 
instrumented for feasibility demonstration and 
examination . 

Although the above represents only a partial listing of 
candidate programs, the list is nevertheless formidable. With 
the hardware in this system as a nucleus , near-term experimen- 
tation in the field of VHF satellite navigation, surveillance, 
and data relay may be cost-effectively conducted. 
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SECTION VI 


PROGRAM COST ESTIMATION 

To maintain the Dual Satellite Navigation System and 
Experiment costs as low as possible four major assumptions have 
been made : 

(1) Every effort has been made to minimize the 
required modification of the GCC. 

(2) A majority of the circuits and the construction 
techniques used in the development of the VHF 
Satellite Navigation Experiment remote terminals 
are applicable to the R/T and A/C terminal 
requirements . 

(3) A suitable computer for airborne use will be made 
available. 

(4) A suitable navigation data display for airborne 
use will be made available as required on a loan 
basis from Goddard Space Flight Center. 

With these considerations it is estimated that the development, 
production, testing, and integration of the system recommended 
in this study will require between 12 and 18 months to complete, 
and the essential tests in the experiment outline will not 
require more than four months to complete. It is further 
estimated that the total program can be completed in three 
phases as follows: 

Phase I 


Ground Control Center 


Engineering Labor 21 man-months 

Redesign of OPLE receivers 
Modification of SELF-CHECK Subsystem 
Design of PHASE SHIFT Subsystem 
Modification of SELF-TEST Subsystem 
Modification of Software 

Material 


Transmitter power amplifier $1,000 
Frequency synthesizer 7,000 
Antenna dip lexer 5,000 
Components and hardware 7,000 



Reference Terminals 


Engineering Labor 

Dual channel design 

Timing and control logic design 

Material 

Antennas (4) 

Components and hardware (3) 

TOTAL COST FOR PHASE I - Approximately 
Phase II 


Aircraft Terminal 


Engineering Labor 

Receiver and demodulator design 
PHASE MEASUREMENT Subsystem design 
Timing and Control logic design 
SELF-TEST Subsystem design 
Software design 

Portable clock integration design 

Material 

Antenna 
Master clock 
Components and hardware 

TOTAL COST FOR PHASE II - Approximately 

Phase III 


System Integration and Check-Out 
Engineering Labor 
Minimal Experiment 
Engineering Labor 

TOTAL COST FOR PHASE III - Approximately 
Total Program Effort 
Engineering Labor 
Material 

TOTAL PROGRAM COST - LESS THAN 


3 man-months 

$4,000 

11,000 

$110K 

30 man-months 


$3,000 
3,000 
20 ,000 

$120K 

16 man-months 

14 man-months 
$ 9 OK 

84 man-months 
$61,000 
$350,000 


The figures itemized above are presented with a confidence 
originating from a complete and thorough understanding of the 
GCC, which is based upon Texas Instruments original effort in 
the design, development, and manufacture of the GCC, the opera- 
tion and maintenance of the GCC throughout the OPLE Experiment, 
and its modification for use in the VHF Satellite Navigation 
Experiment. In addition, the effort involved in these two 
programs has provided us with an appreciation of the capabilities 
and the limitations of the ATS -I and ATS-III with respect to the 
use of the VHF transponders. We are confident, therefore, that 
the program outlined in this report, incorporating the concepts 
developed in the study, will prove the feasibility of on-board 
VHF sidetone ranging navigation. 
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APPENDIX A 


VHF LINK PARAMETERS 


Al. INTRODUCTION 

This appendix summarizes the calculations performed to 
characterize the VHF links associated with the ATS-I and ATS-III 
satellites. It covers the VHF link calculations, satellite 
power sharing, and sicrnal-to-noise (SNR) ratios, and the develop- 
ments presented are sufficiently general to be applicable to 
both the VHF Satellite Navigation Experiment configuration and 
the Dual Satellite Navigation configuration discussed in this 
report . 

A2 . LINK CALCULATION 

The link loss-to-noise power density ratios (L/N ) for the 
various VHF links are summarized in Table A-l . The parameters 
used to obtain these values for the 135.6-MHz links are pre- 
sented in Table A-2 , and those for the 149.22-MHz link are pre- 
sented in Table A-3. The aircraft antenna parameters were taken 
from [19], and the satellite parameters were obtained from [20], 
The multipath fade margins on the aircraft/satellite links are 
derived in Appendix C; the fade margins on the ground-to- 
satellite links are derived from the ionospheric scintillation 
data in Appendix B; the antenna temperatures are derived in 
Section A5 ; and the various noise-power densities (N ) are 
derived from the material presented in Section A6 . 


Table 

A-l. L/N 

for 

0 

the VHF Satellite Links 


Link 

Frequency 

Nominal 

L/M o 

Worst-Case 

L/N 

o 

Units 

GCC-to-ATS-I 

149.22 

MHz 

41 

25.9 

Hz/dBw 

ATS -I-to-A/C 

135.6 

MHz 

28.3 

14.5 

Hz/dBw 

A/C-to-ATS-I 

149.22 

MHz 

26.5 

12.2 

Hz/dBw 

ATS-I-to-GCC 

135.6 

MHz 

43.1 

27.2 

Hz/dBw 

GCC -to- ATS-III 

149.22 

MHz 

41.5 

26.2 

Hz /dBw 

ATS-III-to-A/C 

135.6 

MHz 

30 

16.2 

Hz/dBw 

A/C- to- ATS -III 

149.22 

MHz 

27 

12.5 

Hz/dBw 

ATS-III-to-GCC 

135.6 

MHz 

44.8 

28.9 

Hz/dBw 



Table A-2. Link Parameters for 135.6-MHz Dovm Links 


Link Parameter 

Nominal 

Worst-Case 

Units 

Free Space Loss 

-166 

-167.3 

dB 

Polarization Loss 

-3 

-3 

dB 

Satellite-A/C Fade Margin 

-2 

-7 

dB 

Satellite-A/C Net Propagation Loss 

-171 

-177 .3 

dB 

GCC-Satellite Fade Margin 

-1.3 

-9.4 

dB 

GCC-Satellite Net Propagation Loss 

-170.3 

-179.7 

dB 

ATS -I Transmit Antenna Gain 

9.5 

9.5 

dB 

ATS-I Transmit Losses 

-2.4 

-2.4 

dB 

A/C Antenna Gain 

0 

0 

dB 

A/C Antenna Ellipticity 

-3 

-4.3 

dB 

A/C Receiving Losses (estimate) 

-1 

-1.5 

dB 

ATS-I-to-A/C Net Link Loss (L 1A ) 

-167.9 

-175.4 

dB 

ATS-III Transmit Antenna Gain 

10.5 

10.5 

dB 

ATS-III Transmit Losses 

-1.7 

-1.7 

dB 

ATS-III-to-A/C Net Link Loss (L^) 

-166.2 

-173.7 

dB 

GCC Antenna Gain 

12 

12 

dB 

GCC Antenna Ellipticity 

-.7 

-.7 

dB 

GCC Receiving Losses (estimate) 

-1 

-1.5 

dB 

ATS -I -to-GCC Net Link Loss (L^ G ) 

-152.9 

-162.8 

dB 

ATS-III-to-GCC Net Link Loss (L^q) 

-151.2 

-161.1 

dB 

Aircraft Antenna Effective Noise 

1420 

7100 

°K 

Temperature 
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Aircraft Noise Figure (estimate) 

2.5 

3 

dB 
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Table A-2. Link Parameters for 135.6-MHz Down Links (cont) 


Link Parameter 

Nominal 

Worst-Case 

Units 

Aircraft Received Noise-Power 
Density (N Qa ) 

-196.2 

-189.9 

dBw/Hz 

l ia^ n oa 

28.3 

14.5 

Hz /dBw 

L 3A /N 0A 

30 

16.2 

Hz/dBw 

GCC Antenna Temperature 

1420 

7100 

°K 

GCC Noise Figure 

3 

3.5 

dB 

GCC Received Noise-Power 
Density (N QG ) 

-196 

-190 

dBw/Hz 

L 1G /N 0G 

43.1 

27.2 

Hz /dBw 

L 3G /N 0G 

44.8 

28.9 

Hz/dBw 

Table A-3. Link Parameters 

for 149.22- 

-MHz Up Links 


Link Parameter 

Nominal 

Worst-Case 

Units 

Free Space Loss 

-166.9 

-168.2 

dB 

Polarization Loss 

-3 

-3 

dB 

Satellite-A/C Fade Margin 

-2 

-7 

dB 

Satellite-A/C Net Link Loss 

-171.9 

-178.2 

dB 

GCC-S ate Hite Fade Margin 

-1.3 

-9.4 

dB 

GCC-Satellite Net Propagation Loss 

-171.2 

-180.6 

dB 

ATS -I Receive Antenna Gain 

7.5 

7.5 

dB 

ATS -I Receiving Losses 

-1.4 

-1.4 

dB 

A/C Antenna Gain 

0 

0 

dB 

A/C Antenna Ellipticity 

-2 

-4.3 

dB 

A/C Transmit Losses (estimate) 

-1.5 

-2 

dB 
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Table A-3. Link Parameters for 

149 . 22-MHz 

Up Links 

(cont) 

Link Parameter 

Nominal 

Worst-Case 

Units 

A/C-to-ATS-I Net Link Loss (L^) 

-169.3 

-178.4 

dB 

ATS-III Receive Antenna Gain 

7.5 

7.5 

dB 

ATS-III Receiving Losses 

-1.3 

-1.3 

dB 

A/C-to-ATS-III Net Link Loss (L^) 

-169.2 

-178.3 

dB 

GCC Antenna Gain 

12 

12 

dB 

GCC Transmit Loss (estimate) 

-1 

-1.5 

dB 

GCC Antenna Ellipticity 

-.7 

-.7 

dB 

GCC- to- ATS -I Net Link Loss (L^J 

-159.8 

-164.7 

dB 

GCC-to-ATS-III Net Link Loss (L G3 ) 

- -154.7 

-164.6 

dB 

Satellite Antenna Noise Temperature 

1120 

5520 

°K 

ATS-I Noise Figure 

4 

4 

dB 

ATS -I Noise-Power Density (N Q1 ) 

-195.8 

-190.8 

dBw/Hz 

l ai/ n oi 

26.5 

12.9 

Hz/dBw 

^l^Ol 

41 

25.9 

Hz /dBw 

ATS-III Noise Figure 

3.5 

3.5 

dB 

ATS-III Noise-Power Density (N Q3 ) 

-196.2 

-190.8 

dBw/Hz 

L A3/ N 03 

27 

12.5 

Hz/dBw 

L G3 //N 0 3 

41.5 

26.2 

Hz/dBw 


Carrier-to-noise power density ratios (C/N q ) for the 
various links are presented in Table A-4 for a GCC transmitter 
power of 50 watts (17 dBw) , ATS -I transmitter power of 40 watts 
(16 dBw) (20], ATS-III transmitter power of 50 watts (17 dBw) 
[20] , and an aircraft transponder power of five watts (7 dBw) . 
The weakest link is the aircraf t-to-satellite link, which has a 
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C/N q at least 10 dB lower than the other links for these power 
levels and will, therefore, determine the ultimate C/N q that 
may be achieved in the transponding mode . 


Table A-4. C/N for the VHF Satellite Links 

o 


Link 

Frequency 

Nominal 

C/N 
' o 

Worst-Case 

C/N 

o 

Units 

GCC-to-ATS-I 

149.22 

MHz 

58 

42.9 

dB-Hz 

ATS -I- to- A/C 

135.6 

MHz 

44.3 

30.5 

dB-Hz 

A/C-to-ATS-I 

149.22 

MHz 

33.5 

19.2 

dB-Hz 

ATS-I-to-GCC 

135.6 

MHz 

59.1 

43.2 

dB-Hz 

GCC-to-ATS-III 

149.22 

MHz 

58.5 

43.2 

dB-Hz 

ATS -II I -to- A/C 

135.6 

MHz 

47 

33.2 

dB-Hz 

A/C-to-ATS-I I I 

149.22 

MHz 

34 

19.5 

dB-Hz 

ATS-III-to-GCC 

135.6 

MHz 

61.8 

45.9 

dB-Hz 


A3. SIGNAL-TO-NOISE RATIOS AND SATELLITE POWER SHARING 

A3 . 1 Satellite Power Sharing 

From Section A6 , the satellite power allocated to the 
GCC signal (pilot tone plus ranging tones) is 


-y 


% 


P sg =V [lt| VV ll t S V t( V P g )(H osV 1 (A - 1 > 


where the quantities are defined in Figure A-l. Utilizing the 
nominal values from Table A-3 for the aircraf t-to-ATS-III and 
ATS-III-to-GCC links, together with the power levels specified 
in Section A2 and B = 100 kHz, the satellite power allocated 
to the GCC-transmitled signal is calculated as 

P sg = 40/ [1.142 + 3.55 x 10 _3 m] (A-2) 


where m is the number of transponders using the system simul- 
taneously. For m = 1, P =35 watts; and for m = 5, P g = 34.4 


The satellite power allocated to the remote platforms is 
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P st = P s /[1 + (P g /P t ) (L gs /L ts ) + (B s /P t ) (N os /L ts )] (A-3) 


/T 


- 40/ [ 1 + 322/m] (A-4) 


For m = 1, P gt = 124 mW; and for m = 5, P ^ 611 mW. 

The thermal noise power transmitted by the satellite is 

P sn = V* 1 + <V B s } (L gs /N os } + ( W < L t S / N os^ < A ‘ 5 > 


« 40/ [8.25 + .025 m] 


(A- 6) 


For m = 1, P - 4.83 watts; and for m = 5, P -4.77 watts, 
sn ' ' sn 

The satellite power allocated to the pilot tone is 

K .P , and the GCC pilot tone power is K .P . From Section A6 
sj s' c c gj g 


K ./K . = 1 + (P./P ) (L, /L ) + (B /P ) (N /L )] (A-7) 

gj' sj t' g ts' gs' ' s / g os' gs 7 


1.142 + 3.55 x 10" 3 m 


(A-8) 


For K . = 0.25 and m = 1, K . - 0.286, and for m = 5, K . - 0.29. 
sj go ' sj 

The above equations illustrate that approximately 71 percent of 

the GCC transmitter power, or 35.5 watts, will be available for 

the ranging tones. 

A3 . 2 Transponder Signal-to-Noise Ratios 

For a pilot tone requiring 25 percent of the satellite 
transmitter power, the pilot tone carrier-to-noise power density 
ratio (C/N q ) is found from Section A6 to be 


(C/N ) . = K .P /[N ,/L . + P /B ] 
o j gj sg 7 ot 7 st sn 7 s 


(A-9) 


= 9.55 x 10 /Hz or 39.8 dB-Hz 
The tone package SNR is 


(C/N K = (1 - K .) P /[N ,/L . + P /B ] (A-10) 

7 o tones gj sg 7 ot 7 st sn 7 s J 


= 23.8 x 10 /Hz or 43.8 dB-Hz. 
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A3.3 GCC Signal-to-Noise Ratios 

The C/N^ for the pilot tone received at the GCC is 

(C/N )„• = K -P / [N /L + P /B ] 

' o' go gj gs / og 7 sg sn / s J 


( A-ll) 


= 1.25 x 10 5 /Hz or 50.9 dB-Hz. 

If a transponder bandwidth of 5 kHz is assumed, then 
from Equation (A-9) the tone package C/N at the transponder is 

(C/N) t = (C/N Q ) tj /B t = 39.8 - 37 = 2.8 dB 

The portion of the transponder power that is actually 
signal power is 


P^ = P ^[P • /(P ■ + 0.525 P . )] 

ts t sxg' sig sxg 


(A-12) 


- 0.656 P or 3.28 watts 


The tone package (C/N q ) at the GCC is 

< C /V tones - ' 656 P st/ [ V/ L sg + WV < A - 13 > 

= 10 3 /Hz or 30 dB-Hz. 

A4 . ANTENNA TEMPERATURES 

A4.1 Satellite Antenna Temperature 

If the satellite antenna beamwidth, 9 , is greater than 
the angle subtended by the earth, 0 , as shown in Figure A-2 , 
then the satellite antenna thermal noise temperature will be a 
combination of both the earth and cosmic temperatures. To a 
first approximation, the satellite antenna noise temperature 
may be expressed by 

T - T ( 0 /0_) 2 + T (0 3 - Q 2 )/Ql 
a e e' s c s e's 

= kT + (l-k)T ( A-14) 

s c 
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where 


290 °K 


T = 
e 

k = ( 0 e /0 s )2 (A-15) 

and T , the nominal cosmic temperature, is given by [21] 
c 

T = 290 A 2 

c 

= 2.61 x 10 ? /f 2 MHz ( °K) (A-16) 

The maximum value of T c is five times the value given by 
Equation (A-16) [21] . 

An approximate expression for the beamwidth of the 
satellite antenna in terms of its gain relative to an isotropic 
antenna is [22] 

0 2 = 27000/antilog 10 (G^/10) (A-17) 

When this expression is combined with the value for 0 given 
in Figure A-2, Equation (A-15) becomes 

k = (17.32) 2 /27000/antilog 1Q (G^/10) 

= 0.0111 antilog 1Q (G^g/10) (A-18) 

For a satellite receiving antenna gain of 7.5 dB , k = 0.0624. 
Substitution of this value and f^ = 149.22 into Equation 

(A-14) gives a nominal satellite receiving antenna temperature 

of 1120°K and a maximum temperature of 5520 °K. 

A4.2 Aircraft Antenna Temperature 

The nominal aircraft antenna temperature is given by 
Equation (A-16) . For a receive frequency of 135.6 MHz, Equation 
(A-16) gives a nominal antenna temperature of 1420 °K; and the 
maximum value will be 7100°K. Although the aircraft antenna 
temperature can reach almost 70000°K during precipitation static 
discharge [23] , the 7100°K worst-case value will be used until 
a time distribution of the antenna temperature can be obtained. 


* 


\ 




* 


A-10 



*' 


A5 . 


NOISE POWER DENSITY CALCULATIONS 


?> 


y 


« 


A5 . 1 Receiver Input Noise Temperature 

The thermal noise generated within the receiver along 
with that caused by transmission line losses may be referred 
back to the transmission line input in the form of the effective 
receiver input temperature, T e ff This is expressed in terms 

of the transmission line and receiver parameters by [24] 


T eff = T L (1/L “ X) + T R ( F ~ 1>/L (A- 19) 

where L is the transmission line loss (L = P ^/P . ) , T is 

out / in ' L 

the transmission line temperature, F is the receiver noise 
factor, and T D is the receiver temperature. If it is assumed 
that both T and T are at the "standaVd" temperature of 290°K, 

K 

then Equation (A-19) reduces to 


T eff = 290 (F/L - 1) 


(A-20) 


This is plotted as a function of F dB - L dB in Figure A-3. 

The receiver noise temperature is the sum of Equation 
(A-19) and the antenna temperature. Thus, 


= T . + T 

ant eff 


= T ant + 290 (F/L - 1) 


(A-21) 


A5 . 2 Receiver Noise Power Density 

The receiver noise power is given by 

P = KT B watts ( A-22) 

rn r 

23 

where K is Boltzman's constant (1.38 x 10 J/°K) , T is the 

receiver noise temperature, and B is the receiver bandwidth.. 

The receiver noise power density is given by 

N o = P rn / B = KT r watts/Hz (Joules) (A-23) 


A-ll 



eff 



















4 which is plotted in Figure A-4. N q is plotted as a function of 

antenna temperature, receiver noise figure, and transmission 
line loss in Figure A-5 . 

A6. derivation of link equations 

A6 . 1 Introduction 

The link equations will be derived from the power flow 
diagram shown in Figure A-l. Note that this diagram applies to 
both the VHF Satellite Navigation Experiment and the Dual 
Satellite configuration links. The satellite power sharing 
equations will be derived in Section A6.2, and the transponder 
and GCC receiver carr ier-to-noise power density ratios (C/N ) 
will be derived in Section A6 . 3 . 

A6.2. Satellite Power Sharing 

Satellite Received Power 


The total power received by the satellite is the sum 
of the thermal noise power plus the power received from the GCC 
and the remote platform(s). This may be expressed by 


P 

sr 


P L 

g g s 


+ P t L ts 


+ N 


( A-2 4) 


where the quantities are defined in Figure A-l. Note that the 
remote platform (transponder) power may be the power trans- 
mitted by several platforms, and that the "link loss" term is 
given by 


L = a G t G r ( A-25) 

where a is the free space attenuation; G and G are res- 

-L K 

pectively, the transmitting and receiving antenna gains; and 

i i 

L m and L are, respectively, the transmitting and receiving 
T R 

/ losses. 


Basic Satellite Transmitter Power Allocation 


The satellite power will be allocated to the various 
receiver inputs according to the fraction of the total received 
power that they represent. Thus, the satellite power allocated 
to the GCC signal is 
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Receiver Input Noise Temperature (°K) 
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Receiving Losses Plus Noise Figure (dB) 

Figure A-5. Receiver Noise Power Density vs. Receiving Losses and Noise Figure 














P sg - (P g L g S / P sr> P s 

= V [1 + < V P g> <L ts/ L g S > + < V P g) < N os/ L gs)> < A - 26 > 
the power allocated to the remoted platform(s) is 

P st - (P t L ts/ P sr» P s 

“ P s /[1 + (L gs/ L ts> + ' <V P t> (N os /L ts ) 1 (A - 27) 

and the transmitted noise power (excluding intermodulation 
noise) is 

P sn - p s 

■ V [1 + (P g /B s> (L gs/ N os> + (P t^s» IL t S / N os> 1 < A - 28 > 

Pilot Tone Power Allocation 

A tone is transmitted by the GCC to aid the system 
receivers in removing frequency jitter introduced by the ATS-III 
satellite. For a given pilot tone power transmitted by the 
satellite 

P sj = K sj P s ' 0 1 K sj 1 1 < A - 29 > 

the GCC will be required to transmit a pilot tone power of 


P . = K . P . 

9n gi g ' 


o < k . < 1 

- g: - 


(A -30) 


Substituting (A-29) and (A-30) into Equation (A-26) and solving 
for K . gives 
gi 

K gj “ K sj [1 + < P t /P g> <L ts/V + 'VV'^V 1 <A - 31 > 


A6 . 3 Signal-to-Noise Ratio Derivations 

Transponder Signal-to-Noise Ratio 
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The power received by the transponder is 


P tr = ‘ P sg + <VV P sn ! L st + N t < A ' 32 > 

The carrier-to-noise ratio (C/N) is 

( C /N) t - P sg /[B t (N ot /L at ) + B t /B s ) P sn ] (A-3 3) 

Multiplying both sides of Equation (A-33) by B gives the 
carrier-to-noise power density ratio 


< C ' N o>t - W tN ot' L st + P sn'V 

“ N ot /L st >:> P sn' B s 


(A— 34 ) 
(A-35) 


GCC Signal-to-Noise Ratios 

The pilot tone C/N q can be derived in a manner 
similar to that used to obtain Equation (A-34) and is given by 


(C/N ) . = K .P /[N /L + P /B ] 

O ]g g3 gs o' sg sn / s J 


K .P (L /N ) , N L >> P /B 
gj gs sg x og og sg sn s 


(A- 36) 
(A-37) 


The tone package C/N q is given by 


(C/N ). = ( 1-K .) P /N /L + P /B ] 

' o tone gj gs 7 og' sg sn 7 s J 


(A-3 8) 


- (1-K .) P (L /N ) , N /L >> P /B ( A- 3 9 ) 
gj' gs sg' og ' og ' sg sn' s v ' 

The tone power transponded by a remote terminal is a 
function of the CNR at the transponder. For a given CNR the 
ratio of signal-to-signal plus noise is 

R = S/(S + S/CNR) 


= 1/(1 + 1/CNR ( A-40) 

The actual ranging tone power transmitted by the transponder is 


P. . = RP, 

trt t 


(A-41) 
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The power received at the GCC in the i-th ranging tone '♦ 

channel is 

P = (R/m) P , L + (B /BJ P„ L „ + N (A-42) 

gr 7 ' st sg ' q' s' sn gs g l 

and the i-th ranging tone channel C/N q is 

<C/N o ) rt = (R/m) P st /[N og /L sg + P sn /B s ] (A-43) 

" <«/■»> P st (L sg/ N og> ' N o,/ L sg >> P s n / B s (A ' 44) 

where m is the number of platforms transmitting simultaneously. 
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APPENDIX B 
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IONOSPHERIC EFFECTS 


Bl. INTRODUCTION 

The net effect of the ionospheric electrons is to reduce 
the group velocity of an electromagnetic wave that passes 
through the ionosphere. This effectively increases the path- 
length relative to free space by [25] 


AR 


(1600/oj 2 ) 



meters 


(B-l) 


where gj = 2 tt f and N g is the electron density in electrons/m^ 
along the wave path. AR is a minimum when the satellite is 
nearly overhead, and it increases as the satellite elevation 
decreases because the radio wave then traverses a greater 
portion of the ionosphere . One may account for this by ex- 
pressing Equation (B-l) in the more convenient form 


AR 


4.06 x 10 


-11 


K e I/f 


2 

MHz 


(B-2) 


where K is the satellite elevation factor, f.„. is the 

frequency of the radio tone in MHz, and I = N^dh (sometimes 

referred to as the total electron content) is the integrated 

electron density (i.e., the number of electrons in a column 

2 

of the ionosphere that has a cross-section of one m ) . 

There is some difference among the satellite elevation 
factors used by various authors [26,27,28] as shown by the 
curves in Figure B-l where the biggest difference occurs for 
elevations less than 30 degrees. An average of these curves, as 
plotted in Figure B-l, will be used for the calculations in this 
appendix. 

The range error for a one-way path through the ionosphere 
at 0 degree satellite elevation is plotted as a function of 
frequency in Figure B-2 for three values of I. The curves show 
that the ionospheric range error at VHF frequencies is 100 times 
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Satellite Elevation Angle (Degrees) 

Figure B-l. Satellite Elevation Factor vs. Elevation Angle 








greater than that at L-Band frequencies, and the magnitude is * 

such that it must be examined closely. 

B2 . DIFFERENTIAL PHASE ERROR 

t 

The ionospheric phase delay on the GCC-satellite-aircraft 
link is given by 

e ion " AR up ' (2l,/X down> AR dcwn 

- -< 2 * f u /c) (K ion /f u ) - ( 2 ”V°> (K ion/ f d> 

= _(2TrK ion /c) (1/f u + 1/f d ) (B-3) 

where K. = 4.06 x 10 11 K I and the negative sign indicates 
ion e 

phase delay. 

The differential phase error between the i-th ranging 
tone and the A/R tone is given by 


= -(2 K. /c)[(l/f. + 1/f.J-l/f + 1/f ,)] (B-4) 

ion' w iu id aru ard 

Let f. = f + f + f. and f = f + f Substituting these 

10 o ar i aro o ar ^ 

into Equation (B-4) 


AO. - - ( 2 ttK . /c) [f./f 

ion ion' 1 i' iu 


V f2 iai 


(co./c) [K. /f 2 + K. /f?_] 

i' ion' iu ion' id 


(B-5) 


since f. f 

iu aru 


= f t - f./f- 
iu i' iu 


“ f i and f id f ard = £ id- 


Thus , 

the differential ranging tone phase error caused by the ionos- 
phere represents a decrease in the total path length by an 
amount given by Equation (B-2) . 




*■ 
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The mean ionospheric range error given by Equation (B-2) 
with K £ average taken from Figure B-l is plotted as a function 
of I and e in Figure B-3 for f = 135.6 MHz. The range errors 
for f = 149.22 MHz can also be obtained from Figure B-3 by 
multiplying the ordinate by 0.823. Figure B-3 was used to 
obtain the GCC-satellite-aircraf t ionosphere range errors that 

are given in Table B-l for a niqht-time minimum of 2 x 10 ^ 

2 17 

electrons/m , a "typical" day-time maximum I of 3 x 10 

2 17 

electrons/m , and a maximum day-time I of 8.6 x 10 

2 

electrons/m . Values for a two-way link are also given which 
show that the two-way ranging error can be somewhat greater 
than twice the one-way ranging error, although it generally is 
less . 

The ionospheric range/range errors in Table 1 are suf- 
ficiently large that compensation must be introduced. This 
topic will be discussed in the next subsection. 


Table B-l. Mean Ionospheric Range Error* 


Satellite 

1 = 2 

, n 16 
x 10 

1 = 3 

i n 17 
x 10 

. 1 = 8.6 

in 17 

x 10 

Elevation 
at A/C 

One-way 

Path 

Error 

Two-way 

Path 

Error 

One-way 

Path 

Error 

Two-way 

Path 

Error 

One-way 

Path 

Error 

Two-way 

Path 

Error 

e = 20° 

-0 . 16km 

-0 . 32km 

-2 . 6km 

-5 . 1km 

-6 . 5km 

-13km 

e = 45° 

-0 . 12km 

-0 . 25km 

-2km 

-4km 

-5 . 3km 

-10 . 7km 

e = 60° 

-0 . 11km 

-0 . 2 3km 

-1 . 8km 

-3 . 8km 

-4 . 9km 

-10km 


* Satellite Elevation at GCC = 25° 


B3. IONOSPHERIC RANGE ERROR COMPENSATION 
B3.1 Introduction 

One of the major problems in compensating for 
ionospheric range error is that I is subject to diurnal, 
seasonal, solar cycle, and random day-to-day variations as well 
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Ionospheric Range Error (km) 



10 16 10 17 

Integrated Electron Density (Electrons/Meter 2 ) 


Figure B-3 . One-Way Ionospheric Range Error 
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* as latitudinal and longitudinal variations. DaRosa [29] 

• 17 

cites an example where the predicted value of I was 3 x 10 

2 17 2 

electrons/m , the measured value was 8 x 10 electrons/m , 

17 2 

i and the following day it was 1.4 x 10 electrons /m . A 

summary by RCA [28] of the range of variation of I due to 
various effects is given in Table B-2 . These figures demonstrate 
that real-time measurements are required especially, as DaRosa 
states [29] , if VHF navigation precision better than 5 km is 
desired. 

Either Faraday rotation or phase-path measurements can 
be made to provide the ionospheric path corrections. The ab- 
solute accuracy of Faraday rotation measurements is about 
16 2 

1 x 10 electrons/m [30] , which corresponds to a maximum range 
uncertainty of 135 meters for the VHF frequencies that will be " 
used. Since additional equipment is required to make the 
Faraday rotation measurements, measurement of phase-path 
variations using the ranging equipment are preferable. 


Table B-2. Summary of Ionospheric Variability Factors [28] 


Factor 

Range of Variation 

Diurnal Variation 

3:1 to 13:1 

Seasonal Variation 

2:1, variable 

Solar Variation 

12:1 (11-year cycle) 

Day-to-Day Variation 

+20%, variable 

Latitudinal Variation 

2:1/10°, quite variable 

Short-term rms Variation 

1-2% 


B3.2 GCC-to-Satellite Range Compensation 

The ionospheric ranging error for the GCC-to-satellite 
link can be measured at the GCC from the phase variations of 
the 135.6-MHz satellite signal. As described previously [31], 
independent measurements of range (including ionospheric range 
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errors) and electron density could be used to establish the GCC- 
to-satellite range, and the ionospheric range error would be 
measured with respect to this reference value. This relative 
error measurement would afford somewhat hiqher accuracy [30] , 
but the absolute ranging accuracy would depend on the accuracy 
of the electron content measurement technique that was used. 
Day-time measurements of I with an absolute value of 1 percent 
have been reported [30], but if the ionospheric measurements 
are made at night when the ranging errors are smallest, the 
absolute accuracy of the measurements will also decrease. 

This basic technique can be extended to additional 
locations by placing a transponder at the desired location and 
making an additional electron density measurement during a 
calibration period. However, before determining how many 
additional remote stations are required in order to provide 
ionospheric range compensation over the satellite's area of 
coverage, the spatial and temporal variation of I will first be 
discussed. 

B3.3 Spatial and Temporal Variations of Integrated Electron 

Density 

The temporal variation of I, considering short-term, 
diurnal, and day-to-day variations, as well as diurnal variation 
measured at different locations have been studied by many in- 
vestigators. From their studies it is evident that individual 
measurements during the same period may be compared directly; 
but day-to-day, seasonal, and yearly fluctuations in I are 
great enough so that relative, rather than absolute, variations 
must be compared for different measurement periods. 

Investigators have shown [30 through 37] that the daily 
increase in I begins approximately 0500 local mean time and con- 
tinues to approximately 1200 LMT , although the time when the 
peak is reached varies from day-to-day. The initial rate of 

. 17 2 

increase is rough 0.5 x 10 electrons/m /h, which corresponds 

to a "one way" (GCC-to-satellite and return) VHF ranging error 

of 445 m/h for satellite elevations of 20 degrees. Garriott, 

et al [30] , note that a sudden change in this initial rate of 

increase occurs sometime in the morning in all seasons and all 

phases of solar activity. This change may be either an increase 

or decrease in the rate of growth of I, and the data shows a 

17 

variation in growth rate ranging from approximately 0.16 x 10 

2 17 2 

electrons/m /h to somewhat greater than 2 x 10 electrons/m /h. 

These figures correspond to increases in one-way ranging errors 







1 


of 190 m/h to 1.7 km/h. ; _ j 

^ 1 J* - I 

jL'Thte electron content begins a rapid decrease between 1 f 

• " ti -r, <- ' I ' ‘ • 

1600 an’d 1800 LMT, and reaches- a minimum by, 1900 LMT at mid- . 
latitudes "and sometime after midnight at low latitudes. At ( - 
high latitudes, a nocturnal maximum occurs during the winter at 
local midnight at 50°N, but this effect disappears beyond! 55°N 

[36] . ‘The rate of decay ^of-. I, is:: variable ,£t but i.s -,generaliy r ] . 

1: ’ y ‘ r- — / r- , * V , > t) 

somewhat ffiqher' than> the 1-growth ‘rate ./. Decay -r ate s. as high '"as | Or 
• 1 7 * ■ 2 * *■ J , ,, ’ . ‘ ‘ , . , i t: 

2.5 x TO eleetrons/m /h. ha ve^ beem noted ,owhich corresponds to-" 

• : = i j 1 

a one-way ranging error decrease of approximately 2.25 km/fr 1 . j d 

. 1 ~ ! 

/One ^relative, number/ that/ can be us'ed'ias 1 a basis; of. com- !~ 
f c, . , vv £ / ;j *; j’ , I Jj 

parison is thef'di'urrfal ratio I V /I because dLts ‘doesnexhibit H 

o r. \ i ..max' ;min • . ; ru ; 

a latitudinal. variation. < Garrioft',' et al / [ 3 0"] /" reports that it 

' | V - j 

is relatively constant at approximately 4 .during the. summer • 

months ^at''-ail mid- la t i tu'desi ( 2 0 -.4 Ocvdegreesl); and’, for /all- phases 

ro .» cj o m ©. ^ o o •-> n, - . j ^ o •/ . . . > U j 

of the /solar.; cycle , -but ;it.do"es “vary 'seasonally, and; can' teach.' 

{ t v J , » n q j ' J ^ 1 1 , I 

much higher values in winter. At high latitude's the diurnal , 
ratio ranges, from approximately .1 . 5 in summer to .3 in winter/ 

At equatorial latitudes (0-20°) the diurnal ratio exceeds 10 

[37] , and it has been noted as high as 20 ..[34]; howeyer,.as 
Table B-3 shows, the high equatorial diurnal ratios -are due more 
to low night-time values of I, for which the greatest measure- 
ment uncertainty exists, than to extremely high day-time maxima. 


V-i 


W 


rC 

h 


The model for the variation of I with latitude that is 
presented in Table B-4 was developed from References 28, 30, 
and 32 through 37. Although this variation is more properly 
expressed in terms of geomagnetic latitude, the geographic 
variation of I is required for the ranging calculations. Also, 
the data presently available does not justify a more detailed 
model. The diurnal ratios and values of I were chosen to 
reflect the upper limits of the corresponding values described 
in the references in order to obtain estimates of the upper 
limits of the differences in ionospheric ranging error between 
widely separated stations. 
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Table B-3. Equatorial Diurnal Ratios [34] 


4 


Station 

Period of 
Measurements 

I . 
mm 

Before 

Ground 

Sunrise 

, -2, 
(m ) 

I 

max 

Around 

Midday 

, -2x 

(m ) 

Ratio 

I /I . 

max min 

Zaria , 
Nigeria 

11/14/64 - 
2/13/65 

2 x 10 16 

26 x 10 16 

13 

Bangkok , 
Thailand 

June '64 - 
Sept '64 

1 x 10 16 

20 x 10 16 

20 

Huancayo , 
Peru 

8/15/61 - 

10/31/61 

5 x 10 16 

40 x 10 16 

8 

Huancayo , 
Peru 

1/1/67 - 

3/31/67 

8 x 10 16 

70 x 10 16 

8.7 


Table B-4. Model for Latitude Variation of I 


Latitude 

I (10 17 
I 

max 

o 

electrons/nr) 

I . 
mm 

Ratio 

I /I . 

max min 

0° 

8 

0.67 

12 

20 °N 

5 

0.6 

8 

30 °N 

3 

0.6 

5 

40 °N 

2 

0.5 

4 

£ 

O 

O 

in 

1.6 

.0.4 

4 

60 °N 

1 

0.33 

3 


The longitude variation of I will be based upon the 

earth's rate of rotation of 15 deg/h. Thus, for two stations 

separated by 90 degrees (6 hours in time) , the value of the 

integrated electron density at one station could be I max » and 

the value at the other station could be I • . The values of I 

mm 
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from Table B-4 that will be used to calculate the ionospheric 
ranging error to geographically separated stations are given in 
Table B-5. 

Table B-5. I vs Longitude for Ionospheric Error Model 


Longitude of I 
3 max 


Station 

Longitude 

8 °W 


53°W 


9 8°W 


(Model A) 

i (Model B) ' 

i 

(Model C) 

(Model D) 

(Model E) 

8 °W 

I 

i 

0.5 I 

0.5 I 

j I ■ 

! 

I 

max 

max 

max 

max 

mm 

53°W 

0.5 1 

j I 

I 

I 

1 0.5 1 


max 

max 

max 

max 

max 

9 8 °W 

I . 

| 0.5 I 

0.5 I 

I 

' I 


mxn 

L 

j max 

1 

max 

max 

max 

1 

1 


The values of I given in Table B-5 were used to cal- 
culate the ionospheric range error from the GCC located at 33°N, 
9 8°W to stations at 0°, 20°N, 30°N, 40°N, 50°N, and 60 °N at 
longitudes of 8°W, 53°W, and 98°W. The satellite elevations 
used were obtained from an assumed position of 50°W for the 
satellite . 

The ranging errors that were calculated from the model 
in Table B-5 are given in Tables B-6 through B-10, in which 
AR^ is the one-way (GCC-to-satellite-to-remote station) iono- 
spheric range error, AR 2 is the two-way (one-way plus return) 
ionospheric range error, and AG is the GCC-to-satellite-to-GCC 
ionospheric range error. The one-way error ranges in value from 
-0.4 km to -4.1 km and has a mean value of -1.7 km and a 
standard deviation of 0.8 km. The two-way error ranges in value 
from -0.8 to -7.4 km and has a mean value of -3.35 km and a 
standard deviation of 1.6 km. 

,*> One possible compensation method is to subtract AG from 

AR^, since this may be readily measured at the GCC. The re- 
sultant error decreases to a range of 0 to 2.8 km and has a 
mean of -0.2 km and a standard deviation of about 0.8 km. 
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Table B-6. Ionospheric Ranging Errors Predicted by Model 
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Table B-7. Ionospheric Ranging Errors Predicted by Model 
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Table B-8. Ionospheric Ranging Errors Predicted by Model 
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Table B-9. Ionospheric Ranging Errors Predicted by Model 
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Table B-10. Ionospheric Ranging Errors Predicted by Model 
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Although this procedure does decrease the error and its mean 
value , the standard deviation remains the same so the uncertainty- 
in the error value is unchanged. 

Another possible compensation method is to measure the 
two-way ionospheric path error to a given location and subtract 
half this value from the one-way error. This procedure can re- 
duce the magnitude of the maximum ionospheric error to less than 
0.4 km. The mean error is -0.03 km, and the standard deviation 
is 0.12 km so a substantial improvement has been made. However, 
this improvement comes at the expense of the large number of 
remote stations required to obtain accurate ionospheric range 
error estimates at points within the satellite's region of 
coverage . 

The maximum value, mean value, and standard deviation 
of AR^ and AR 2 is tabulated in Table B-ll for various regions 
to see if regional coverage with a single compensation value 
appears feasible. The figures show that the higher latitudes 
might provide the largest reqion in which a single ionospheric 
range error estimate might apply. The mean value of 0.5AR 2 
at 5 0 °N in Table B-6 was subtracted from the AR^ values for 
40 °N , 5 0 °N and 60°N in Table B-6, and a similar procedure was 
followed in Tables B-7 through B-10. The maximum resultant 
error was -0.6 km, and except for the GCC longitude the 
resultant error was substantially less than that obtained by 
subtracting AG from AR^. 

Since the mean ionospheric range errors can be several 
kilometers, a compensation method is necessary. The most 
promising scheme is to obtain real-time measurements from remote 
stations and retard the phase of the ranging tones by the amount 
derived from the ionospheric range error measurements. A number 
of different phase delays will be required because a single 
value cannot serve all regions. The number of different phase 
delays that will be required depends upon the spatial variation 
of I at a given moment, the satellite elevation at the receiving 
station, and the desired ranging accuracy. One or two compensa- 
tion values might suffice for the period during which I . 

applies for the entire satellite region of coverage, but several 

values will be necessary otherwise. 

Since the ionosphere changes relatively slowly the phase 
compensations need not be updated continuously, but only every 
few minutes; the range error estimation system capacity need not 
be very large. 



Table B-ll. Ionospheric Range Errors for Various Regions 
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W B3.4 Ionospheric Scintillation 

Irregularities in the ionospheric structure produce 
amplitude fluctuations, or scintillations. The fluctuations at 
3 mid-latitudes are not as large as those at equatorial and high 

latitudes [10] nor do they occur as frequently. To be con- 
sistent with the discussion of multipath effects presented in 
Appendix C, it is convenient to consider the scintillation 
fluctuation effect in terms ov the fade depth of the received 
signal, where 

Fade Depth = 20 log 1Q (1 - m^) 


and 


i 


\ 


m . = Scintillation Index 
sx 


(P 

max 


P . ) /(P + P . ) 

min ' max mm 


Aarons and Whitney [11] investigated scintillation 
fading distribution at mid-latitudes. They showed that 90 per- 
cent of the time the fade depth was less than 1.3 dB (m . = 0.14) , 

b -L 

and the 99 percentile value was -9.4 dB (m . = 0.67) . Both the 

S 1 

GCC-to-satellite and satellite-to-remote station links can ex- 
perience ionospheric scintillations of these magnitudes, simul- 
taneously, but it is possible that only one of the links might 
experience scientillation fading or that the fade values on the 
two links will differ because of the irregular boundary of the 
region of ionospheric irregularity [10], 

The duration of the mid-latitude scintillations ranged 
from a few minutes to a few hours with scintillation bursts of 
3 to 14 minutes duration occurring near the daytime ionization 
peaks [11]. The distribution of the fading periods is 
summarized in Table B-12 . 


Table B-12. Mid-Latitude Scintillation Fading Period [11] 
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The net effect of ionospheric scintillation is to pro- 
duce an amplitude fading and a corresponding phase shift similar 
to that produced by specular multipath interference, as discussed 
in Appendix C. The phase shift in this case, however, will be a 
function of the change in I. 



APPENDIX C 


MULTIPATH EFFECTS 


Cl INTRODUCTION 

The multipath effects that will be considered for the 
satellite-to-aircraft links arise from interference between the 
direct path signal and an earth-reflected signal as shown in 
Figure C-l. As a result of this interference, the amplitude of 
the received signal will vary so that the effective link 
signal-to-noise ratio will be reduced at times, and the signal 
phase will vary so that a range error will be introduced. 


Since the ranging information is transmitted on a set 
of spaced tones, which may be considered as separate unmodulated 
carriers, the multipath effects may be calculated by using 
the vector model in Figure C-2. Here the direct signal is 
represented by a unit vector, a specular (phase coherent) 
reflected signal is represented by a vector with an amplitude 
m (0 < m < 1) and a time-varying (caused by aircraft motion) 

b b 

phase 4> s , and a diffuse (random phase) reflected component is 
represented by a vector with an rms amplitude m d (0 < m^ < 1). 
Both types of reflected signals will usually be present, 
although one may be sufficiently stronger than the other so that 
it will essentially determine the amplitude and phase variations 
of the received signal. 

C2 REFLECTED SIGNAL PARAMETERS 

The amplitude and phase of the reflected signal are 
functions of several factors such as signal polarization, 
reflecting surface characteristics, the antenna pattern, and 
satellite elevation angle. Since these factors will vary over 
any given flight path, a statistical description of the reflected 
signal is required. Bergemann and Kucera [7, 9, 23] have 
reported signal fading distributions and fading spectra obtained 
from 135.6 MHz-transmissions from ATS-I and ATS-III to aircraft 
equipped with a Dome and Margolin SATCOM antenna. 

The amplitude fading distribution obtained by Bergemann 
and Kucera were found to differ mainly with satellite elevation 
angle. Although some difference between overland and overwater 
fading was noted [7], the distribution obtained as a function 






Figure C-2. Vector Multipath Model 


of satellite elevation angle generally applied for all types 
of terrain. Bergemann and Kucera's data, as well as that of 
other investigators [38, 39] indicate that the fading is 
primarily due to specular multipath interference for satellite 
elevation angles less than 30 degrees, while diffuse fading 
is more common for higher satellite elevation angles. 

Kucera [9] and Jordan [39] have found the spectral 
components of specular fading to be generally less than 0.3 Hz, 
but their measurements differed on the bandwidth of the 
diffuse spectrum. Kucera noted that except for cases of very 
strong multipath interference or satellite spin modulation, the 
fading spectrum did not extend beyond 2.5 to 3 Hz while Jordon 
reported spectral components extending somewhat beyond 10 Hz. 
Kucera 1 s data also showed that the diffuse signal was composed 
of relatively few components whose spectral envelope' has 
approximately a gaussian shape; however, the sum of four to 
six randomly phase sine waves will closely approach the 
Rayleigh distribution [40] that is commonly used to describe 
the diffuse signal, and a narrow-band filter will give a 
uniform shape to the diffuse spectrum so that the Nakagami-Rice 
model [41, 42, 43] can be used to describe the amplitude and 
phase distributions of the received signal when diffuse 
multipath interference is present. 

Bergemann and Kucera noted no evidence of selective 
fading on transmissions consisting of 25 tones spaced one kHz 
apart. The fading period was approximately the same for each 
tone, but the tone fades did not occur simultaneously in time 
because of the frequency difference between each tone. 

The reflected signal parameters that will be used in sub- 
sequent analyses have been derived from Bergemann and Kucero's 
data [44] and are given in Table C-l, in which M, = 20, m.,. 

d log 10 d 


Table C-l. Reflected Signal Parameters 
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C3 


DEMODULATION PROCESS 


In the modulation process the ranging tones are trans- 
lated from baseband to radio frequency with respect to a 
reference (A/R) tone. The demodulation process shown in Figure 
C-3 performs the inverse function by filtering the A/R tone 
from the received tones and using it to coherently translate 
the ranging tones back to baseband for phase comparison with 
the transmitted ranging tones. 


From Figures C-2 and C-3 the input into the i-th phase 
comparator is 


whe re 


e.(t) = A. A /2P . cos [to. t + A4> . + A0.] 

1 1 ar 1 L l l 


= [(1 + m s coscj> s + m dc ) 

2 V2 

+ (m s sin<j) s + m ds ) ] 


A4> 


i = -2TT(1/A i - 1/X ar )d = 0) i i 


A0^ = arctan 


m s sin »si + m dsi 
1 + m s cos4, simacl| 


- arctan ' 


m s sin *sar + m dsar 


1 + m cosd) + m, 

s T sar dcsar 


(C-l) 


(C-2) 

(C-3) 


(C-4) 


and P^ is the received power in the direct path signal, x d 
is the direct path time delay, and the in-phase and quadrature 
diffuse components have gaussian distributions with 


and 


m 


dc 



m 


ds 


zero 


2 

m, 

dear 


2 2 

m. . = m, 
dsi dsar 



(C-5) 


From Equation (C-l) it can be seen that the filtered A/R tone 

input to the mixer should be limited (A = 1) so additional 

sr 

multipath amplitude modulation of the ranging tone is not 
introduced. Considering only the specular multipath terms 
in Equation (C-4) , it is seen that I s minimized because 
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Figure C-3. Ranging Tone Detection Process 
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|; however, when considering only the diffuse com- 
ponents in Equation (C-4) , it is seen that the variance of the 
phase noise increases rather than decreases. This implies that 
if a phase-locked- loop (PLL) filter such as that in Figure C-4 
is used to filter the A/R tone, the multipath phase error will 
be minimized if the PLL filter passes the slower specular phase 
variations and filters the diffuse phase variations. 


Two types of phase detectors will be considered. The 
first is the linear phase detector shown in Figure C-5 (a) , and 
the second is the coherent phase detector shown in Figure 
C-5 (b) . The output of the linear phase detector will be 


e 


lin 


A4> . + A 0 . 
l 1 


(C-6) 


while the output of the coherent detector will be 

e = sin (A<£ . + A6.) (C-7) 

c 1 i 

In each case, the phase detector is preceded by a limiter to 
eliminate the amplitude fluctuation effects, and the output is 
filtered to reduce the multipath phase error. 

C4 MULTIPATH RANGING ERRORS 

The specular multipath phase error, although diminished 
because of the phase difference, presents the most difficult 
problem because it fluctuates at a rate (< 0.3 Hz) which is 
far less than the output filter bandwidth. By utilizing various 
trigonometric identities, the specular multipath phase error in 
Equation (C-4) can be expressed by [44] 


AO 


(C-8) 


arctan 


-{m /2 sin(7rAr/T.) 

S 1 


[cos ( <p -7tAt/T. ) +4m cos (vAt/T. ) ] J 

SclJT X S 1 


l-m^+ (m /2) cos (ttAt/T . ) [cos ( <j> -ttAt/T. ) + 4m cos (ttAt/T . ) ] 

S S IS a. i IS 1 . 


where Ax is the differential delay between the direct and re- 
^ fleeted signals and T^ is the period of the i-th ranging tone. 

Since <p will vary over a 2 tt range due to aircraft motion, 
the magnitude of the maximum value of A0 s and the rms value 
were calculated for a 2tt range of values for <j> and 
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Figure C-4. A/R Tone Filter 
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plotted as functions of Ax/T, in Figures C-6 and C-7. The 
curves are repetitive over unit ranges (e.g., 0 to 1, 1 to 2, 

20 to 21)of Ax/T^. 

The reduction in the magnitude of the maximum specular 
multipath phase error allowed by the differential process can 
be noted by comparing Figures C-6 and C-8. For M s = -14 dB 

(where M„ = 20. m ) the maximum error for a single tone 
log 10 s 

is 11.5 degrees while the maximum differential error 
is 6.8 degrees. The maximum value of the corresponding rms 
differential error is 5.2 degrees. 

For satellite elevation angles between 10 and 30 degrees, 
the differential multipath time delay will vary from 10 ys 
to 55 ys. The rms and maximum differential multipath phase 
errors for 426-Hz, 707-Hz, and 941-Hz ranging tones were 
obtained from figures C-6 and C-7 for an rms reflected signal 
amplitude of = -14 dB and a "worst-case" (one percent) value 
of M g = -5 dB. The corresponding range errors were determined 
from Figure C-9, and the results are given in Table C-2 for 
each ranging tone frequency (f^). For subsequent calculations, 
a 1-km rms multipath ranging error will be assumed together with 
a "worst-case" value of 4 km. 


Table C-2. Specular Multipath Ranging Errors 


f . 
1 


Ax/T. 

0 (rms ) 

R (rms) 

| © s | (max) 

R s (max) 





Ax = lOys, 

M = -14 dB 
s 


426 

Hz 

.00426 

<0.1° 

<190m 

0.12° 

2 35m 

707 

Hz 

. 00707 

0.14° 

162m 

0.2° 

2 35m 

941 

Hz 

.00941 

0.162° 

145m 

0.27° 

240m 





Ax = 55ys, 

M = -14 dB 
s 


426 

Hz 

.0234 

0.39° 

760m 

0.65° 

1 . 2 7km 

707 

Hz 

.0384 

0.63° 

730m 

1.1° 

1 . 2 8km 

941 

Hz 

.0518 

0.82° 

720m 

1.45° 

1 . 2 8km 



Magnitude of Maximum Differential Specular Multipath Phase Error (Degrees) 



At/T 1 

Figure C-6 . Magnitude of Maximum Differential Specular 
Multipath Phase Error 
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RMS Differential Specular Multipath Phase Error (Degrees) 



Figure C-7. RMS Differential Specular Multipath Phase Error 
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Magnitude of Multipath Phase Error (Degrees) 


Figure C-9. Multipath Ranging Error 
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Table C-2. Specular Multipath Ranging Errors (Cont.) 


f . 
1 


Ax/T i 

0 g (rms) 

R (rms) 
s 

| 6 s 1 (max) 

R (max) 
s 




At = 10ys, 

M = -5 dB 
s 


426 

Hz 

.00426 

0.35° 

6 80m 

0.43° 

830m 

707 

Hz 

.00707 

0.58° 

6 80m 

0.7° 

810m 

941 

Hz 

. 00941 

0.78° 

690m 

0.95° 

8 40m 




At = 55ys, 

M = -5 dB 
s 


426 

Hz 

.0234 

1.95° 

3 . 8km 

2.4° 

4 . 7km 

707 

Hz 

. 0389 

3.1° 

3.6km 

3.9° 

4 . 6km 

941 

Hz 

.0518 

4.2° 

3 . 7km 

5.2° 

4 . 6km 


For satellite elevation angles of 30 to 50 degrees and 
50 to 70 degrees the reflected siqnal was assumed to be diffuse 
with relative rms amplitudes of -14 dB and -17 dB, respectively. 
The output signal-to-noise ratio (SNR) is given by 

v = - M d - 3 dB + I d0 (C-9 ) 

where the 3 dB factor accounts for the A/R tone diffuse noise 
and I^g is an improvement factor due to filtering. If the A/R 

tone filter bandwidth is one Hz and the phase detector bandwidth 
is also one Hz, then the SNR will improve somewhat due to 
reduction in the diffuse noise at the filter outputs. If a 
gaussian shape is assumed for the diffuse noise spectrum, 
then I^g = 4.4 dB if the diffuse bandwidth (94 percent of the 

power ) is 3 Hz, and I^g = 8.2 dB if the diffuse bandwidth (B n ) 

is 5 Hz [44]. When the SNR is determined from Equation (C-9), 

the rms diffuse phase error (0^) may be obtained from Figure 

C-10, and the ranging error (R^) can be determined from Figure 

C-9, The results are tabulated in Table C-3 for each ranging 

tone frequency (f^). 
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Figure C-10. RMS Phase Error For NaJcag ami -Rice Distribution 
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Table C-3. Diffuse Multipath Ranging Errors 



| 

I 

B 

n 

= 3 H 

z 


B = 5 H 
n z 


f . 
1 


0 d 

R d 


0 d 

R d 




M d = -14 

dB 



426 

Hz ' 

7° 

13.6km 


4.4° 

8 . 7km 

707 

Hz 

7° 

8 . 1km 


4.4° 

5.1km 

941 

Hz ■ 

7° 

6 . 2km 


4.4° 

3.9km 




t" 
i — i 

i 

ii 

T 

dB 



426 

Hz ' 

4.9° 

9 . 5km 


3.3° 

6 . 5km 

707 

Hz 1 

4.9° 

5 . 7km 

1 

3.3° 

3 . 9km 

941 

Hz I 

4.9° 

5 . 3km 


3.3° 

2 . 95km 


From Table C-3 it is seen that the diffuse multipath 
ranging error will be less than the specular ranging error. 

The major factor that accounts for this somewhat unusual 
situation is the relatively narrow diffuse noise bandwidth 
indicated by the Collins spectral data [7, 8, 9], Further 
multipath spectral measurements to verify this are needed. 

Since the magnitude of the reflected signal depends 
directly upon that of the transmitted signal, any reduction 
of the magnitude of the reflected signal cannot occur 
through variation of the transmitter power, but rather through 
more attenuation of the reflected signal by the aircraft antenna. 
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